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Abstract 



This thesis presents the analysis of the clustering of galaxies in the 6dF Galaxy Survey 
(6dFGS). At large separation scales the baryon acoustic oscillation (BAO) signal is de- 
tected which allows a measurement of the distance ratio, r s (zd)/Dv{z e g) = 0.336 ± 0.015 
(4.5% precision), where r s (zd) is the sound horizon at the drag epoch Zd and Dv{z e s) is 
the absolute distance to the effective redshift of the survey, given by z e s = 0.106. The 
low effective redshift of 6dFGS makes it a competitive and independent alternative to 
Cepheids and low- z supernovae in constraining the Hubble constant. The value of the 
Hubble constant reported in this work is Ho = 67.0 ± 3.2 km s -1 Mpc" 1 (4.8% precision) 
which depends only on the WMAP-7 calibration of the sound horizon and on the galaxy 
clustering in 6dFGS. Compared to earlier BAO studies at higher redshift, this analysis is 
less dependent on other cosmological parameters. This thesis also includes forecasts for 
the proposed TAIPAN all-southern-sky optical galaxy survey and the radio WALLABY 
survey. TAIPAN has the potential to constrain the Hubble constant with 3% precision 
using the BAO technique. 

Modelling the 2D galaxy correlation function of 6dFGS, £(r p , it), allows a measure of the 
parameter combination f(z e s)as(z e s) = 0.423 ± 0.055, where / ~ &m{z) is the growth 
rate of cosmic structure and a$ is the r.m.s. of matter fluctuations in 8/i _1 Mpc spheres. 
The effective redshift of this analysis is z c g = 0.067. Such a measurement allows to 
test the relationship between matter and gravity on cosmic scales by constraining the 
growth index of density fluctuations, 7. The 6dFGS measurement of fas combined with 
WMAP-7, results in 7 = 0.547±0.088, consistent with the prediction of General Relativity 
(7GR ~ 0.55). Because of the low effective redshift of 6dFGS this measurement of the 
growth rate is independent of the fiducial cosmological model (Alcock-Paczynski effect). 
Using a Fisher matrix analysis it can be predicted, that the WALLABY survey will be 
able to measure fas with a precision of 4 - 10%, depending on the modelling of non-linear 
structure formation. This is comparable to the predicted precision for the best redshift 
bins of the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that low- 
redshift surveys have a significant role to play in future tests of dark energy and modified 
gravity. 

The last chapter of this thesis studies the stellar-mass dependence of galaxy clustering in 
the 6dF Galaxy Survey. The near-infrared selection of 6dFGS allows more reliable stellar 
mass estimates compared to optical bands used in other galaxy surveys. Using the Halo 
Occupation Distribution (HOD) model, this analysis investigates the trend of dark matter 



halo mass and satellite fraction with stellar mass by measuring the projected correlation 
function, w p (r p ). This is the first rigorous study of halo occupation as a function of stellar 
mass at low redshift using galaxy clustering. The findings of this analysis are, that the 
typical halo mass (Mi) as well as the satellite power law index (a) increase with stellar 
mass. This indicates, (1) that galaxies with higher stellar mass sit in more massive dark 
matter halos and (2) that these more massive dark matter halos accumulate satellites faster 
with growing mass compared to halos occupied by low stellar mass galaxies. Furthermore 
there seems to be a relation between Mi and the minimum dark matter halo mass (M m ; n ) 
of Mi ~ 22M m i n , in agreement with similar findings for SDSS galaxies, but explored 
for the first time as a function of stellar mass. The satellite fraction of 6dFGS galaxies 
declines with increasing stellar mass from 21% at M ste n ar = 2.6 x 10 10 /i~ 2 M to 12% at 
Mgtciiar = 5.4 x 10 10 /i~ 2 Mq indicating that high stellar mass galaxies are more likely to 
be central galaxies. Finally the 6dFGS results are compared to two different semi-analytic 
models derived from the Millennium Simulation, finding some disagreement. The 6dFGS 
constraints on the satellite fraction as a function of stellar mass can be used to place new 
constraints on semi-analytic models in the future, particularly the behaviour of luminous 
red satellites. 
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Chapter 1 
Introduction 



With the introduction of General Relativity in 1915 by Albert Einstein Einstein 



1915 , it became possible for the first time to describe the gravitational Universe 



with one set of equations. Nevertheless it was not until later in the century that 
the modern model of cosmology emerged. Although complemented by many the- 



oretical achievements like inflation Guth, 1981 and the explanation of the origin 



of heavy elements by Big Bang nucleosynthesis Alpher, Bethe k Gamow , 1948 , it 
were largely observational breakthroughs which led to what we now call the stan- 
dard model of cosmology. The most important of these were the discovery of the 



Hubble expansion Hubble, 1929 , the discovery of dark matter Zwicky, 1937 Kahn 



k Woltjer, 1959 Freeman et al. , 1970; Rubin k Ford, 1970 , the discovery of the 



Cosmic Microwave Background (CMB) jPenzias k Wilson 1965 , the discovery of 
fluctuations in the CMB Smoot et al. 11992 Bennett et al. 1996 and the discovery 



of the accelerating expansion of the Universe Perlmutter et al. , 1998 Riess et al. 



1998 as well as the transition to precision cosmology Spergel et al. , 2003 



In this chapter we give a brief introduction to the theoretical basis of modern 
cosmology. We also discuss the tools of observational cosmology which are used in 
this work to test the current standard model, namely baryon acoustic oscillations 



(BAO) and redshift space distortions. For a more detailed discussion we refer to Kolb 



k Turner| [1994] ; |Peacock| [1999] ; |Dodelson| [2003]; |Carroll| [2004a]; |Lyth k Liddle 
20T0]. 
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Supernova Cosmology Project 
Suzuki, et al., Ap.J. (2011) 



Union2. 1 SN la 
Compilation 




Figure 1.1: Cosmological constraints indicating dark energy (Qa) comprising 73% 
of the total energy density and (dark) matter (fl m ) making up 27%. (credit: [Suzuki 
et al.||201ip 



According to current best estimates (see e.g. Komatsu et al. 2011 Suzuki et al. 
2011 Sanchez et al. 2012 and Figure [TTT] for the latest results), approximately 73% 
of our Universe is made up of dark energy, about 27% is in the form of matter, but 
only 20% of this matter fraction is in the form of baryonic matter, while the rest is so 
called cold dark matter (CDM). Recent constraints on dark energy as a cosmological 
constant (Qa) and matter (Q m = cold dark matter + baryonic matter) are shown 



in Figure |1.1| We will discuss the cosmological probes used in this diagram within 
this introduction, especially the CMB and BAO techniques. 

It will be crucial for the understanding of the concepts used in this work to 
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CHAPTER 1. INTRODUCTION 



be familiar with the basic ideas of General Relativity and the expanding Universe. 



Hence we start with a brief introduction to General Relativity in section fET including 
a discussion of the solutions of the Einstein equations that are believed to describe 
our Universe. In section 11.21 we will introduce the statistical tools used in this 
analysis, namely the correlation function and power spectrum. In section 1.3 we 



introduce inflation followed by a discussion of the fundamental equations for the 



evolution of the initial matter fluctuations. In section 1.4 we discuss the impact of 
peculiar velocities on the distribution of galaxies, and how we can detect this signal 



with galaxy redshift surveys. In section 1.5 we introduce the halo model, which 



connects the matter and galaxy clustering. In section |1.6| we discuss techniques of 
modern observational cosmology, focusing on Baryon Acoustic Oscillations and their 
signature in the CMB and the local galaxy distribution. In section 1.7\ we discuss 
alternative models for dark energy, some of which will be tested in this analysis. In 
section 11.81 we introduce the 6dFGS dataset which is used in this work. Section 11.91 
will give a overview and a motivation for this thesis and introduce the following 
three chapters. 



1.1 General Relativity 

Since gravity is the dominant force on cosmic scales, the evolution of the Universe 
is described by the Einstein field equations, which relate the mass-energy content 
of the Universe to the geometry of space-time, including a geometrical cosmological 
constant (we use units where the speed of light is c = 1): 

G^ v = R^u - -jQiivR = 8-kGT^ + g^A, (1.1) 

where is the Einstein tensor describing the geometry of the Universe, R^ v is the 
Ricci tensor, R is the Ricci scalar, g^ u is the metric tensor, A is the cosmological 
constant and G is Newton's gravitational constant. The matter distribution is char- 
acterised by the energy-momentum tensor T Ml/ , which in the case of a perfect fluicj^ 
leads to T 00 = p which is the energy density and Ta = pi which are the pressure 

X A perfect fluid is one that can be completely specified by two quantities, the rest-frame energy 
density p, and an isotropic rest-frame pressure p that specifies the pressure in every direction. 
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1.1. GENERAL RELATIVITY 



components for i = 1,2,3. The indices /i = 0,1,2,3 correspond to the space-time 
coordinates x M = x 1 , x 2 , x 3 ) = (ct,x,y,z). 

In order to find solutions for these equations we need several simplifying assump- 
tions. The most important of these is the Cosmological Principle, which states that 
on the largest scales our Universe is homogeneous and isotropic, i.e. there are no spe- 
cial locations or directions. Modern cosmology provides very strong support for this 



assumption from the observation of the almost constant CMB temperature Smoot 



et al. , 1991 Bennett et al. , 1996 as well as large scale structure measurements at 



low redshift Hogg et al. , 2005; Scrimgeour et al. 2012 . Using these assumptions, 



the most general space-time interval can be written as: 



ds 2 = dt 2 - a 2 {t)[d X 2 + S 2 k { X ){d6 2 + sin 2 6d<f> 2 )}. 



;i.2) 



This is the Friedmann-Lemaitre- Robertson- Walker (FLRW) metric Friedman, 1922 



Lemaitre 1931 Robertson 1935; Walker 1937 with t being the proper time mea- 
sured by a co-moving observer, a(t) is the scale factor, 9 and <fi are the spherical 
coordinates, and k is known as the curvature parameter. Furthermore 



s k (x) 



sinx if k — +1 (closed) 
X if k = (flat) 

sinhx if k — — 1 (open), 



1.3) 



with x being a dimensionless radial coordinate related to physical coordinate dis- 
tance by r = ax- We call a Universe open if k = — 1, closed if k — +1 and flat 
(Euclidean) if k = 0. As we will see later, k = seems to be the case favoured by 
current observational data. Using the scale factor a(t) we can define the cosmological 
redshift as z = ao/a(t) — 1, where ao is the scale factor today. 

We can now move on and solve the Einstein field equations to obtain a set of 
relations which describe the dynamics of the scale factor a(t). Einstein's equations 
relate the evolution of this scale factor, a(t), to the pressure and energy of the matter 
in the Universe 

H 2 = -[ - ) = — p-^ + i, (1.4) 
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which is derived from the 00-component of Einstein's field equations, and 

a AttG , „ . A ,„ . 

- = -— (p + 3p) + - (1.5) 

which is derived from the trace of Einstein's field equations. The energy density p 
and pressure p correspond to all particle species, relativistic or not. The equations 
above also include the definition of the Hubble parameter 

Hit) ee a - (1.6) 
a 

where the dot denotes a derivative with respect to time t. 

Whether the curvature of the Universe is closed (k = +1), flat (k = 0) or open 
(k = —1) depends on the energy density of the components it contains. The density 
required to make the Universe flat is called the critical density p c and is found by 
using the first of the Friedmann equations and setting A and the normalised spatial 
curvature, k, to zero: 

p c = = 2.78/r 1 lO 11 M / i 3 Mpc- 3 , (1.7) 

where h is defined as 

~ 100km/s/Mpc' 1 ' ' 

The density parameter is then defined as: 



(L9) 



The evolution of the density as the Universe expands can be studied by the adiabatic 
equation 

d(aW) da? 

dt F dt v 1 

which can be derived from the Friedmann equations above. Characterising each 
component by its equation of state 

w=- (1.11) 
P 
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and making the ansatz p oc a n , one obtains 



p(a) = p a- 3{w+1) oc { 



a" 4 if 


it) 


1 

3 


radiation & relativistic species 


gT 3 if 


w 


= 


dust (or pressure-less matter) 


a' 2 if 


w 


1 

3 


curvature 


const. if 


w 


= -1 


A & inflation, 



;i.i2) 



where po denotes the density value at the present time. 

Using the density parameter and their evolution discussed above, and substi- 



tuting them into the Friedmann equation (eq. 1.4) we can write 

H(a) 2 



tt2 ud 7 
H 



fi 7 a _4 + fi M a~ +£4a~ 2 + fi A , (1-13) 



where we treat curvature also as a density defined by 



* = -i^- (L14) 



Eq. |1 .13 is the version of the Friedmann equation which is most useful for the 



study of cosmology, because the different £Vs are directly related to observational 
quantities. We will use this equation frequently in this thesis. 

Since radiation density (including neutrinos, if they are relativistic) follows p ~ 
a -4 and matter density p ~ a~ 3 , we can infer that there must have been a time 
when the Universe was dominated by radiation and the transition from radiation 
domination to matter domination is given by 



P 



r 



a eq = — ~ 3 x 10" 4 . (1.15) 

Pm 

For a > a eq the expansion was matter-dominated. Only very recently (a > 0.5) dark 
energy became the dominating energy component of the Universe and began to play 
a major role in the expansion of the Universe. In fact, since dark energy does not 
appear to depend on time, while all other energy components decrease with time, 
dark energy will only become more dominant in the future. 
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CHAPTER 1. INTRODUCTION 



1.1.1 Cosmic distances 



Within this formalism it is easy to understand how the Hubble law follows directly 
from the cosmological principle and the FLRW-metric. Following our nomenclature, 
the proper distance is defined by r = a(t)x, whereas its velocity is given by v = 
dr/dt = ax = ar/a. Putting this in the definition of the Hubble parameter we find 



H - 
-no — ~ 

a 



o 



vx, (1-16) 



which is Hubble's law and H is the value of the Hubble parameter today 

The co-moving distance Dc of an object that emits light at conformal time t e 
and is observed today (t ) is 

„ f to cdt f 1 cda f z cdz 

Dc= —= -^773= / 7773> (1-17) 



l ac a*H(a) Jo H(z) ' 
where we used dt = da/(aH(a)). Using the combined density parameter 

E(z) = [fi 7 a" 4 + n m a- 3 + n k a~ 2 + ft A ] 1/2 , (1-18) 

the distance equation can be re-written as 

f z cdz c f z dz , 

*=Lwr*Lwr (L19) 

The co-moving distance Dc is the most common distance used in this analysis. 
However we also employ the transverse co-moving distance: 

D M (z) = -^r^=S k (y^h) , (1.20) 



h q y/rr k 



where we used eq. 1.3, x — H Dc(z)/c and k = +1, k = and k — — 1 for the 
cases of positive curvature (f2 fc > 0), no curvature (f^ = 0) and negative curvature 
(Q k < 0), respectively. From this we can directly derive the angular diameter 
distance: 

D A (z) = 2^£> (1.21) 
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and the luminosity distance is denned as 

D L (z) = (l + z)D M (z). (1.22) 

The angular diameter distance does not increase indefinitely as z — > oo, instead it 
turns over at z ~ 1. Consequently, objects with higher redshift appear larger in 
angular size. 



1.2 Clustering statistics 



The Universe we live in was seeded by quantum fluctuations, which were pushed 
to cosmic scales during inflation and grew by the subsequent evolution to form the 



highly nonlinear structures we observe today Kolb & Turner, 1994; Peacock, 1999 



As a consequence of this quantum mechanical origin, the structures are stochastic 
with random initial conditions. Due to this fact we cannot hope to develop a theory 
that exactly reproduces the Universe we observe today. Rather, we should consider 
our Universe as one representation of an ensemble of possible Universes. Therefore, 
we need to introduce statistical quantities, which can be used to compare theoretical 
predictions with the observed data. 



1.2.1 Two-point correlation function and power spectrum 

We can define a dimensionless over-density or density contrast 

*<r) = (1.23) 

p 

where p(r) is the density at a specific position r and p is the average density. The 
Two point correlation function is now defined as 

£(r) = (6(r + x)8(r)). (1.24) 

The correlation function determines the probability of finding two objects at sepa- 
ration r in excess of the probability one would expect for a random distribution. 
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It will prove convenient to build up the actual density field from a superposition 
of modes that describe the behaviour on a certain scale. We therefore consider a 
finite box of volume V with periodic boundary conditions and expand the density 
contrast in terms of its Fourier components, allowing one to write the correlation 
function as 

= ^3 y d 3 k\5(k)\ 2 exp[-iAr], (1.25) 

with the wavenumber k. The ensemble average of the squared amplitude of the 
Fourier component of a field 

P(k) = (\5(k)f), (1.26) 

is called a power spectrum. Using spherical symmetry the relation between the 
correlation function and the power spectrum reduces to a Hankel transform 

= (2^3 / d 3 kP(k)k 2 j (kr), (1.27) 

where jo(kr) = sin(kr)/kr refers to the spherical Bessel function of order 0. 



1.2.2 Filtering of the density field and the linear bias scheme 

As we are not only interested in the local properties of perturbations, but also in 
averages over a certain volume, we can convolve the density field with a filter Wr 
of scale R (often called window function). This convolution in real space translates 
into a simple multiplication in Fourier space. The variance of the smoothed density 
field is given by 

where we use the Fourier transform of the top hat filter with radius R 

W R (k) = [sm(kR) - cos(kR)} . (1.29) 

The scale R of the filter is related to a typical mass by the relation 

Air 

M = —R 3 p. (1.30) 
o 
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1.3. INITIAL CONDITIONS AND THE EVOLUTION OF THE PRIMORDIAL 
POWER SPECTRUM 



Historically R = 8/i -1 Mpc has been chosen as the smoothing scale at which to 
quantify the variance of the density field and is denoted by as- This parameter 
normalises the matter power spectrum (P(k) oc erf) and we could expect to derive 
this parameter from the amplitude of the measured power spectrum. However a 
galaxy redshift survey actually measures the galaxy power spectrum P g (k), which 
is not necessarily identical to the matter power spectrum P m (k). On linear scales 
we can expect that the galaxy power spectrum and the matter power spectrum are 
related by a linear factor 

P g {k) = b 2 P m (k), (1.31) 
where b is the linear galaxy bias. This model has been verified using N-body simu- 



lations (e.g. Coles 1993 Scherrer & Weinberg 1998). A more physically motivated 
picture of the relation between matter and galaxies is the halo model, which will be 
discussed in section 11.51 



While the assumption of a linear bias will break down for small scales, it also 
means that the galaxy bias b and as are degenerate and the normalisation of the 
galaxy power spectrum alone cannot constrain both of these two parameters. We 
will discuss this further in chapter [3] and we will show possibilities for breaking this 
degeneracy using extra information in the galaxy clustering. 



1.3 Initial conditions and the evolution of the pri- 
mordial power spectrum 

Temperature fluctuations in the Universe today are of the order of 10 8 Kelvirj^J while 
the CMB shows fluctuations of the order of 10~ 5 Kelvin. The growth from these 
very small fluctuations to today's fluctuations can (to some extent) be calculated by 
perturbation theory. To do this we need to solve the coupled Einstein and Boltzmann 
equations. The Einstein equations tell us how perturbations in the metric influence 
the different energy components, while the Boltzmann equations tell us how these 



quantities evolve with time. The Boltzmann equation is (see e.g. Kolb & Turner 
2 Approximate temperature in the core of the sun 
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1994b 



df(x,p,t) 
dt 



C(f), 



1.32) 



where f(x, p, t) is the phase-space density of the different components of the Universe 
(photons, neutrinos, electrons, collision-less dark matter and protons). The left hand 
side of this equation is the collision-less part, while the right hand side, C(f), has 
to contain all the appropriate scattering physics (in the early Universe, the photons 
are coupled to the electrons, which in turn are coupled to the protons). To solve 
this equation we need to perturb the distribution function /, using the perturbed 
metric coefficients, together with their equation of motion derived from the Einstein 
equations (see e.g. 



Peebles & Yu 1970 Wilson k Silk 1981 Bond & Efstathiou 



1984 Zlosnik 2011). The result is a system of coupled differential equations for the 



distribution functions of the non-fluid components (photons, neutrinos and collision- 
less dark matter) together with the density and pressure of the collisional baryon 
fluid and gravitational waves. Given a set of initial conditions, these equations 
are sufficient to describe the evolution of cosmological perturbations. To do this 
we have to integrate the whole system numerically and evolve it to the present. 



Conveniently this can be done using a publicly available code, e.g. CMBeasy Seljak 



& Zaldarriaga 1996 or CAMB Lewis & Bridle 2002 . Some alternative methods 



using analytic solutions and fitting functions have been suggested Bardeen et al. 



1986 Hu & Sugiyama 1996; Eisenstein & Hu 1998 Montanari & Durrer 2011 



but modern cosmology depends on such high precision that most of the time the 
full approach is needed. In this analysis we use the software package CAMB and 



refer the interested reader to Zlosnik 2011 for a detailed derivation and theoretical 
discussion. 



1.3.1 Inflation 

The inflation theory was primarily introduced as an explanation of the two main 
problems of the standard cosmological model, the flatness problem and the horizon 
problem. The horizon problem describes the observational fact that CMB photons 
observed from two opposite directions of the sky have the same mean temperature 
(with overall fluctuations of the order of 10~ 5 Kelvin), although those regions were 
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never in casual contact. 

The flatness problem can easily be discussed by looking at a modified version 
of the Friedmann equation 

(«S " 1)P« 2 = ~- (1-33) 

Flatness is expressed by (0^ — 1) = 0. The right hand side of this expression 
contains only constants, and therefore the left hand side must remain constant 
throughout the evolution of the Universe. As the Universe expands the scale factor 



a increases and the density p decreases (see eq. 1.12). For most of the lifetime of the 
Universe it is matter or radiation dominated and hence p decreases with a -3 and or 4 
respectively. Therefore the factor pa 2 will decrease. Since the time of the Planck era, 
shortly after the Big Bang, this term has decreased by a factor of around 10 60 , and 
so (fi^t — 1) m ust have increased by a similar amount to retain the constant value of 
their product. Nevertheless all measurements are still consistent with (fi^t — 1) = 0- 
Why did the Universe pick flatness with such high precision? 

To avoid the flatness problem, inflation introduces a period in which the term 
p decreases more slowly than a 2 and hence — 1) moves towards zero. This 
also means that regions which are not causally connected today could have been 
connected before the epoch of inflation. It is interesting to note that the future of 
the Universe, which will be dark energy dominated, will show a similar behaviour. 

For a FLRW Universe a rapid expansion follows the criteria 

a _ 3(w+ i) > a - 2 ^ W <_1. (1.34) 

This condition can be obtained by introducing a scalar field <fi minimally coupled 
to gravity. The equation of motion for such a field is given by the Klein-Gordon 
equation 

+ 3#0+^^=O, (1.35) 
d(p 

where V is the potential of the scalar field. The energy density and the pressure of 
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the scalar field are 



1 . 



V 
P 



1 ■ 



(1.36) 



Using the inflation condition w 
to satisfy 



p/p < —1/3 we find that the potential V(4>) has 



2 < v {<!>). 



1.37) 



The Universe during inflation is composed of this uniform scalar field. The scalar 
field energy decays in a quantum process into the energy components we see today 
and since it decays differently in different parts of the Universe, it seeds the density 
fluctuations which then gravitationally grow over time. 

Observations of the Cosmic Microwave Background are now good enough to 



begin to constrain the type of inflational potential (see e.g. Leach & Liddle 2003 



Kinney et al. 2008 Finelli 2010). Here we are going to use the horizon flow pa- 



rameters ei = —d\nH(N)/dN and 62 = dln|ei|/diV jSchwarz, Terrero-Escalante 



& Garcia, 2001 , where N is the number of e-foldings before the end of inflation. 
Assuming that the inflationary phase is driven by a potential of the form V(4>) oc cff 
(chaotic inflation) we can find the following relation between the horizon flow pa- 



rameters, the power-law index, p, and the number of e-folds, N Leach & Liddle 



2003 



62 



N 



-ei, 
p 

P (1_ 
4 U 



- 1 



1.38) 
1.39) 



Figure 1.2 shows one of the latest constraints on these parameters using data from 



WMAP7 Komatsu et al. 2011 and large scale structure data from the BOSS- 



CMASS sample Sanchez et al. , 2012 . The dot-dashed lines correspond to infla- 



tionary models with p = 1, 2 and 4. The blue contours represent the constraints 
from WMAP7 alone. Large scale structure data alone cannot constrain e\ or €2 but 
such data can help to break degeneracies in the CMB. Therefore adding data from 
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Figure 1.2: Constraints on the horizon flow parameters e\ and €2 (see text) from 
WMAP7 combined with the BOSS-CMASS sample (red contours) and WMAP7 
alone (blue contours). This analysis assumes a power law inflational potential 
{V{(j)) oc P ). (credit: |Sanchez et aT] [2012]) 



the BOSS-CMASS sample improves the constraints by about 30%. Using this data 



together with an assumed upper limit for the number of e-folds of N < 60 Dodelson 



& Hui 2003; Liddle & Leach, 2003 , Sanchez et al. 2012 report a limit of p < 1.2 



at the 95% confidence level, imposing a constraint on the inflationary potential used 



in this analysis (see also Komatsu et al. 2009, 2011). Future CMB experiments 
will improve these constraints and hence we can expect to learn much more about 
inflation in the coming decade. 



1.3.2 The Universe dominated by radiation 

Directly after inflation, the Universe entered a phase dominated by the radiation 
density which lasted untill z ~ 3400. Even though the matter energy density sur- 
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passed the radiation density at this point, the Universe remained optically thick to 
radiation until a cosmological redshift of z ~ 1100, when the Universe was about 
380 000 years old. The events during this short period are essential for this thesis 



and will be discussed in detail (see also Hu & Sugiyama 1996 Eisenstein, Seo & 



White 2007). 



The small over-densities seeded by inflation represent not only an excess of 
matter but also an excess of photons and thus such regions are subject to photon 

In the radiation dominated era, the photons are 



Peebles & Yu, 1970 



pressure 

tightly coupled to baryons via Thomson scattering. As a result, these over-pressured 
density peaks initiate sound waves, travelling with the sound speed 



dp c 2 
dp ~ 3(1 + R)' 



1.40) 



with R = 3pb/4p 7 , where p 7 is the photon density and p b is the baryon density. 
The cold dark matter does not interact with photons and hence its perturbations 



grow via gravitational collapse (see Figure 1.3). Therefore CDM perturbations begin 



to dominate the total density perturbations connected to the gravitational potential 
0. The baryon-photon fluid oscillates within these potential wells caused by the 
CDM. 

The evolution of the baryon-photon sound wave can be expressed as a plane 
wave perturbation of wavenumber k. Ignoring the expansion of the Universe and 
assuming that the gravitational potential <fi is caused by CDM only, the evolution 



of the perturbation S b in mode k in the tight-coupling limit follows Peebles & Yu 



1970 Doroshkevich, Zel'dovich & Syunyaev, 1978 Hu & White, 1996 



[l + R)5 b + ^5 b 



[l + R)(-k 2 c 2 - 



;i.4i) 



This equation represents a harmonic oscillator with mass (1 + R) and frequency 
lo = kcj a/3(1 + R) = kc s Hu & Sugiyama 



1996 



. The solution is of the form 



5 b + (1 + R)(j) oc cos(kc s r) 



1.42) 
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with r being the conformal time r = f dt/a. We define 

P P c (t) 

r s (t) = / c s dr) = a Q / -^dt, (1.43) 

Jo Jo °W 

which we call the sound horizon at time t, since it represents the co-moving distance 
sound has traveled by time t. This leads to 

5 b + (1 + R)<j> oc cos(kr s (t)). (1.44) 

For photons we have 

9 + <f) + R(f)(xcos(kc s T), (1.45) 

where 6 = 5T/T and 6 + is the temperature fluctuation of photons after they 
climb up the gravitational potential. 

At redshifts around z ~ 1100 the tight coupling between the baryons and pho- 
tons breaks down as the Universe starts to recombine and the cross section for 
Thomson scattering drops dramatically. At this point the sound-wave stops and the 
physical structure becomes imprinted into the distribution of matter as well as into 
the distribution at photons. From now on the photons travel through the cosmos 
with little interaction and form the Cosmic Microwave Background. 

The oscillation pattern present at decoupling is imprinted in the temperature 
and density distribution of the photons and baryons. However, since the number 
density of photons remains higher after decoupling, than the number density of 
baryons, the small percentage of photons that happen to interact with baryons 
remains sufficient to drive the pressure wave in the baryons until z ~ 1030. The 
epoch of baryon decoupling is called the drag epoch, and the sound horizon r s (t) for 
baryons is ~ 4% larger than for photons. 

The sound horizon for photons is given by r s (t dec ) and marks the wavenumber 
of the first compression at k — n/r s (tdcc)- Subsequent compressions are given by 

kr s (t dec ) = mix (1.46) 

with m = 1,2,3,.... Therefore we see strong structure in the CMB anisotropy at 
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the multipoles 

£ = kd A {t dec) = 

corresponding to these scales. Here 

£ A = n 



_d A (t 



dec) 



r s (i 



ml. 



dec I 



(1.47) 



d A {t 



dec ) 



r s (t 



dec I 



7T 



(1-48) 



is the acoustic scale in multipole space and 



0, = 



r s {t 



dec J 



d A {t 



1.49) 



dec J 



is the sound horizon angle, i.e. the angle at which we see the sound horizon on 
the last scattering surface. Because of these acoustic oscillations, the CMB angular 
power spectrum Cg has a structure of acoustic peaks at sub-horizon scales. We will 
discuss these features further in section 1X761 when we discuss the CMB observations. 

While photon fluctuations before the last scattering surface are dominated by 
density fluctuations, the growing mode of the baryon fluctuations after the drag 
epoch is dominated by velocity fluctuations, i.e. velocity overshoot, not by the 



density term Sunyaev & Zeldovich, 1970 Press & Vishniac, 1980 . As a result, the 



acoustic oscillations in the baryon component are displaced from those in the CMB 
by a phase shift of tt/2 on large scales. 

After the drag epoch, pressure-less baryons quickly fall into the potential wells 
of cold dark matter fluctuations due to gravity. The oscillatory features are diluted 
relative to the CMB because of the small abundance of baryons relative to dark 
matter. 

Whilst the baryon acoustic oscillations are the dominant feature in the CMB 
power spectrum, they are difficult to detect in the matter power spectrum. While 
in Fourier space we have an oscillatory feature, in real-space we expect a single peak 



at the sound horizon scale Eisenstein, Seo & White, 2007 . 
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o 




Figure 1.3: This plot shows the evolution of a k = 1 Mpc -1 mode for different energy 
components. This example represents a cold + hot dark matter model (see text for 
details), (credit: |Ma &: Bertschinger, 1995 with additional labels) 



1.3.3 Growth of perturbations 

After inflation has seeded perturbations in the distribution of matter, the dark 
matter perturbations start to grow through gravitational collapse, while the baryonic 
matter is dominated by its interaction with photons as discussed above. On scales 
larger than the horizon, matter perturbations are driven by the evolution of the 
perturbations in the dominant component of the Universe, which is oc a 2 in the 
radiation-dominated era and oc a in the matter-dominated era. As the Hubble radius 
grows in the expanding Universe, it encompasses larger and larger perturbations. 
Within the horizon, perturbations are subject to different effects which determine 
their growth rate. 



Figure 1.3 shows how the clustering of the different energy components evolves 



with time for a mode at k = 1 Mpc 1 Ma & Bertschinger 1995 . The solid line 



shows the cold dark matter which grows as log(a) in the radiation-dominated area 
(—5 < log 10 (a) < —4.5) and as a in the matter-dominated area log 10 (a) > —4.5. 
Very small wavenumbers k (large scales) enter the horizon in the matter-dominated 
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area and never encounter the suppression of growth during the radiation-dominated 
time. The transition of scales which enter the horizon during radiation domination 
to the scales which enter in the matter-dominated era is given by k eq ~ 0.0 lh Mpc -1 



and marks a peak in the power spectrum (see Figure 1.4) 



The radiation components show strong oscillations caused by the photon-baryon 



fluid. After the decoupling, the baryons (dashed-dotted line in Figure 1.3) fall into 
the potential wells of cold dark matter and soon have a similar clustering amplitude. 
The photons however do not participate in the clustering after decoupling. 

Massless neutrinos behave similarly to photons and do not participate in cluster- 
ing, while massive neutrinos (dashed line) can cluster, depending on their momen- 



tum. In Figure 1.3, the neutrinos are at first, relativistic and behave like massless 
neutrinos, but at some point they become non-relativistic and start to cluster. The 
smaller clustering amplitude of massive neutrinos, which is in conflict with observa- 
tions, is the strongest argument against hot dark matter. 

Furthermore the growth of cold dark matter perturbations during the radiation- 
dominated era is essential for the perturbation amplitude we see today. Without 
dark matter, the epoch of galaxy formation would occur substantially later in the 
Universe than is observed. 

A convenient way to express how the amplitude of a certain mode changes over 
time is the transfer function 

d(k,a e ) 

Here we use a e to denote the scale factor at horizon entry, given by a e ex k~ l . Using 
the transfer function we can express the power spectrum at any time a as a product 
of the initial power spectrum and the transfer function: 

P(k, a) = T 2 (k, a)P(k, a e ). (1.51) 

The transfer function has an asymptotic behaviour of T(k) ~ 1 for small k and 
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k [hMpc' 1 ] k [hMpc 1 ] 

Figure 1.4: Matter power spectrum and transfer function at redshift zero. The power 
spectrum follows the initial power spectrum k Us untill k eq where its power is sup- 
pressed because such modes enter the horizon during the radiation-dominated era. 
The transfer function shows the suppression and is related to the power spectrum 
by P(k) oc k n °T{k). 



T{k) ~ k 2 for large k, with a turning point at 

1 



k 



rh(a eq ), (1-52) 



where r^(a eg ) is the Hubble radius at the time of radiation-matter equality. The 
behaviour T(k) ~ k~ 2 at large k is because the large k modes correspond to small- 
scale fluctuations which entered the horizon at the radiation domination era. At 
that time, the horizon size grows proportionally to the scale factor r/ l (a) ~ a. Since 
a fluctuation enters the horizon when kr^ ~ 1, one has a e oc k~ x . Therefore the 
small-scale fluctuations are suppressed by (a e /a eq ) 2 oc k~ 2 . 



Figure 1.4 shows a matter power spectrum (left) and transfer function (right) 
at redshift zero. On large scales (small k) we have the primordial power spectrum 
which behaves as k ns since the transfer function is 1. The power spectrum peaks 
at matter-radiation equality. The position of the peak (at k eq ) depends on the 
horizon size at matter radiation equality and follows k eq oc VL m h 2 . This means that 
all modes with k > k eq enter the horizon when the Universe was still radiation- 
dominated, where the perturbations can only grow logarithmically, while the modes 
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with k < k eCL enter the horizon when the Universe is matter dominated, where the 
perturbations are able to grow linearly. Hence the small-scale modes are suppressed 
relative to the larger-scale modes. Decreasing Q m will cause the matter radiation 
equality to happen later, the horizon size at matter radiation equality to be larger 
and the corresponding k eq to be smaller (the power spectrum peak moves to the 
left). 



1.4 Redshift space distortions 



Although it is possible to measure true distances for galaxies using standard candles, 
this is very time consuming. Redshifts are much easier to acquire than distances and 
the distortion in the mapping of galaxies in redshift space compared to real space 



allows a statistical measurement of peculiar velocities Kaiser, 1987 Hamilton, 1992 



Fisher et al. , 1994 



The redshift-space position, s, of a galaxy, is related to its real-space position, 



r , by 



r + v z (r)z, 



1.53) 



where v z (f) = u(f)/(aH(z)) is the line-of-sight peculiar velocity with the z-axis 
being the line of sight. For a uniform z-independent mean galaxy density, the exact 
Jacobian for the transformation from real-space to redshift-space is 



dh 
d 3 r 



1 + 



v z (r) 



1 + 



dv z (f) 
dz 



;i.54) 



The galaxy over-density field in redshift-space can be obtained by imposing mass 
conservation, (1 + S s )d 3 s = (1 + 5 r )d 3 r. Furthermore we assume that v z is much 



smaller than the distance to the pair (plane parallel approximation; but see Papai 
& Szapudi, 2008 and the analysis in chapter [2] and [3]) and hence the term v z /z can 
be neglected. We get 

dv z (r) 



5, = S r 



dz 



1.55) 



If we assume an irrotational velocity field, we can write dv z /dz = — (d/dz) V 9, 
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where 9 = — V • v, and V is the Laplacian operator. In Fourier space, (d/dz) 2, V 2 = 
(k z /k) 2 = [i 2 ,, where /i& is the cosine of the line of sight angle, so we have 

S s (k) = 5 r {k) +fi 2 6(k). (1.56) 

Assuming the velocity field comes from linear perturbation theory, we have 

9(k) = -V • v = 6 = 5-^- = Sf, (1.57) 

am a 

where / = din 6/ din a is called the growth rate. We now have a linear relation 
between the redshift-space density field and the real space density field 

S e {k) = {l + fi 2 f)S r {k). (1.58) 

For a population of galaxies with a linear bias b, the linear redshift-space power 
spectrum, P g (k) is then proportional to the linear, real-space matter power spectrum 
Pm\k) 

P*(k, fi) = (b + ffi 2 ) 2 P m (k) = b 2 (l + Pfi 2 ) 2 P^(k), (1.59) 

where we have defined (3 = f/b. Although the real space power spectrum is isotropic, 
this is not true of the power spectrum in redshift-space. Plane waves along the line of 
sight are compressed in redshift-space and so their amplitudes are increased relative 
to those perpendicular to the line of sight. To see this effect in real-space, we can 
define a 2-dimensional correlation function £(r p , tt) with the two separations 

7T = t^l, (1.60) 



r p = yj\h\ 2 -ir 2 , (1.61) 

where we have defined the two vectors s = (si + s 2 )/2 and h = s i — s 2 through the 
positions of two galaxies s*x and s 2 . Now we can write |Hamilton 1992 



£(r p , tt) = Us)V (fi) + U^PM + U^Wfi), (1-62) 
where the Vi are Legendre polynomials, fi = cos(6) and 9 is the angle between r 
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and the line of sight direction ir. The correlation function moments & are given by 

2(3 f (3 2 ~ 

'4(3 4(3 2 ' 



Us) 



3 + 7 



8/3* 

35 



£(r) + -£(r) + -£(r) 



(1.63) 
(1.64) 
(1.65) 



and 



l(r) 



£(r')r' 2 dr' 



£(r')r' dr'. 



(1.66) 
1.67) 



We will discuss more sophisticated models for the 2-dimensional correlation function 
in chapter [3] including wide-angle effects. 

Figure |1.5| shows one of the first detailed measurements of the two-dimensional 



correlation function, £(7T,r p ), using the 2-degree Field Galaxy Redshift Survey Pea- 



cock et al., 2001 . The overall non-circular structure of the signal is caused by the 



linear redshift space distortion effect, which we just discussed. The enhanced signal 
at small r p is the so called "finger of God" effect, which is non-linear in nature and 
hence difficult to model. 



1.4.1 Testing models of gravity 



While galaxy surveys measure (3, the parameter interesting for cosmology is the 
growth rate / = h(3. In linear theory, the growth rate can be predicted for a given 
set of cosmological parameters and allows tests of General Relativity on cosmic 



scales Guzzo et al. , 2008 



Peebles 1980 showed, that / = f{VL m ) and is well approximated by / 



neglecting dark energy. Including dark energy as a cosmological constant, a small 
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Figure 1.5: The 2D correlation function of the 2-degree Field Galaxy Redshift Sur- 
vey. The non-circular shape allows a measurement of the infall parameter /3, which 
can be predicted by linear theory. While the plot uses a to denote the separation 



perpendicular to the line-of-sight, we used r p in this introduction, (credit: Peacock 
et al.||200ip 



correction is needed Lahav et al. 1991 



/(0 TO , ft,/ 



70 



(1.68) 



This clearly shows that the growth rate / mainly depends on the matter density 
Q m (z) and is only weakly dependent on the cosmological constant. This makes 
sense, since the uniform distribution of dark energy does not participate itself in the 
growth of perturbations and influences structure growth mainly through its impact 
on the expansion of the background Universe. 
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Figure 1.6: The growth rate / as a function of redshift z. The plot includes measure- 
ments from 2dFGRS, VVDS and 2SLAQ together with predictions ACDM (solid red 



line) and several alternative models, (credit: Guzzo et al. 2008). 



Linder 2005 and Linder & Cahn 2007 generalised this relation to / ~ ^m( z 



with the growth index 7, which for a ACDM Universe including General Relativ- 
ity is predicted to be 7gr ~ 0.55. This allows us to directly test the standard 
model through the measurement of redshift space distortions in the galaxy distribu- 
tion. Geometrical probes such as supernovae and baryon acoustic oscillations cannot 
alone distinguish between dark energy as a new negative pressure component of the 
Universe or a modification of laws of gravity. Redshift space distortions therefore 
offer a powerful way to complement geometrical probes by directly testing the laws 
of gravity on the largest scales possible. 



Figure 1.6 shows the first test of General Relativity using the growth rate / 



from Guzzo et al. 2008 . The plot includes three data points from different galaxy 



surveys which all seem to be in agreement with the prediction of ACDM. Since this 
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r p (h _1 Mpc) 



Figure 1.7: The projected correlation function of a volume limited sample from 



SDSS together with the best fitting HOD model, (credit: Zehavi et al. 2004) 



first analysis, new data have been collected, but up to now all such measurements 



seem to be in agreement with the standard model (see e.g. Rapetti et al. 2009 



Samushia et~al] |2011] ; |Blake et al.| |2011a] ; |Reid et al.| |2012] ; |Hudson fc Turnbull 



2012 ; Rapetti et al. 2012 ; Basilakos 2012 and this thesis). 



1.5 The halo model 



Jing, Mo & Borner 


1998 


Ma & Fry 


2000 


Pea- 



cock & Smith 2000; Seljak 2000 Berlind & Weinberg 2001 Cooray & Sheth 2002) 



to galaxy clustering assumes that large-scale structure can be described by a distri- 
bution of dark matter halos while galaxies themselves form within these halos. Such 
models have been tested against n-body simulations and seem to be able to capture 
the galaxy-mass relation. 
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The contribution to the power spectrum (or correlation function) can be written 
as a sum of correlations of galaxies that occupy two different halos and the correlation 
of galaxies within a single halo 



p B {k) = pi h (k)+if{k). 



1.69) 



The exact definition of the one-halo and two-halo terms will be discussed in chap- 
ter n 

The halo model is based on external inputs, such as the dark matter halo profile 
and the halo mass function, which come from n-body simulations. The halo model 
then needs to provide a strategy for how to populate dark matter halos. Usually 
the assumption is that above a certain mass threshold the halo will have a central 
galaxy, and the probability to have a satellite galaxy grows with mass. 

The first measurements of galaxy clustering in the correlation function showed 



that it follows a power law Totsuji & Kihara, 1969 Peebles, 1973; Hauser & Peebles 



1973, 1974 Peebles, 1974 . The halo occupation distribution (HOD) model however 



predicts clear deviations from a power law, especially at the transition scale, between 
the one-halo and two halo terms, above which galaxy pairs sit in different dark matter 



halos. Figure |1.7| shows the correlation function of a volume limited sample of the 
SDSS survey together with the best fitting HOD model. The HOD model is able to 
explain the structure seen in the correlation function, most obviously the transition 
from the one-halo term to the two halo term at around 2/i _1 Mpc. We will discuss 
such an analysis using data from the 6dFGS in chapter [4} 



1.6 Observational cosmology 

In order to test the geometry of the Universe we have to measure the redshift- 
distance relation, since this relation is directly connected to the different components 
of the Universe. For this purpose one needs objects whose intrinsic luminosity or 
size is known. The first class of objects are called "standard candles", while the 
second are known as "standard rulers". The currently best known objects to qualify 
as standard candles are Type la supernovae, while a very good standard ruler is 
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Figure 1.8: The CMB temperature distribution as measured by WMAP7. To make 
the very small temperature fluctuations visible, the mean CMB temperature of 
Tomb — 2.75 Kelvin as well as the CMB dipole and the Galactic emission have 
been subtracted, (credit: NASA) 



provided by the sound horizon, i.e. the distance a sound wave travelled between the 
Big Bang and the baryon drag epoch. The standard ruler technique is an important 
part of this thesis and will be discussed in more detail in the following sections. 
The standard candle technique is only used to compare cosmological constraints 
from different models and we refer the reader to the literature for a more detailed 



discussion of this method (e.g. Perlmutter 2003). 



1.6.1 The Cosmic Microwave Background 

We have already discussed the physics before recombination, after which the photons 
travel with minimal interaction until today. When the photons decoupled from 
the baryons, they travelled out of the small matter perturbations which existed at 
that time. Although these perturbations were small, they caused a gravitational 
redshift or blueshift of the photons [Sachs fc Wolfe 1967 , which is still visible 



in the distribution of photons today (see Figure 1.8). By observing the photon 
distribution we can therefore learn about the very early matter distribution, which 
includes baryon acoustic oscillations. 

The CMB temperature anisotropy is a function over a sphere. We can separate 
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Figure 1.9: The CMB angular power spectrum from WMAP7 |Komatsu et al.[|2011 
and SPT 



Keisler et al. , 2011 



These two measurements together allow the detection 
of up to 8 peaks as well as CMB lensing. At large £ the data deviate from the simple 
model (dashed line) because of emission from background sources like clusters and 



starburst galaxies, (credit: Keisler et al. 2011) 



out the contributions of different angular scales by doing a multipole expansion, 

5T(e,</>)=^2Y im (e,<f>)a bn , (1.70) 

£m 

where the sum runs over the multipole £ = 1,2, ...oo and m = —£, ...,£, giving 2£ + 1 
values of m for each £. The functions Ye m (0,<j)) are the spherical harmonics, which 
form an orthonormal set of functions over the sphere, so that we can calculate the 
multipole coefficients ai m from 

a lm = jY; m (6 1 <f ) )5T(e i( f ) )dn. (1.71) 

Since each a im represents a deviation from the average temperature, their expecta- 
tion value is zero, 

(a tm ) = (1.72) 
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and the quantity we are interested in is the variance (\ae m \ 2 ) which gives the typical 
size of the a£ m . The ai m are independent of m (in an isotropic Universe), since m 
represents just the orientation, which then allows us to define 

C t = (\a em \ 2 ) = J^iM 2 )- ( L73 ) 

m 

For Gaussian perturbations, the values of Ce contain all the statistical information 
about the CMB temperature anisotropy. The function (of integers I > 1) is called 
the angular power spectrum, and this is the quantity which usually is compared to 
observations (analogous to the power spectrum P(k) of density perturbations). 

The angular power spectrum £(£ + l)Ce/2ir as measured by WMAP7 and SPT 



together with the best fitting ACDM model is shown in Figure |1.9| as a function of 
multipole L 



1.6.2 Baryon Acoustic Oscillations as a standard ruler 

The important consequence of the linear gravitational instability theory is that 
all the features present in the initial matter fluctuation spectrum should survive 
throughout cosmic evolution. The baryon acoustic oscillations which are responsible 
for the characteristic features in the CMB power spectrum should be visible in the 
low redshift matter distribution. 

The co-moving size of an object or a feature at redshift z in the line-of-sight (tt) 
and transverse (r p ) direction is related to the observed redshift size Az and angular 
size A9 by the Hubble parameter H(z) and the angular diameter Da(z): 

cAz 

(1.74) 

r p = (l + z)D A {z)A6. 

Thus the measurement of the observed dimensions, Az along the line-of-sight direc- 
tion and AO in the transverse direction, give 7rH(z) and t p /Da{z). When the true 
physical scale of the object or feature, ir and r p , are known, measurements of the 
observables, Az and A9, give estimates of H(z) and Da(z). Such an object is called 
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Figure 1.10: The large scale correlation function of the BOSS-CMASS sample. 



(credit: Anderson et al. 2012) 



a "standard ruler r 



In the case that the signal-to-noise is not sufficient to separate the two signals 
parallel and perpendicular to the line-of- sight, we can constrain a combination of 



both given by Eisenstein et al., 2005 Padmanabhan & White, 2008 



Dv(z) 



;i + zfD\{ z ) 



cz 
Hjz)\ 



1/3 



;i.75) 



Currently we don't know of any process which can alter the sound horizon scale since 
decoupling by more than 1% and hence such a measurement has a low systematic 
uncertainty 



Eisenstein et al., 2007; Crocce & Scoccimarro, 2008 Sanchez et al. 



2008 Matsubara, 2008b; Seoetal., 2008 Smith et al. , 2008 Padmanabhan & White 



2009 



The first measurements of the baryon acoustic peak were performed using the 
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SDSS-LRG sample Eisenstein et al. 2005 and the 2dFGRS sample Cole et al. 



2005 . More recently similar measurements haven been made in other galaxy surveys 



like 6dFGS (this thesis), WiggleZ Blake et al. 2011c and BOSS Anderson et al. 



2012 . The correlation function of the BOSS-CMASS sample is shown in Figure 1.10 
The BAO peak is clearly visible at around r m 100/i _1 Mpc. 

Currently this technique is limited by the statistical precision of the BAO signal 



detection in the galaxy distribution, but ongoing surveys like BOSS Schlegel, White 



& Eisenstein, 2009 aim for 1% constraints on H(z) and Da(z). 
We will discuss the 6dFGS BAO analysis in chapter [2} 



1.6.3 Combining different cosmological probes 

Different cosmological probes are sensitive to different cosmological parameters. For 
example, supernova (SN) data are particularly good for constraints on the dark en- 
ergy equation of state parameter especially if this parameter is varied as a function 



of redshift (see Figure 1.11). The reason for this is that SN la data sample many 
different redshifts and hence can detect variations in w(z) more precisely. However 
SN data are now limited by systematic uncertainties, which are difficult to overcome, 
while BAO data are still statistically limited. With BAO data rapidly improving, 
we expect the BAO data to eventually dominate over SN data, while the Cosmic 
Microwave Background will remain unchallenged as the most powerful cosmological 
probe. Compared to SN, BAO data are more sensitive to the curvature parameter, 
because they have a very large lever arm from low redshift to the era of decou- 
pling. Different cosmological probes can easily be combined using the product of 
the likelihoods or the sum of the \ 2 



2 _ 2 1 2 1 2 

Xtot — XCMB i XbAO i* XSN- 



1.76) 



with 



X 



^^data ^modcl^ ^ 



(J- 



1.77) 



The likelihood contours of such a combination of different probes are shown in 
Figure 1.11 for w(z) and in Figure [Li] for and O m . 
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Figure 1.11: Constraints in w(z) from BAO, CMB and SN data. Because SN data 
sample the redshift distribution with more data-points, the SN data contribute 
strongly to this constraint. BAO data are expected to become dominant in the 
future, because SN samples have reached strong systematic limitations, which are 



difficult to overcome, while BAO data are still statistically limited, (credit: Blake 
et al.||2011cb 



It is reassuring that the standard model of cosmology has now been confirmed 
by several different independent probes. For example BAO data alone now require 
dark energy Blake et al. 2011b and since CMB lensing has been detected, the CMB 



data alone also points towards dark energy van Engelen et al. , 2012 . High precision 
constraints can be obtained by combining any two of the three main cosmological 
probes, SN la, BAO or CMB. Hence, although the CMB remains the most power- 
ful cosmological probe, even excluding the CMB data, the remaining observational 



probes still give strong evidence for the standard model of cosmology Percival et 
aH|20T0 ■ 
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1.7 Alternative models of dark energy 



The interpretation of the accelerated expansion of the Universe as vacuum energy 
has several theoretical weaknesses. For example, it is not clear why we happen 
to live at a time when dark matter and dark energy have approximately the same 
energy density, while the Universe was dominated by only one energy component 
during most of its lifetime and will be in its future (coincidence problem). 

It is also not clear why the vacuum energy is so small. Quantum field theory 
predicts a level of vacuum energy 120 orders of magnitude larger than observed. It 
was always clear that the Quantum field theoretical prediction would not lead to a 
Universe which would allow life to develop and hence some unknown process must 
cancel out vacuum energy. While this situation should have been already unsettling, 



the discovery of the accelerated expansion of the Universe in 1998 Perlmutter et al. 



1998 Riess et al., 1998 and its consistency with a cosmological constant implied 



that vacuum energy did not cancel out completely, but an imbalance of one part in 
10 120 remains. 

This situation has led to the suggestion of many alternative models. In this 
section we will outline the three main categories: modified gravity models, scalar 
field models and inhomogeneous models, grouped by which part of the standard 
model they propose to modify. We will focus on the basic idea of these models and 



refer to recent reviews of this topic for a more comprehensive discussion Copeland, 



Sami fc Tsujikawa" 2006 Durrer Maartens , 2008 Frieman, Turner fc Huterer[ 



2008 Tsujikawa 2010 Li, Wang & Wang, 2011 Clifton, 2011 . 



1.7.1 Scalar Fields: Quintessence 

Scalar fields naturally arise in particle physics and can act as candidates for dark 
energy. Quintessence represents a scalar field that is coupled with gravity. Given 
a particular potential, quintessence can describe the late-time acceleration of the 
expansion of the Universe. The action for quintessence is given by 

(1.78) 



S 



;(V0) 2 -V(0) 
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where g = det(<7 M „), with the metric tensor g^ u , (V0) 2 = g^ u d^(j)d u (j) and V(<f)) is the 
potential of the field. In a flat FLRW Universe this action varies as 



3H(f) + 



dV 







;i.79) 



with respect to 0, where homogeneity of the scalar field is assumed. The stress- 
energy tensor has a form identical to that of an ideal fluid with pressure and density 



identical to the one we found for the case of inflation (see eq. 1.36) 



P=l<P 2 -V(<f>). 
The equation of state parameter for the field <fi is given by 



;i.8o) 



p 
p 



1 2 <P 2 + VW 
2 - v{4>) ' 



;i.8i) 



which suggests that — 1 < w < 1. If the time evolution is slow w ~ — 1, and the 
field behaves like a slowly varying vacuum energy. 

Because of the lack of observational constraints, there is no reason to favour 
one form of the potential with respect to others. For example, the original scenario 
proposed by Ratra & Peebles 1988 was a potential of the form V(<f)) = M 4+a /(f) a , 



whereas Caldwell & Linder 2005 proposed a quintessence theory consisting of two 



classes of thawing and freezing models, depending on whether the field accelerates 
or decelerates with time. Another approach involves a modification of the kinetic 



term, as in k-essence models Armendariz-Picon, Mukhanov & Steinhardt, 2001 



Dynamical models can generally deliver some answers to the question of the 
nature of dark energy, but do not yet have the capability to provide a complete so- 
lution. Some classes of models have the capability to solve the coincidence problem, 
but the smallness of the cosmological constant (or minimum of the potential in this 
case) means that the fine tuning problem remains. 
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1.7.2 Modified Gravity: f(R) 



Here we discuss models which modify the geometrical part of the Einstein equation. 
The f(R) modified gravity theories replace the standard Einstein-Hilbert action by 
an arbitrary function of the Ricci scalar R. In Riemannian geometry, R describes 
the simplest curvature invariant of a Riemannian manifold. It assigns a single real 
number to each point on the manifold which is determined by the intrinsic geometry 



of the manifold near that point Lobo , 2008 



The general action for a modified gravity field is given by 



S 



1 



167rG 



d 4 x^[f(R) + C r , 



1.82) 



where C m is the matter Lagrangian density. Variation of this action with respect to 
the metric g^ u yields the field equation 



FR 



V^V U F + g^UF = 8-kGT 



m 



1.83) 



where F = df(R)/dR, □ = <9 M <9 M = d 2 /dt 2 — V 2 is the d'Alembert operator and T t 
is the matter stress-energy tensor. This equation can be written as 

1 



R 9fiu = 8nGTf u , 



1.84) 



where the new term = T^ u + T™/F describes the effective stress-energy ten- 
sor. All of the effects of the modification of gravity are now contained within the 
curvature part of the tensor T° . 

For a FLRW metric, the generalised Friedmann equations have an identical form 
to eqs. 1.4 and 1.5 with p e ff = p m +p c and pg/y = Pm + Pc- The pressure and density 



of the curvature field are given by: 

1 



Pc 
Pc 



8%GF 
1 

8nGF 



1 



2-Rf + RF' + R F" ~ 7,[f ~ RF] 



l -[f-Rf]-£ R F> 



1.85) 
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Given the right choice of f{R), these models can reproduce the late-time acceleration 
of the Universe. 

While f(R) models are very general, they can provide viable models to compare 



to observations, and many f(R) theories have already been ruled out Bertotti, less 



& Tortora, 2003; Amendola & Tsujikawa, 2008 Hu & Sawicki, 2007 



1.7.3 Inhomogeneous Models 



Using the Friedmann equations for a Universe which is clearly not homogeneous 
on small scales could introduce a bias into cosmological parameter constraints. We 
have to consider that non-linear processes could have a back-reaction effect on the 
background cosmology. Furthermore it is not obvious how we have to perform a 
covariant and gauge invariant averaging over the inhomogeneous Universe to arrive 
at the correct FLRW background. 

It has been claimed in recent years that for certain inhomogeneous models, 



such effects are large enough to mimic an accelerating Universe (see e.g. Buchert 



2000, 2001; Maartens 2011 Ellis 2011). This would be a satisfying resolution of the 



coincidence problem without the need for any dark energy field. However in many 
cases it also means that the Cosmological Principle is not correct. 

In perturbation theory, we introduce small perturbations to a homogeneous 
metric. In an inhomogeneous model, it is assumed that the perturbations are too big 
to just treat them as linear perturbations and hence the exact solutions of Einstein's 
field equations need to be used. The best known example of such an exact solution 



is the Lemaitre-Tolman-Bondi metric Lemaitre 1933; Tolman, 1934 Bondi 1947 



which describes a spherical cloud of dust (finite or infinite) that is expanding or 
collapsing under gravity. A spherical under-density of sufficient size could explain 
the discrepancy between nearby and distant supernovae luminosities without dark 



energy Pascual- Sanchez 1999; Vanderveld, Flanagan & Wasserman, 2006; Enqvist 



2008 



These proposals have been disputed, and so far we can not say that there is a 
convincing demonstration that acceleration could emerge naturally from nonlinear 
effects during structure formation. Even if such effects are not the cause of the 
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accelerating expansion of the Universe, we should note that back-reaction effects 
could significantly affect our estimations of cosmological parameters, even if they do 



not lead to acceleration Li & Schwarz, 2007 



1.8 The 6dF Galaxy Survey 

This thesis presents new constraints on the cosmological model by analysing data 
from the 6-degree Field Galaxy Survey, a near- infrared selected (JHK) redshift 
survey of about 125 000 galaxies across almost the complete southern sky, with 
secondary samples selected in bj and r F . The J,H and K bands exclude a ±10° band 
along the Galactic plane to minimise extinction and foreground source confusion, 
while b j and r> exclude a 20° band. 

The total sky coverage is about 17000deg 2 and the median redshift of the 
sample is z = 0.05. The galaxies were selected from the Two-Micron All-Sky Survey 



- Extended Source Catalog [2MASS XSC; Jarrett et al. , 2000 . The spectra were 
taken with the 6-degree Field (6dF) multi-fibre instrument on the UK Schmidt 
Telescope from 2001 to 2006. The mean completeness of 6dFGS is 92% but varies 



with sky position and magnitude (see Jones et al. 2006 for details). 



For our analysis we mostly used the K-band selected sample, since it is by far 
the largest sample in 6dFGS. For the BAO analysis in chapter [2] we used a faint 
magnitude limit of K — 12.95, while for the redshift space distortion analysis in 
chapter [3] we used a limit of K = 12.75. The motivation for the higher limit in 
the BAO analysis was to maximise the effective volume by selecting more faint, 
high redshift galaxies. For the HOD analysis in chapter [4] we used the J-band 
selected sample, because this band is particularly robust against background noise 
and hence provides the best stellar mass estimates. The details of the samples used 
for the different analysis techniques will be explained at the beginning of each of the 
following chapters. 

A subset of early-type 6dFGS galaxies (approximately 10 000) have measured 
line-widths that will be used to derive Fundamental Plane distances and peculiar 
velocities. This will lead to the largest peculiar velocity dataset ever produced and 



38 



CHAPTER 1. INTRODUCTION 




Dk - if, 




■li 11/ Ml *i) aJ In. U. IL - J3 



Figure 1.12: Density of 6dFGS target sources (per square degree) on the sky; key 
supercluster over-densities are labelled, (b) Full 6dFGS field coverage (filled discs) 
and unobserved target fields (open circles), (c) Redshift completeness for K < 12.65. 



All panels show equal-area Aitoff projections, (credit: | Jones et aL 2009) 
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will allow bulk flow measurements as well as tests of General Relativity. 



1.8.1 Comparison to Other Surveys 



The 6dF Galaxy Survey covers approximately ten times the area of the 2dF Galaxy 
Redshift Survey (2dFGRS; |Colless et aj]|200j) and has more than twice the areal 



coverage of the Sloan Digital Sky Survey seventh data release (SDSS DR7; Abazajian 



et al. 2009). In terms of secure redshifts, 6dFGS has around half the number of 
2dFGRS and one-sixth those of SDSS DR7 (r < 17.77). The co-moving volume 
covered by 6dFGS is about the same as 2dFGRS at their respective median redshifts, 
and around 30% that of SDSS DR7. A sub-sample of about 60 000 luminous red 



galaxies (LRGs) jEisenstein et al. 2001 in SDSS covers a volume of about lh~ 3 Gpc 3 
at redshift z ~ 0.35, with the main focus on cosmological measurements like baryon 
acoustic oscillations. 

In terms of fibre aperture size, the larger apertures of 6dFGS (6.7") give a 
projected diameter of 4.8/i _1 kpc at the median redshift of the survey, covering 40% 
more projected area than SDSS at its median redshift, and more than three times 
the area of 2dFGRS. 



The WiggleZ Dark Energy Survey Drinkwater et al., 2010 measured the red- 



shifts of about 250 000 galaxies over a volume comparable to the SDSS-LRG sample 
but at higher redshift of around z ~ 0.6. This survey focuses on the measurement of 
dark energy through the detection of the baryon acoustic oscillations, similar to the 
still ongoing Baryon Oscillation Spectroscopic Survey (BOSS) [Schlegel, White fc 



Eisenstein 


2009 


stein et al. 


2011 



BOSS has a mean redshift of z ~ 0.5 and will collect 1.5 million 
galaxy redshifts over a volume of about 6h~ 3 Gpc 3 . It will also measure baryon 
acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of 



150 000 quasars at redshifts z > 2.2 Ross et al. 2011 . 
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1.9 Overview and motivation for this thesis 



In the next chapter (chapter[2]) we will discuss the baryon acoustic oscillation analysis 
of 6dFGS. The detection of the acoustic feature allows an absolute distance mea- 
surement to the effective redshift of the survey z c g = 0.1. At such a low redshift, 
the distance is independent of most cosmological parameters, except the Hubble 
constant H . Therefore our analysis allows a robust measurement of this parameter 
similar to the distance ladder technique. The systematic uncertainties of the BAO 
technique to derive the Hubble constant are very different to the distance ladder 
technique and hence allow an independent crosscheck. It is also possible to break 
degeneracies in the WMAP7 dataset by combining it with our measurement. 

The redshift of 6dFGS is much lower than other galaxy surveys which reported 
BAO detections so far and hence delivers a new data point on the BAO Hubble 
diagram. Combining BAO measurements at different redshifts allows us to detect 



dark energy by just using BAO data, without combining it with the CMB Blake et 



al., 2011b 



The error in our BAO measurement is sample variance limited and hence bet- 
ter constraints are possible with larger surveys. At the end of chapter [2] we include 
predictions for two future low redshift surveys, the radio galaxy survey, WALLABY, 
which will be conducted using the Australian SKA pathfinder (ASKAP) telescope 
in Western Australia and the TAIPAN survey, which is planned for the UK Schmidt 
telescope at Siding Spring Observatory. With an effective volume three times that 
of 6dFGS, TAIPAN will be able to constrain the Hubble constant to 3% using the 
BAO technique, significantly improving upon the 6dFGS result. 

In chapter [3] we use the amplitude of redshift-space distortions to measure the 
parameter combination fa 8 , where / = din D/ din a is the linear growth rate and 
erg is the r.m.s. of matter fluctuations in 8/i _1 Mpc spheres. The expected growth 
rate can be derived from General Relativity, allowing a test of this theory on cosmic 



scales. As we explained in section [1.4.1 such tests might reveal more information on 
the nature of dark energy. With 6dFGS, we can probe a new redshift range, which 
has not been investigated by other galaxy surveys so far. The simplest ACDM 
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parameterisation of the growth rate / is given by / = Q^z), where the growth 
index, 7, is the parameter we are interested in. Within this simple parameterisation 
we can see that we are most sensitive to 7 when Q m (z) is small. The minimum 
of Q m (z) is at z = with Q m (z) ~ 1 at redshift z > 1. This means that low 
redshift constraints on the growth rate are the most powerful in measuring 7. This 
is not necessarily true for other parameterisations of the growth rate, but other 
parameterisations usually depend on more parameters and as long as we haven't 
detected any deviation from ACDM, and without a convincing alternative model of 
gravity, it is practical to follow the simplest parameterisation. 

A further advantage of low redshift constraints on the growth rate is their rela- 
tive independence from the Alcock-Paczynski (AP) effect. The AP effect originates 
from the assumption of a fiducial cosmological model to turn redshifts into distances. 
If this initial assumption is wrong, this will introduce a distortion into the 2D cor- 
relation function, which is very similar to the linear redshift-space distortion signal. 
While the AP effect is interesting as an additional tool to test cosmological models, 
it leads to a degeneracy with /erg. In 6dFGS the redshift distance relation is almost 
completely independent of the fiducial cosmological model and hence the AP effect 
is negligible. 

Galaxy redshift surveys alone cannot constrain the growth rate / because there 
is a degeneracy between the galaxy bias, b, and c 8 . This is the reason why most 
galaxy surveys report constraints on fa 8 or instead of the actual parameter of 
interest, /. In our analysis we will combine our measurement of fag with WMAP7, 
but we also investigate alternative methods to directly constrain / without addi- 
tional datasets. 

Again we include forecasts for constraints on /c 8 expected from the WALLABY 
and TAIPAN surveys. Although WALLABY has a smaller effective volume, com- 
pared to TAIPAN, it will deliver better constraint on /erg, since it has a much smaller 
galaxy bias. 

Chapter |4] uses the halo model to investigate the relation between dark matter 
clustering and galaxy clustering. The question is how can we relate the observable 
properties of galaxies to the not-observable properties of dark matter halos? Studies 
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have shown that stellar mass is tightly related to dark matter halo mass (e.g. Wang 



et al. 2007| ). Hence we chose to investigate galaxy clustering as a function of 



this observable in 6dFGS. The near-infrared selection of 6dFGS allows very reliable 



estimates of stellar masses, more robust against background noise Bell & De Jong 



200T| |Drory et al.[ [20041 |Kannappan fc Gawiser] [20071 |Longhetti fc Saracco] [2009 



Grillo et al. , 2008 Gallazzi & Bell , 2009 compared to other galaxy surveys which rely 



on optical bands. We will discuss this aspect in more detail in chapter [4} We divide 
6dFGS into four volume-limited sub-samples with thresholds in stellar mass and 
redshift and calculate the projected correlation functions w p {r p ). We then fit these 
correlation functions with our HOD model. Our analysis allows us to investigate 
the effective dark matter halo mass for a given stellar mass as well as the satellite 
fraction. 

We compare our results to semi-analytic catalogues derived from the Millennium 
Simulation. Semi-analytic catalogues attempt to model baryonic effects, such as gas 
cooling and supernova feedback, with analytical descriptions placed on top of a 
dark matter only N-body simulation. The analytical descriptions are calibrated 
to reproduce different observables, like the galaxy luminosity function or stellar 
mass function. By comparing observational results with semi-analytic models we 
can investigate whether the different properties of these models are a satisfying 
description of all processes influencing the formation and evolution of galaxies. Such 
information can than be used to improve upon these semi-analytic models, with the 
goal of furthering our understanding of galaxy formation. 

We summarise and conclude in chapter [5] 
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Chapter 2 

Baryon Acoustic Oscillations and 
the Local Hubble Constant 



Beutle r et al.\ 
MNRAS 416, 3017B (2011) 



Abstract 

We analyse the large-scale correlation function of the 6dF Galaxy Survey (6dFGS) and 
detect a Baryon Acoustic Oscillation (BAO) signal. The 6dFGS BAO detection al- 
lows us to constrain the distance-redshift relation at z e g = 0.106. We achieve a dis- 
tance measure of Dy(z e s) = 456 ± 27 Mpc and a measurement of the distance ratio, 
r s (zd)/ Dy(z e fi) = 0.336 ±0.015 (4.5% precision), where r s (zd) is the sound horizon at the 
drag epoch Zd- The low effective redshift of 6dFGS makes it a competitive and independent 
alternative to Cepheids and low- z supernovae in constraining the Hubble constant. We 
find a Hubble constant of Hq = 67.0 ± 3.2 km s _1 Mpc -1 (4.8% precision) that depends 
only on the WMAP-7 calibration of the sound horizon and on the galaxy clustering in 
6dFGS. Compared to earlier BAO studies at higher redshift, our analysis is less depen- 
dent on other cosmological parameters. The sensitivity to Hq can be used to break the 
degeneracy between the dark energy equation of state parameter w and Hq in the CMB 
data. We determine that w = —0.97 ± 0.13, using only WMAP-7 and BAO data from 
both 6dFGS and |Percival et al.| (20To) . 

We also discuss predictions for the large scale correlation function of two future wide-angle 
surveys: the WALLABY blind HI survey (with the Australian SKA Pathfinder, ASKAP), 
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and the proposed TAIPAN all-southern-sky optical galaxy survey with the UK Schmidt 
Telescope (UKST). We find that both surveys are very likely to yield detections of the 
BAO peak, making WALLABY the first radio galaxy survey to do so. We also predict 
that TAIPAN has the potential to constrain the Hubble constant with 3% precision. 



2.1 Introduction 

The current standard cosmological model, ACDM, assumes that the initial fluctu- 
ations in the distribution of matter were seeded by quantum fluctuations pushed 
to cosmological scales by inflation. Directly after inflation, the Universe is radia- 
tion dominated and the baryonic matter is ionised and coupled to radiation through 
Thomson scattering. The radiation pressure drives sound-waves originating from 



over-densities in the matter distribution Peebles & Yu , 1970 Sunyaev & Zeldovich 



1970 Bond & Efstathiou, 1987. At the time of recombination (z* ~ 1090) the 
photons decouple from the baryons and shortly after that (at the baryon drag epoch 
Zd Rj 1020) the sound wave stalls. Through this process each over-density of the 
original density perturbation field has evolved to become a centrally peaked per- 



turbation surrounded by a spherical shell Bashinsky & Bertschinger 2001, 2002 



Eisenstein, Seo & White, 2007 . The radius of these shells is called the sound 
horizon r s . Both over-dense regions attract baryons and dark matter and will be 
preferred regions of galaxy formation. This process can equivalently be described in 
Fourier space, where during the photon-baryon coupling phase, the amplitude of the 
baryon perturbations cannot grow and instead undergo harmonic motion leading to 
an oscillation pattern in the power spectrum. 

After the time of recombination, the mean free path of photons increases and 
becomes larger than the Hubble distance. Hence from now on the radiation re- 
mains almost undisturbed, eventually becoming the Cosmic Microwave Background 
(CMB). 

The CMB is a powerful probe of cosmology due to the good theoretical under- 
standing of the physical processes described above. The size of the sound horizon 
depends (to first order) only on the sound speed in the early Universe and the age 
of the Universe at recombination, both set by the physical matter and baryon den- 
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sities, Q m h and VL^h Eisenstein & Hu, 1998 . Hence, measuring the sound horizon 



in the CMB gives extremely accurate constraints on these quantities Komatsu et 



al. 2011 . Measurements of other cosmological parameters often show degeneracies 



in the CMB data alone Efstathiou & Bond, 1999 , especially in models with extra 



parameters beyond flat ACDM. Combining low redshift data with the CMB can 
break these degeneracies. 

Within galaxy redshift surveys we can use the correlation function, £, to quantify 
the clustering on different scales. The sound horizon marks a preferred separation of 
galaxies and hence predicts a peak in the correlation function at the corresponding 
scale. The expected enhancement at s = r s is only A£ 10~ 3 6 2 (1 + 2/3/3 + /3 2 /5) 
in the galaxy correlation function, where b is the galaxy bias compared to the mat- 
ter correlation function and /3 accounts for linear redshift space distortions. Since 
the signal appears at very large scales, it is necessary to probe a large volume 
of the Universe to decrease sample variance, which dominates the error on these 



scales Tegmark, 1997; Goldberg & Strauss, 1998 Eisenstein, Hu & Tegmark, 1998 



Very interesting for cosmology is the idea of using the sound horizon scale as a 



standard ruler Eisenstein, Hu & Tegmark, 1998 Cooray et al. 2001 Seo & Eisen- 



stein 



2003; Blake & Glazebrook, 2003 



A standard ruler is a feature whose absolute 
size is known. By measuring its apparent size, one can determine its distance from 
the observer. The BAO signal can be measured in the matter distribution at low 
redshift, with the CMB calibrating the absolute size, and hence the distance-redshift 



relation can be mapped (see e.g. Bassett & Hlozek 2009 for a summary). 



The Sloan Digital Sky Survey (SDSS; York et al. 2000), and the 2dF Galaxy 



Redshift Survey (2dFGRS; Colless et al. 2001) were the first redshift surveys which 



have directly detected the BAO signal. Recently the WiggleZ Dark Energy Survey 
has reported a BAO measurement at redshift z = 0.6 Blake et al. 2011b . 



Eisenstein et al. 2005 were able to constrain the distance-redshift relation to 
5% accuracy at an effective redshift of z b q = 0.35 using an early data release of 
the SDSS-LRG sample containing w 47000 galaxies. Subsequent studies using the 
final SDSS-LRG sample and combining it with the SDSS-main and the 2dFGRS 
sample were able to improve on this measurement and constrain the distance-redshift 
relation at z c g = 0.2 and z e fr = 0.35 with 3% accuracy Percival et al. 2010 . Other 
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studies of the same data found similar results using the correlation function £(s) 



|Martinez et al. 2009 C 


laztanaga, Cabre & Hui 2009 Labini et al. 


2009 


; Sanchez 


et al. 


2009; 


Kazin et al. , 


2010 , the power spectrum P(k) |Cole et al. 


2005 


Tegmark 


et al. 


2006 


Huetsi 2006 


Reid et al. , 2010 , the projected correlation function w(r p ) 


of photometric redshift samples Padmanabhan et al. , 2007 Blake et al. , I 


2007| and 



a cluster sample based on the SDSS photometric data Huetsi, 2009 . Several years 



earlier a study by Miller, Nichol & Batuski 2001 found first hints of the BAO 
feature in a combination of smaller datasets. 

Low redshift distance measurements can directly measure the Hubble constant 
Hq with a relatively weak dependence on other cosmological parameters such as the 
dark energy equation of state parameter w. The 6dF Galaxy Survey is the biggest 
galaxy survey in the local Universe, covering almost half the sky. If 6dFGS could be 
used to constrain the redshift-distance relation through baryon acoustic oscillations, 
such a measurement could directly determine the Hubble constant, depending only 
on the calibration of the sound horizon through the matter and baryon density. The 
objective of the present paper is to measure the two-point correlation function on 
large scales for the 6dF Galaxy Survey and extract the BAO signal. 

Many cosmological parameter studies add a prior on Hq to help break degen- 
eracies. The 6dFGS derivation of Hq can provide an alternative source of that 
prior. The 6dFGS if - measuremen t can also be used as a consistency check of other 
low redshift distance calibrators such as Cepheid variables and Type la supernovae 



(through the so called distance ladder technique; see e.g. Freedman et al. , 2000 



Riess et al. , 2011 ). Compared to these more classical probes of the Hubble constant, 
the BAO analysis has an advantage of simplicity, depending only on Q m h 2 and Vt^h 2 
from the CMB and the sound horizon measurement in the correlation function, with 
small systematic uncertainties. 

Another motivation for our study is that the SDSS data after data release 3 
(DR3) show more correlation on large scales than expected by ACDM and have 
no sign of a cross-over to negative £ up to 200/i _1 Mpc (the ACDM prediction is 
140/i" 1 Mpc) |Kazin et alj |2010j . It could be that the LRG sample is a rather 
unusual realisation, and the additional power just reflects sample variance. It is 
interesting to test the location of the cross-over scale in another red galaxy sample 
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at a different redshift. 



This paper is organised as follows. In Section [272] we introduce the 6dFGS sur- 

we 



vey and the i^-band selected sub-sample used in this analysis. In Section 2.3 



explain the technique we apply to derive the correlation function and summarise our 



error estimate, which is based on log-normal realisations. In Section 2.4 we discuss 
the need for wide angle corrections and several linear and non-linear effects which 
influence our measurement. Based on this discussion we introduce our correlation 
function model. In Section 12.51 we fit the data and derive the distance estimate 
Dy(z e g). In Section 2.6 we derive the Hubble constant and constraints on dark 



energy. In Section 2.7 we discuss the significance of the BAO detection of 6dFGS. 



In Section |2.8| we give a short overview of future all-sky surveys and their power to 
measure the Hubble constant. We conclude and summarise our results in Section [2791 



2.2 The 6dF galaxy survey 



2.2.1 Targets and Selection Function 



The galaxies used in this analysis were selected to K < 12.9 from the 2MASS 



Extended Source Catalog [2MASS XSC; Jarrett et al. 2000 and combined with 



redshift data from the 6dF Galaxy Survey [6dFGS; Jones et al. 2009 . The 6dF 
Galaxy Survey is a combined redshift and peculiar velocity survey covering nearly 
the entire southern sky with |6| < 10°. It was undertaken with the Six-Degree Field 
(6dF) multi-fibre instrument on the UK Schmidt Telescope from 2001 to 2006. The 
median redshift of the survey is z = 0.052 and the 25% : 50% : 75% percentile 



completeness values are 0.61 : 0.79 : 0.92. Papers by Jones et al. 2004, 2006, 2009 



describe 6dFGS in full detail, including comparisons between 6dFGS, 2dFGRS and 
SDSS. 

Galaxies were excluded from our sample if they resided in sky regions with 
completeness lower than 60 percent. After applying these cuts our sample contains 
75 117 galaxies. The selection function was derived by scaling the survey complete- 
ness as a function of magnitude to match the integrated on-sky completeness, using 
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Figure 2.1: Redshift distribution of the data (black solid line) and the random 
catalogue (black dashed line). The weighted distribution (using weights from eq. 2.6 ) 
is shifted to higher redshift and has increased shot noise but a smaller error due to 
sample variance (blue solid and dashed lines). 



mean galaxy counts. This method is the same adopted by Colless et al. 2001 for 



2dFGRS and is explained in Jones et al. 2006 in detail. The redshift of each object 
was checked visually and care was taken to exclude foreground Galactic sources. 
The derived completeness function was used in the real galaxy catalogue to weight 
each galaxy by its inverse completeness. The completeness function was also ap- 
plied to the mock galaxy catalogues to mimic the selection characteristics of the 
survey. Jones et al. (in preparation) describe the derivation of the 6dFGS selection 
function, and interested readers are referred to this paper for a more comprehensive 
treatment. 
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2.2.2 Survey volume 



We calculated the effective volume of the survey using the estimate of |Tegmark 

n(x)P 



1997 



eff 



d x 



(2.1) 



l + n(x)P _ 

where n(x) is the mean galaxy density at position x, determined from the data, and 
Pq is the characteristic power spectrum amplitude of the BAO signal. The parameter 
Po is crucial for the weighting scheme introduced later. We find that the value of 
Po = 40 000/i~ 3 Mpc 3 (corresponding to the value of the galaxy power spectrum at 
k ~ 0.06/i Mpc -1 in 6dFGS) minimises the error of the correlation function near the 
BAO peak. 

Using P = 40 000/i~ 3 Mpc 3 yields an effective volume of 0.08/i -3 Gpc 3 , while 
using instead Po — 10 000/i~ 3 Mpc 3 (corresponding to A; ~ 0A5h Mpc" 1 ) gives an 
effective volume of 0.045/?r 3 Gpc 3 . 

The volume of the 6dF Galaxy Survey is approximately as large as the vol- 
ume covered by the 2dF Galaxy Redshift Survey, with a sample density similar to 



SDSS-DR7 Abazajian et al. 2009 . Percival et al. 2010 reported successful BAO 



detections in several samples obtained from a combination of SDSS DR7, SDSS- 
LRG and 2dFGRS with effective volumes in the range 0.15 - 0.45/i~ 3 Gpc 3 (using 



Po = 10 000/i Mpc ), while the original detection by Eisenstein et al. 2005| used 
a sample with Kff = 0.38/i" 3 Gpc 3 (using P = 40 000/i" 3 Mpc 3 ). 



2.3 Clustering measurement 

We focus our analysis on the two-point correlation function. In the following sub- 
sections we introduce the technique used to estimate the correlation function and 
outline the method of log-normal realisations, which we employed to derive a co- 
variance matrix for our measurement. 
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2.3.1 Random catalogues 

To calculate the correlation function we need a random sample of galaxies which 
follows the same angular and redshift selection function as the 6dFGS sample. We 
base our random catalogue generation on the 6dFGS luminosity function of Jones 
et al. (in preparation), where we use random numbers to pick volume- weighted red- 
shifts and luminosity function- weighted absolute magnitudes. We then test whether 
the redshift-magnitude combination falls within the 6dFGS i^-band faint and bright 
apparent magnitude limits (8.75 < K < 12.9). 



Figure 2.1 shows the redshift distribution of the 6dFGS fT-selected sample (black 
solid line) compared to a random catalogue with the same number of galaxies (black 
dashed line). The random catalogue is a good description of the 6dFGS redshift 
distribution in both the weighted and unweighted case. 



2.3.2 The correlation function 

We turn the measured redshift into co-moving distance via 

D ^=k[m (2 ' 2) 

with 

E{z)=[^{l + zf + nf{l + zf 

+ ^ d (l + ^) 3 ( 1+ - fld ))] 1/2 , (2.3) 

where the curvature f2jj d is set to zero, the dark energy density is given by f2^ d = 
1 — and the equation of state for dark energy is w M = — 1. Because of the very 
low redshift of 6dFGS, our data are not very sensitive to w or any other higher 
dimensional parameter which influences the expansion history of the Universe. We 
will discuss this further in Section 12.5.31 

Now we measure the separation between all galaxy pairs in our survey and count 
the number of such pairs in each separation bin. We do this for the 6dFGS data 
catalogue, a random catalogue with the same selection function and a combination 
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Figure 2.2: The large scale correlation function of 6dFGS. The best fit model is 
shown by the black line with the best fit value of VL m h 2 = 0.138 ± 0.020. Models 
with different fl m h 2 are shown by the green line (VL m h 2 = 0.12) and the blue line 
(Q m h 2 = 0.15). The red dashed line is a linear CDM model with Vt^h 2 = (and 
fl m h 2 = 0.1), while all other models use the WMAP-7 best fit value of Vt b h 2 = 
0.02227 jKomatsu et al. , 2011 . The significance of the BAO detection in the black 
line relative to the red dashed line is 2. 4a (see Section 2.7). The error-bars at 



the data points are the diagonal elements of the covariance matrix derived using 
log-normal mock catalogues. 



of data-random pairs. We call the pair-separation distributions obtained from this 
analysis DD(s), RR(s) and DR(s), respectively. The binning is chosen to be from 
10/i -1 Mpc up to 190/i 1 Mpc, in 10/z 1 Mpc steps. In the analysis we used 30 
random catalogues with the same size as the data catalogue. The redshift correlation 



function itself is given by Landy & Szalay 1993 



3 



datal 



1 + 



DD(s) 
RR(s) 



n r 



DR(s) 
'RR(s) 



n d 



(2.4) 
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where the ratio n r jn^ is given by 



(2.5) 



and the sums go over all random (N r ) and data (Nd) galaxies. We use the inverse 



density weighting of Feldman, Kaiser & Peacock 11994] : 



a 



wAx) 



l + n(f)P ' 



(2.6) 



with Pq = 40 000/i 3 Mpc 3 and Cj being the inverse completeness weighting for 



6dFGS (see Section 2.2.1 and Jones et al., in preparation). This weighting is designed 
to minimise the error on the BAO measurement, and since our sample is strongly 
limited by sample variance on large scales this weighting results in a significant 
improvement to the analysis. The effect of the weighting on the redshift distribution 



is illustrated in Figure 2.1 



Other authors have used the so called ^-weighting which optimises the error 



over all scales by weighting each scale differently [e.g. Efstathiou, 1988; Loveday et 



al. , 1995 . In a magnitude limited sample there is a correlation between luminosity 



and redshift, which establishes a correlation between bias and redshift Zehavi et 



al. , 2005 . A scale-dependent weighting would imply a different effective redshift for 



each scale, causing a scale dependent bias. 



Finally we considered a luminosity dependent weighting as suggested by |Per 

However the same authors found that explicitly 



cival, Verde & Peacock 2004 



accounting for the luminosity- redshift relation has a negligible effect for 2dFGRS. 
We found that the effect to the 6dFGS correlation function is <C la for all bins. 



Hence the static weighting of eq. 2.6 is sufficient for our dataset. 
We also include an integral constraint correction in the form of 

fdata(s) = £data(s) + *C, 



(2.7) 



where ic is defined as 



ic 



Y. s RR(s)£model(s) 

E s RR(s) 



{2.1 
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Figure 2.3: Correlation matrix derived from a covariance matrix calculated from 
200 log-normal realisations. 



The function RR(s) is calculated from our mock catalogue and £ m odei(s) is a corre- 
lation function model. Since ic depends on the model of the correlation function we 
have to re-calculate it at each step during the fitting procedure. However we note 
that ic has no significant impact to the final result. 



Figure 2J2 shows the correlation function of 6dFGS at large scales. The BAO 
peak at ~ 105/i _1 Mpc is clearly visible. The plot includes model predictions of 



different cosmological parameter sets. We will discuss these models in Section 2.5.2 



2.3.3 Log-normal error estimate 



To obtain reliable error-bars for the correlation function we use log-normal realisa- 
tions Coles & Jones, 1991; Cole et al. , 2005; Kitaura et al. 2009 . In what follows 



we summarise the main steps, but refer the interested reader to Appendix [A] in 
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which we give a detailed explanation of how we generate the log-normal mock cat- 
alogues. In Appendix [B] we compare the log-normal errors with jack-knife estimates. 

Log-normal realisations of a galaxy survey are usually obtained by deriving a 
density field from a model power spectrum, P(k), assuming Gaussian fluctuations. 
This density field is then Poisson sampled, taking into account the window function 
and the total number of galaxies. The assumption that the input power spectrum 
has Gaussian fluctuations can only be used in a model for a density field with over- 
densities < 1. As soon as we start to deal with finite r.m.s. fluctuations, the 
Gaussian model assigns a non-zero probability to regions of negative density. A 
log-normal random field LN(x), can avoid this unphysical behaviour. It is obtained 
from a Gaussian field G(x) by 

LN(x) = exp[G(x) - a 2 G /2] - 1 (2.9) 

where o G is the variance of the field. LN(x) is positive-definite but approaches 
1 + G(x) whenever the perturbations are small (e.g. at large scales). Calculating 
the power spectrum of a Poisson sampled density field with such a distribution will 
reproduce the input power spectrum convolved with the window function. As an 
input power spectrum for the log-normal field we use 

P nl (k) = AP hn (k) exp{-(k/h) 2 ] (2.10) 

where A = b 2 (l + 2/3/3 + /3 2 /5) accounts for the linear bias and the linear redshift 
space distortions. P\\ n {k) is a linear model power spectrum in real space obtained 



from CAMB Lewis, Challinor & Lasenby, 2000 and P n \(k) is the non-linear power 



spectrum in redshift space. Comparing the model above with the 6dFGS data gives 
A = 4. The damping parameter fc* is set to k* = 0.33/i Mpc -1 , as found in 6dFGS 
(see fitting results later). How well this input model matches the 6dFGS data can 



be seen in Figure 2.9 



We produce 200 such realisations and calculate the correlation function for each 
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of them, deriving a covariance matrix 



n=l 

Here, £ n (sj) is the correlation function estimate at separation and the sum goes 
over all N log-normal realisations. The mean value is defined as 

1 - 

f(*) = ^X>(*). ( 2 - 12 ) 

n=l 

The case i = j gives the error (ignoring correlations between bins, of = Ca). In the 
following we will use this uncertainty in all diagrams, while the fitting procedures 
use the full covariance matrix. 

The distribution of recovered correlation functions includes the effects of sample 
variance and shot noise. Non-linearities are also approximately included since the 
distribution of over-densities is skewed. 



In Figure 2.3 we show the log-normal correlation matrix calculated from the 



covariance matrix. The correlation matrix is defined as 

m = -£^=> ( 2 - 13 ) 

where C is the covariance matrix (for a comparison to jack-knife errors see ap- 
pendix [b| . 



2.4 Modelling the BAO signal 

In this section we will discuss wide-angle effects and non-linearities. We also intro- 
duce a model for the large scale correlation function, which we later use to fit our 
data. 
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2.4.1 Wide angle formalism 



The model of linear redshift space distortions introduced by Kaiser 1987 is based 



on the plane parallel approximation. Earlier surveys such as SDSS and 2dFGRS are 
at sufficiently high redshift that the maximum opening angle between a galaxy pair 
remains small enough to ensure the plane parallel approximation is valid. However, 
the 6dF Galaxy Survey has a maximum opening angle of 180° and a lower mean 
redshift of z ~ 0.1 (for our weighted sample) and so it is necessary to test the 
validity of the plane parallel approximation. The wide angle description of redshift 



space distortions has been laid out in several papers Szalay et al. 1997; Szapudi 



2004 Matsubara 2004 Papai & Szapudi, 2008; Raccanelli et al. 2010 , which we 



summarise in Appendix [Cj 

We find that the wide-angle corrections have only a very minor effect on our 
sample. For our fiducial model we found a correction of A£ = 4 • 10 -4 in amplitude 
at s = lOOh- 1 Mpc and A£ = 4.5 • 10~ 4 at s = 200/i" 1 Mpc, (Figure |^4| in the 
appendix). This is much smaller than the error bars on these scales. Despite the 
small size of the effect, we nevertheless include all first order correction terms in 
our correlation function model. It is important to note that wide angle corrections 
affect the correlation function amplitude only and do not cause any shift in the 
BAO scale. The effect of the wide-angle correction on the unweighted sample is 
much greater and is already noticeable on scales of 20h~ 1 Mpc. Weighting to higher 
redshifts mitigates the effect because it reduces the average opening angle between 
galaxy pairs, by giving less weight to wide angle pairs (on average). 



2.4.2 Non-linear effects 

There are a number of non-linear effects which can potentially influence a measure- 
ment of the BAO signal. These include scale-dependent bias, the non-linear growth 
of structure on smaller scales, and redshift-space distortions. We discuss each of 
these in the context of our 6dFGS sample. 

As the Universe evolves, the acoustic signature in the correlation function is 
broadened by non-linear gravitational structure formation. Equivalently we can say 
that the higher harmonics in the power spectrum, which represent smaller scales, 



58 



CHAPTER 2. BARYON ACOUSTIC OSCILLATIONS AND THE LOCAL 

HUBBLE CONSTANT 



are erased Eisenstein, Seo & White, 2007 



The early Universe physics, which we discussed briefly in the introduction, is 
well understood and several authors have produced software packages (e.g. CMB- 
FAST and CAMB) and published fitting functions [e.g Eisenstein & Hu, |1998 to 
make predictions for the correlation function and power spectrum using thermody- 
namical models of the early Universe. These models already include the basic linear 
physics most relevant for the BAO peak. In our analysis we use the CAMB software 
package Lewis, Challinor & Lasenby, 2000 . The non- linear evolution of the power 
spectrum in CAMB is calculated using the halofit code Smith et al. , 2003 . This 



code is calibrated by n-body simulations and can describe non-linear effects in the 
shape of the matter power spectrum for pure CDM models to an accuracy of around 



5 - 10% Heitmann et al. , 2010 . However, it has previously been shown that this 
non-linear model is a poor description of the non-linear effects around the BAO 
peak Crocce & Scoccimarro, 2008 . We therefore decided to use the linear model 
output from CAMB and incorporate the non-linear effects separately. 

All non-linear effects influencing the correlation function can be approximated 
by a convolution with a Gaussian damping factor exp[— (rA^/2) 2 ] jEisenstein et al 



2007 Eisenstein, Seo & White, 2007 , where k* is the damping scale. We will use 
this factor in our correlation function model introduced in the next section. The 
convolution with a Gaussian causes a shift of the peak position to larger scales, 
since the correlation function is not symmetric around the peak. However this shift 
is usually very small. 

All of the non-linear processes discussed so far are not at the fundamental scale 
of 105/i -1 Mpc but are instead at the cluster-formation scale of up to lO/i^Mpc. The 
scale of \05h~ 1 Mpc is far larger than any known non-linear effect in cosmology. This 
has led some authors to the conclusion that the peak will not be shifted significantly, 



but rather only blurred out. For example, Eisenstein, Seo & White 2007 have 
argued that any systematic shift of the acoustic scale in real space must be small 
(< 0.5%), even at z = 0. 



However, several authors report possible shifts of up to 1% Guzik & Bernstein 



2007 


Smith et al. 


2008 


2007 



used re-normalised perturbation theory (RPT) and found percent-level shifts in the 
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BAO peak. In addition to non-linear evolution, they found that mode-coupling 
generates additional oscillations in the power spectrum, which are out of phase with 
the BAO oscillations predicted by linear theory. This leads to shifts in the scale of 
oscillation nodes with respect to a smooth spectrum. In real space this corresponds 
to a peak shift towards smaller scales. Based on their results, [Crocce fc Scoccimarro 



1 2008 propose a model to be used for the correlation function analysis at large scales. 
We will introduce this model in the next section. 



2.4.3 Large-scale correlation function 



To model the correlation function on large scales, we follow [Crocce fc Scoccimarro 
2008 1 and Sanchez et al. 2008 and adopt the following parametrisatiorj^J 



model v 



B(s)b 2 



««)*G(r) + tf(r) 



ggg) 
ds 



(2.14) 



Here, we decouple the scale dependency of the bias B(s) and the linear bias b. G{r) is 
a Gaussian damping term, accounting for non-linear suppression of the BAO signal. 
£(s) is the linear correlation function (including wide angle description of redshift 
space distortions; eq. 5.8 in the appendix). The second term in eq 2.14 accounts 



for the mode-coupling of different Fourier modes. It contains d£,(s)/ds, which is the 
first derivative of the redshift space correlation function, and £i(r), which is defined 

as 

1 f°° 

£i(r) = — 2 / dk kP hn (k)ji(rk), (2.15) 



with ji(x) being the spherical Bessel function of order 1. Sanchez et al. 2008 used an 
additional parameter Amc which multiplies the mode coupling term in equation|2.14 



We found that our data are not good enough to constrain this parameter, and hence 
adopted A MC = 1 as in the original model by Crocce & Scoccimarro 2008 . 



In practice we generate linear model power spectra Pn n {k) from CAMB and 



1 note that r — s, the different letters just specify whether the function is evaluated in redshift 
space or real space. 
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convert them into a correlation function using a Hankel transform 

1 f°° 

£( r ) = / dk k 2 PUk) Jo (rk) } (2.16) 
^ Jo 

where jo(x) = sin(x)/x is the spherical Bessel function of order 0. 

The *-symbol in eq. 2.14| is only equivalent to a convolution in the case of a 



3D correlation function, where we have the Fourier theorem relating the 3D power 
spectrum to the correlation function. In case of the spherically averaged quantities 
this is not true. Hence, the *-symbol in our equation stands for the multiplication 
of the power spectrum with G(k) before transforming it into a correlation function. 
G{k) is defined as 

G{k) = exp [-(k/k*) 2 ] , (2.17) 

with the property 

G(k) -> as k ->■ oo. (2.18) 



The damping scale k* can be calculated from linear theory Crocce & Scoccimarro 



2006 Matsubara[ 2008 by 

1 



k-M 



dk P hn (k) 



1-1/2 







(2.19) 



where P\i n (k) is again the linear power spectrum. ACDM predicts a value of k* ~ 
0.17/i Mpc -1 . However, we will include fc* as a free fitting parameter. 

The scale dependance of the 6dFGS bias, B(s), is derived from the GiggleZ 
simulation (Poole et al., in preparation); a dark matter simulation containing 2160 3 
particles in a lh~ l Gpc box. We rank-order the halos of this simulation by Vmax 
and choose a contiguous set of 250 000 of them, selected to have the same clustering 
amplitude of 6dFGS as quantified by the separation scale r , where C( r o) — 1- 111 
the case of 6dFGS we found r = 9.3/i -1 Mpc. Using the redshift space correlation 
function of these halos and of a randomly subsampled set of ~ 10 6 dark matter 
particles, we obtain 



B(s) = 1 + (s/O.mh^M-pc) L332 , (2.20) 
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Table 2.1: This table contains all parameter constraints from 6dFGS obtained in 
this paper. The priors used to derive these parameters are listed in square brackets. 
All parameters assume VL^Ii 2 = 0.02227 and in cases where a prior on Vt m h 2 is 



used, we adopt the WMAP-7 Markov chain probability distribution Komatsu et 
al. 2011] . A(z e s) is the acoustic parameter defined by Eisenstein et al. 2005| (see 



equation 2.27 in the text) and R(z e g) is the distance ratio of the 6dFGS BAO 
measurement to the last-scattering surface. The most sensible value for cosmological 
parameter constraints is r s (zd) / Dv(z e s), since this measurement is uncorrelated with 
VL m h 2 . The effective redshift of 6dFGS is z e ff = 0.106 and the fitting range is from 
10 - 190/i- 1 Mpc. 

Summary of parameter constraints from 6dFGS 



[VL m h 2 prior] 



Q m h 2 0.138 ±0.020 (14.5%) 

D v (z eS ) 456 ± 27 Mpc (5.9%) 

D v (z eS ) 459 ± 18 Mpc (3.9%) 

r s (z d )/D v (z eff ) 0.336 ±0.015 (4.5%) 

R{zJ) 0.0324 ± 0.0015 (4.6%) 

A(z cS ) 0.526 ± 0.028 (5.3%) 

n m 0.296 ± 0.028 (9.5%) \Vt m h 2 prior] 

H 67.0 ± 3.2 (4.8%) [fi m h 2 prior] 



which describes a 1.7% correction of the correlation function amplitude at separation 
scales of 10/i -1 Mpc. To derive this function, the GiggleZ correlation function 
(snapshot z = 0) has been fitted down to 6/i _1 Mpc, well below the smallest scales 
we are interested in. 



2.5 Extracting the BAO signal 

In this section we fit the model correlation function developed in the previous section 
to our data. Such a fit can be used to derive the distance scale Dv(z e s) at the 
effective redshift of the survey. 
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2.5.1 Fitting preparation 

The effective redshift of our sample is determined by 



^eff 



N b N b 

2Nj l 



Zi + z. 



■3)1 



(2.21) 



where iV& is the number of galaxies in a particular separation bin and Wi and Wj 



are the weights for those galaxies from eq. 2J3 We choose z e ff from bin 10 which 
has the limits 100/i -1 Mpc and llO/i" 1 Mpc and which gave z e g = 0.106. Other 
bins show values very similar to this, with a standard deviation of ±0.001. The 
final result does not depend on a very precise determination of z e g, since we are not 
constraining a distance to the mean redshift, but a distance ratio (see equation 



2.24 



later). In fact, if the fiducial model is correct, the result is completely independent 
of z e s- Only if there is a ^-dependent deviation from the fiducial model do we need 
z e s to quantify this deviation at a specific redshift. 



Along the line-of-sight, the BAO signal directly constrains the Hubble constant 
H(z) at redshift z. When measured in a redshift shell, it constrains the angular 
diameter distance D^(z) |Matsubara 2004 . In order to separately measure Da (z) 
and H(z) we require a BAO detection in the 2D correlation function, where it will 
appear as a ring at around 105/ir 1 Mpc. Extremely large volumes are necessary 
for such a measurement. While there are studies that report a successful (but very 
low signal-to-noise) detection in the 2D correlation function using the SDSS-LRG 



Gaztanaga, Cabre & Hui 


2009 


Chuang & Wang 


2011 



et al.|2010 , our sample does not allow this kind of analysis. Hence we restrict ourself 
to the ID correlation function, where we measure a combination of Da{z) and H(z). 
What we actually measure is a superposition of two angular measurements (R.A. 
and Dec.) and one line-of-sight measurement (redshift). To account for this mixture 



of measurements it is common to report the BAO distance constraints as Eisenstein 



et al., 2005; Padmanabhan & White, 2008 



D v (z) 



:i + z?d) 



cz 



1/3 



H E(z) 



(2.22) 
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where Da is the angular distance, which in the case of Qf. = is given by Da{z) = 
D c (z)/(l + z). 

To derive model power spectra from CAMB we have to specify a complete 
cosmological model, which in the case of the simplest ACDM model (Qk = 0, w — 
— 1), is specified by six parameters: oj c ,oJb,n s ,T, A s and h. These parameters are: 
the physical cold dark matter and baryon density, (u c = Q c h 2 , Ub = Qth 2 ), the scalar 
spectral index, (n s ), the optical depth at recombination, (r), the scalar amplitude 
of the CMB temperature fluctuation, (A s ), and the Hubble constant in units of 
100 km s^Mpc" 1 (h). 

Our fit uses the parameter values from WMAP-7 Komatsu et al. 2011 : Vl^h 2 = 



0.02227, r = 0.085 and n s = 0.966 (maximum likelihood values). The scalar am- 
plitude A s is set so that it results in as = 0.8, which depends on Vt m h 2 . However 
erg is degenerated with the bias parameter b which is a free parameter in our fit. 
Furthermore, h is set to 0.7 in the fiducial model, but can vary freely in our fit 
through a scale distortion parameter a, which enters the model as 

£modcl(s) = Codcl(" S )- (2-23) 

This parameter accounts for deviations from the fiducial cosmological model, which 



we use to derive distances from the measured redshift. It is defined as Eisenstein 



et al., 2005 Padmanabhan & White, 2008 



Dy(zes) 
D^ d (z eS ) 



The parameter a enables us to fit the correlation function derived with the fidu- 
cial model, without the need to re-calculate the correlation function for every new 
cosmological parameter set. 

At low redshift we can approximate H(z) ~ H , which results in 

a « -3-. (2.25) 



Compared to the correct equation 2.24 this approximation has an error of about 3% 
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Figure 2.4: Likelihood contours of the distance Dy(z e s) against Q m h 2 . The corre- 
sponding values of a are given on the right-hand axis. The contours show 1 and 
2a errors for both a full fit (blue solid contours) and a fit over 20 — 190/i -1 Mpc 
(black dashed contours) excluding the first data point. The black cross marks the 
best fitting values corresponding to the dashed black contours with (D v ,Q m h 2 ) = 
(462, 0.129), while the blue cross marks the best fitting values for the blue contours. 
The black solid curve corresponds to a constant fl m h 2 Dv(z e s) (Dy ~ h~ x ), while the 
dashed line corresponds to a constant angular size of the sound horizon, as described 
in the text. 



at redshift z = 0.1 for our fiducial model. Since this is a significant systematic bias, 
we do not use this approximation at any point in our analysis. 

2.5.2 Extracting Dy(z e s) and r s (zd) / Dy(z e s) 

Using the model introduced above we performed fits to 18 data points between 
lOh^ 1 Mpc and 190/i -1 Mpc. We excluded the data below 10/i -1 Mpc, since our 
model for non-linearities is not good enough to capture the effects on such scales. 
The upper limit is chosen to be well above the BAO scale, although the constraining 
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Figure 2.5: The distance measurement Dy{z) relative to a low redshift approxima- 



tion. The points show 6dFGS data and those of Percival et al. 2010 . 



contribution of the bins above 130/i 1 Mpc is very small. Our final model has 4 free 
parameters: fl m h 2 , b, a and k*. 

The best fit corresponds to a minimum \ 2 of 15.7 with 14 degrees of freedom (18 



data-points and 4 free parameters). The best fitting model is included in Figure 2.2 
(black line). The parameter values are fl m h 2 = 0.138 ± 0.020, b = 1.81 ± 0.13 and 
a = 1.036 ± 0.062, where the errors are derived for each parameter by marginalising 
over all other parameters. For we can give a lower limit of fc* = 0.19/i Mpc -1 
(with 95% confidence level). 



We can use eq. |2.24| to turn the measurement of a into a measurement of the 
distance to the effective redshift Dv(z e s) = aD^(z c s) = 456 ± 27 Mpc, with a 
precision of 5.9%. Our fiducial model gives D^(z e g) = 440.5 Mpc, where we have 



followed the distance definitions of Wright 2006) throughout. For each fit we de- 

■ u 545 /b, which we need to calculate the wide angle 



rive the parameter (3 = il m ( : ) 
corrections for the correlation function. 

The maximum likelihood distribution of k* seems to prefer smaller values than 
predicted by ACDM, although we are not able to constrain this parameter very well. 
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This is connected to the high significance of the BAO peak in the 6dFGS data (see 
Section 2.7). A smaller value of k* damps the BAO peak and weakens the distance 
constraint. For comparison we also performed a fit fixing k* to the ACDM prediction 
of k* ~ 0.17/i Mpc -1 . We found that the error on the distance Dy(z c s) increases 
from 5.9% to 8%. However since the data do not seem to support such a small value 
of k* we prefer to marginalise over this parameter. 

The contours of Dy(z e e) — Q m h 2 are shown in Figure 2.4 together with two 
degeneracy predictions jEisenstein et al. 2005 . The solid line is that of constant 
Q m h 2 Dy(z e f[), which gives the direction of degeneracy for a pure CDM model, where 
only the shape of the correlation function contributes to the fit, without a BAO peak. 
The dashed line corresponds to a constant r s (zd)/Dy(z e e), which is the degeneracy 
if only the position of the acoustic scale contributes to the fit. The dashed contours 
exclude the first data point, fitting from 20 - 190/i -1 Mpc only, with the best fitting 
values a = 1.049 ± 0.071 (corresponding to D v (z eS ) = 462 ± 31 Mpc), Vt m h 2 = 
0.129 ± 0.025 and b = 1.72 ± 0.17. The contours of this fit are tilted towards 
the dashed line, which means that the fit is now driven by the BAO peak, while 
the general fit (solid contours) seems to have some contribution from the shape of 
the correlation function. Excluding the first data point increases the error on the 
distance constraint only slightly from 5.9% to 6.8%. The value of VL m h 2 tends to be 
smaller, but agrees within la with the former value. 

Going back to the complete fit from 10 - 190/i -1 Mpc, we can include an external 
prior on VL m h 2 from WMAP-7, which carries an error of only 4% (compared to the 
~ 15% we obtain by fitting our data). Marginalising over VL m h 2 now gives Dy(z e e) = 
459 ± 18 Mpc, which reduces the error from 5.9% to 3.9%. The uncertainty in Q m h 2 
from WMAP-7 contributes only about 5% of the error in Dy (assuming no error in 
the WMAP-7 value of Vt m h 2 results in D v (z eS ) = 459 ± 17 Mpc). 

In Figure 2.5 we plot the ratio D v (z) / D l ° w ~ z (z) as a function of redshift, where 



jjiow _ cz jjj Q At sufficiently low redshift the approximation H(z) ~ Hq is 
valid and the measurement is independent of any cosmological parameter except the 



Hubble constant. This figure also contains the results from Percival et al. 2010 . 

Rather than including the WMAP-7 prior on Q m h 2 to break the degeneracy 
between Q m h 2 and the distance constraint, we can fit the ratio r s (zd)/Dy(z e s), 
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where r s (zd) is the sound horizon at the baryon drag epoch Zd- In principle, this is 
rotating Figure 2A so that the dashed black line is parallel to the x-axis and hence 
breaks the degeneracy if the fit is driven by the BAO peak; it will be less efficient if 
the fit is driven by the shape of the correlation function. During the fit we calculate 
r s (zd) using the fitting formula of Eisenstein & Hu 1998 . 



The best fit results in r s (zd)/Dy(z c s) = 0.336 ± 0.015, which has an error of 
4.5%, smaller than the 5.9% found for Dy but larger than the error in Dy when 
adding the WMAP-7 prior on Q m h 2 . This is caused by the small disagreement in 



the Dy — Q m h 2 degeneracy and the line of constant sound horizon in Figure 2.4 The 
X 2 is 15.7, similar to the previous fit with the same number of degrees of freedom. 



2.5.3 Extracting A(z c g) and R(z e s) 



We can also fit for the ratio of the distance between the effective redshift, z e g, and 



the redshift of decoupling [z* = 1091; Eisenstein et al. , 2005 



R{z, 



Dv{z eS ) 



(i + z*)D A {z*y 



(2.26) 



with (1 + z*)D A (z*) being the CMB angular comoving distance. Beside the fact 
that the Hubble constant Hq cancels out in the determination of R, this ratio is also 



more robust against effects caused by possible extra relativistic species Eisenstein & 
White, 2004 . We calculate D A (z*) for each Q m h 2 during the fit and then marginalise 
over tt m h 2 . The best fit results in R = 0.0324 ± 0.0015, with \ 2 = 15.7 and the 
same 14 degrees of freedom. 

Focusing on the path from z = to z e g = 0.106, our dataset can give interesting 



constraints on Q m . We derive the parameter Eisenstein et al. , 2005 



A(Z CS ) = lOOZMZeff 



Vtt m h 2 

cz eS 



(2.27) 



which has no dependence on the Hubble constant since Dy oc h 1 . We obtain 
A(z cS ) = 0.526 ± 0.028 with x 2 /d.o.f . = 15.7/14. The value of A would be identical 
to \/T2 m if measured at redshift z — 0. At redshift z c s = 0.106 we obtain a deviation 
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from this approximation of 6% for our fiducial model, which is small but systematic. 
We can express A, including the curvature term Q k and the dark energy equation 
of state parameter w, as 



A(z) 



/a 



E(z) 1 / 3 



sinh( v / n^x(2 ; )) 
1 2 /3 



2/3 



2/3 



'\n k \z 



n k >o 

Q k = 
tt k < 



(2.28) 



with 



and 



X (z) = D c {z) 



Ho 



dz' 
Eiz 7 ' 



(2.29) 



E{z) = [n m {l + zf + ft fc (l + 



+ n A (i + z) 



3(l+w)\] i/ 2 



(2.30) 



Using this equation we now linearise our result for Q m in Q k and w and get 



tt m = 0.287 + 0.039(1 + w) + 0.039fi fc ± 0.027. 



(2.31) 



For comparison, Eisenstein et al. 2005 found 



Q m = 0.273 + 0.123(1 + w) + 0.137fi fc ± 0.025 



(2.32) 



based on the SDSS LRG DR3 sample. This result shows the reduced sensitivity of 
the 6dFGS measurement to w and fl k - 



2.6 Cosmo logical implications 

In this section we compare our results to other studies and discuss the implications 
for constraints on cosmo logical parameters. We first note that we do not see any 
excess correlation on large scales as found in the SDSS-LRG sample. Our correla- 
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Table 2.2: wCDM constraints from different datasets. Comparing the two columns 
shows the influence of the 6dFGS data point. The 6dFGS data point reduces the 
error on w by 24% compared to WMAP-7+LRG which contains only the BAO data 
points of Percival et al. 2010 . We assume flat priors of 0.11 < Q m h 2 < 0.16 and 



marginalise over VL m h . The asterisks denote the free parameters in each fit. 



parameter 


WMAP-7+LRG 


WMAP-7+LRG+6dFGS 




69.9±3.8(*) 


68.7 ± 2. 8(*) 


Q m 


0.283 ± 0.033 


0.293 ±0.027 


Q A 


0.717 ±0.033 


0.707 ±0.027 


w 


-1.01 ±0.17(*) 


-0.97±0.13(*) 



Table 2.3: Parameter constraints from WMAP7±BAO for (i) a flat ACDM model, 
(ii) an open ACDM (oACDM), (iii) a flat model with w = const. (wCDM), and 
(iv) an open model with w = constant (owCDM). We assume flat priors of 0.11 < 
VL m h 2 < 0.16 and marginalise over VL m h 2 . The asterisks denote the free parameters 
in each fit. 



parameter 


ACDM 


oACDM 


wCDM 


owCDM 




69.2 ± 1.1(*) 


68.3 ± 1.7(*) 


68.7±2.8(*) 


70.4±4.3(*) 




0.288 ±0.011 


0.290 ±0.019 


0.293 ±0.027 


0.274 ± 0.035 




(0) 


-0.0036 ±0.0060(*) 


(0) 


-0.013 ±0.010(*) 




0.712 ±0.011 


0.714 ±0.020 


0.707 ±0.027 


0.726 ±0.036 


w 


(-1) 


(-1) 


-0.97±0.13(*) 


-1.24±0.39(*) 



tion function is in agreement with a crossover to negative scales at lAOh Mpc, as 
predicted from ACDM. 



2.6.1 Constraining the Hubble constant, H { 



o 



We now use the 6dFGS data to derive an estimate of the Hubble constant. We use 
the 6dFGS measurement of r s (zd)/Dy (0.106) = 0.336 ± 0.015 and fit directly for 
the Hubble constant and Q m . We combine our measurement with a prior on Q m h 2 



coming from the WMAP-7 Markov chain results Komatsu et al. , 2011 . Combining 
the clustering measurement with VL m h 2 from the CMB corresponds to the calibration 



70 



CHAPTER 2. BARYON ACOUSTIC OSCILLATIONS AND THE LOCAL 

HUBBLE CONSTANT 




Figure 2.6: The blue contours show the WMAP-7 Q m h 2 prior Komatsu et al. , 2011 



The black contour shows constraints from 6dFGS derived by fitting to the measure- 
ment of r s (zd)/Dy(z e fi). The solid red contours show the combined constraints 
resulting in H = 67.0 ± 3.2 km s^Mpc" 1 and fi m = 0.296 ± 0.028. Combining the 
clustering measurement with Q m h 2 from the CMB corresponds to the calibration of 
the standard ruler. 



of the standard ruler. 

We obtain values of Hq = 67.0 ± 3.2 km s _1 Mpc _1 (which has an uncertainty 



of only 4.8%) and Q m = 0.296 ± 0.028. Table 2.1 and Figure 2.6 summarise the 
results. The value of Q m agrees with the value we derived earlier (Section 2.5.3). 

To combine our measurement with the latest CMB data we use the WMAP-7 
distance priors, namely the acoustic scale 



[I + z.) 



r s (z*) 



(2.33) 
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Figure 2.7: The blue contours shows the WMAP-7 degeneracy in H and w Komatsu 



et al. , 2011 , highlighting the need for a second dataset to break the degeneracy. The 
black contours show constraints from BAO data incorporating the r s (zd)/ D v (z e s) 

and 6dFGS. The solid red contours show the 
= -0.97 ± 0.13. Excluding the 6dFGS data 

-1.01 ±0.17. 



measurements of Percival et al. 2010 



combined constraints resulting in w = 
point widens the constraints to the dashed red line with w 



the shift parameter 



R = 100 ^ Qmh2 n + z AD A (z' 



(2.34) 



and the redshift of decoupling z* (Tables 9 and 10 in Komatsu et al. 2011[ ). This 
combined analysis reduces the error further and yields H = 68.7±1.5 km s _1 Mpc _1 
(2.2%) and VL m = 0.29 ± 0.022 (7.6%). 

Percival et al. 2010 determine a value of H = 68.6 ± 2.2 km s~ 1 Mpc~ 1 using 
SDSS-DR7, SDSS-LRG and 2dFGRS, while |Reid efaL] [20lo] found H = 69.4 ± 
1.6 km s _1 Mpc _1 using the SDSS-LRG sample and WMAP-5. In contrast to these 
results, 6dFGS is less affected by parameters like Qk and w because of its lower 
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redshift. In any case, our result of the Hubble constant agrees very well with earlier 
BAO analyses. Furthermore our result agrees with the latest CMB measurement of 
Hq = 70.3 ± 2.5 km s^Mpc" 1 



Komatsu et al. , 2011 



The SHOES program Riess et al. , 2011 determined the Hubble constant using 
the distance ladder method. They used about 600 near-IR observations of Cepheids 
in eight galaxies to improve the calibration of 240 low redshift (z < 0.1) SN la, and 
calibrated the Cepheid distances using the geometric distance to the maser galaxy 
NGC 4258. They found H = 73.8 ± 2.4 km s^Mpc -1 , a value consistent with the 
initial results of the Hubble Key project Freedman et al. [Hq = 72 ±8 km s~ 1 Mpc~ 1 ; 



2000 but 1.7a higher than our value (and 1.8c higher when we combine our dataset 



with WMAP-7). While this could point toward unaccounted or under-estimated 
systematic errors in either one of the methods, the likelihood of such a deviation by 
chance is about 10% and hence is not enough to represent a significant discrepancy. 
Possible systematic errors affecting the BAO measurements are the modelling of 
non-linearities, bias and redshift-space distortions, although these systematics are 
not expected to be significant at the large scales relevant to our analysis. 

To summarise the finding of this section we can state that our measurement 
of the Hubble constant is competitive with the latest result of the distance ladder 
method. The different techniques employed to derive these results have very differ- 
ent potential systematic errors. Furthermore we found that BAO studies provide 
the most accurate measurement of Hq that exists, when combined with the CMB 
distance priors. 



2.6.2 Constraining dark energy 

One key problem driving current cosmology is the determination of the dark energy 
equation of state parameter, w. When adding additional parameters like w to ACDM 
we find large degeneracies in the WMAP-7-only data. One example is shown in 



Figure 2.7 WMAP-7 alone can not constrain Hq or w within sensible physical 
boundaries (e.g. w < —1/3). As we are sensitive to H , we can break the degeneracy 
between w and Hq inherent in the CMB-only data. Our assumption of a fiducial 
cosmology with w — — 1 does not introduce a bias, since our data are not sensitive 
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to this parameter and any deviation from this assumption is modelled within the 
shift parameter a. 

We again use the WMAP-7 distance priors introduced in the last section. In ad- 



dition to our value of r s (zd) / Dy(0. 106) = 0.336±0.015 we use the results of Percival 



et al. 2010 , who found r s (z d )/D v (0.2) = 0.1905 ± 0.0061 and r s (z d )/D v (0.35) 



0.1097±0.0036. To account for the correlation between the two latter data points we 
employ the covariance matrix reported in their paper. Our fit has 3 free parameters, 
fl m h 2 , H and w. 

The best fit gives w = -0.97 ±0.13, H = 68.7 ±2.8 km s^Mpc" 1 and Vt m h 2 = 
0.1380 ± 0.0055, with a x 2 /d.o.f. = L3 A Table [2~2] and Figure [2T7] summarise the 
results. To illustrate the importance of the 6dFGS result to the overall fit we also 
show how the results change if 6dFGS is omitted. The 6dFGS data improve the 
constraint on w by 24%. 

Finally we show the best fitting cosmological parameters for different cosmolog- 



ical models using WMAP-7 and BAO results in Table 2.3 



2.7 Significance of the BAO detection 

To test the significance of our detection of the BAO signature we follow |Eisenstem 



et al. 2005 and perform a fit with a fixed ilb = 0, which corresponds to a pure 
CDM model without a BAO signature. The best fit has \ 2 = 21-4 with 14 degrees 



of freedom and is shown as the red dashed line in Figure 2^ The parameter values 
of this fit depend on the parameter priors, which we set to 0.7 < a < 1.3 and 
0.1 < Vt m h 2 < 0.2. Values of a much further away from 1 are problematic since 



eq. 2.24 is only valid for a close to 1. Comparing the best pure CDM model with 
our previous fit, we estimate that the BAO signal is detected with a significance 
of 2.4er (corresponding to A% 2 = 5.6). As a more qualitative argument for the 



detection of the BAO signal we would like to refer again to Figure 2A_ where the 
direction of the degeneracy clearly indicates the sensitivity to the BAO peak. 

We can also use the log-normal realisations to determine how likely it is to find 
a BAO detection in a survey like 6dFGS. To do this, we produced 200 log-normal 
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Figure 2.8: The number of log-normal realisations found with a certain a/ A^ 2 , 
where the Ax 2 is obtained by comparing a fit using a ACDM correlation function 
model with a no-baryon model. The blue line indicates the 6dFGS result. 



mock catalogues and calculated the correlation function for each of them. We can 
now fit our correlation function model to these realisations. Furthermore, we fit a 
no-baryon model to the correlation function and calculate A% 2 , the distribution of 



which is shown in Figure 2^8. We find that 26% of all realisations have at least a 2a 
BAO detection, and that 12% have a detection > 2 Act. The log-normal realisations 
show a mean significance of the BAO detection of 1.7±0.7er, where the error describes 
the variance around the mean. 



Figure 2.9 shows the 6dFGS data points together with all 200 log-normal real- 
isations (grey). The red data points indicate the mean for each bin and the black 



line is the input model derived as explained in Section 2.3.3 This comparison shows 
that the 6dFGS data contain a BAO peak slightly larger than expected in ACDM. 

The amplitude of the acoustic feature relative to the overall normalisation 
of the galaxy correlation function is quite sensitive to the baryon fraction, /& = 
Qb/Q m |Matsubara, 2004 . A higher BAO peak could hence point towards a larger 
baryon fraction in the local Universe. However since the correlation function model 



75 



2. 7. SIGNIFICANCE OF THE BAO DETECTION 



250 




) r i i i I i i i I i i i I i i i I i i i I i i i I i i : I I : i I I I, I 

20 40 60 80 100 120 140 160 180 200 
S [IT 1 Mpc] 



Figure 2.9: The different log-normal realisations used to calculate the covariance 
matrix (shown in grey). The red points indicate the mean values, while the blue 
points show actual 6dFGS data (the data point at 5/i _1 Mpc is not included in the 
fit). The red data points are shifted by 2h~ x Mpc to the right for clarity. 



seems to agree very well with the data (with a reduced \ 2 °f 1-12) and is within 
the range spanned by our log-normal realisations, we can not claim any discrepancy 
with ACDM. Therefore, the most likely explanation for the excess correlation in the 
BAO peak is sample variance. 



In Figure 2.10 we show the distribution of the parameter a obtained from the 
200 log-normal realisations. The distribution is well described by a Gaussian with 
X 2 /d.o.f. = 14.2/20, where we employed Poisson errors for each bin. This confirms 
that a has Gaussian distributed errors in the approximation that the 6dFGS sample 
is well-described by log-normal realisations of an underlying ACDM power spectrum. 
This result increases our confidence that the application of Gaussian errors for the 
cosmological parameter fits is correct. The mean of the Gaussian distribution is at 
0.998 ± 0.004 in agreement with unity, which shows, that we are able to recover the 
input model. The width of the distribution shows the mean expected error in a in 
a ACDM universe for a 6dFGS-like survey. We found a = 0.057 ± 0.005 which is in 
agreement with our error in a of 5.9%. Figure 2.10 also contains the distribution 
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Figure 2.10: This plot shows the distribution of the parameter a derived from the 
200 log- normal realisations (black). The distribution is well fit by a Gaussian with 
a mean of \x = 0.998 ± 0.004 and a width of a = 0.057 ± 0.005. In blue we show 
the same distribution selecting only the log-normal realisations with a strong BAO 
peak (> 2a). The Gaussian distribution in this case gives a mean of 1.007 ± 0.007 
and a = 0.041 ± 0.008. 

of a, selecting only the log- normal realisations with a strong (> 2a) BAO peak 
(blue data). We included this selection to show, that a stronger BAO peak does 
not bias the estimate of a in any direction. The Gaussian fit gives x 2 /d.o.f. = 5/11 
with a mean of 1.007 ± 0.007. The distribution of a shows a smaller spread with 
a = 0.041 ± 0.008, about 2a below our error on a. This result shows, that a survey 
like 6dFGS is able to constrain a (and hence Dy and Hq) to the precision we report 
in this paper. 

2.8 Future all sky surveys 

A major new wide-sky survey of the local Universe will be the Wide field ASKAP L- 
band Legacy All-sky Blind surveY (WALLABY)[^J This is a blind HI survey planned 

2 http://www. atnf.csiro.au/research/ WALLABY 
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Figure 2.11: Redshift distribution of 6dFGS, WALLABY and two different versions 
of the proposed TAIPAN survey. See text for details. 



for the Australian SKA Pathfinder telescope (ASKAP), currently under construction 
at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. 

The survey will cover at least 75% of the sky with the potential to cover 47r of 
sky if the Westerbork Radio Telescope delivers complementary northern coverage. 
Compared to 6dFGS, WALLABY will more than double the sky coverage including 
the Galactic plane. WALLABY will contain ~ 500 000 to 600 000 galaxies with 
a mean redshift of around 0.04, giving it around 4 times greater galaxy density 
compared to 6dFGS. In the calculations that follow, we assume for WALLABY a 
Att survey without any exclusion around the Galactic plane. The effective volume 
in this case turns out to be 0.12/i -3 Gpc 3 . 

The TAIPAN survej]^] proposed for the UK Schmidt Telescope at Siding Spring 
Observatory, will cover a comparable area of sky, and will extend 6dFGS in both 
depth and redshift (z ~ 0.08). 



The redshift distribution of both surveys is shown in Figure 2.11, alongside 



3 TAIPAN: Transforming Astronomical Imaging surveys through Polychromatic Analysis of Neb- 
ulae 
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Figure 2.12: Predictions for two versions of the proposed TAIPAN survey. Both 
predictions assume a 2n steradian southern sky-coverage, excluding the Galactic 
plane (i.e. \b\ > 10°). TAIPAN 1 contains 406 000 galaxies while TAIPAN2 contains 



221000, (see Figure 2.11). The blue points are shifted by 2h~ l Mpc to the right for 
clarity. The black line is the input model, which is a ACDM model with a bias of 1.6, 
(3 = 0.3 and k* = 0.17h Mpc -1 . For a large number of realisations, the difference 
between the input model and the mean (the data points) is only the convolution 
with the window function. 



6dFGS. Since the TAIPAN survey is still in the early planning stage we consider 
two realisations: TAIPAN1 (406 000 galaxies to a faint magnitude limit of r = 17) 
and the shallower TAIPAN2 (221000 galaxies to r = 16.5). We have adopted the 
same survey window as was used for 6dFGS, meaning that it covers the whole 
southern sky excluding a 10° strip around the Galactic plane. The effective volumes 
of TAIPAN 1 and TAIPAN2 are 0.23/i~ 3 Gpc 3 and 0.13/z 3 Gpc 3 , respectively. 

To predict the ability of these surveys to measure the large scale correlation 
function we produced 100 log-normal realisations for TAIPAN1 and WALLABY 



and 200 log-normal realisations for TAIPAN2. Figures 2.12 and 2.13 show the 
results in each case. The data points are the mean of the different realisations, and 
the error bars are the diagonal of the covariance matrix. The black line represents 
the input model which is a ACDM prediction convolved with a Gaussian damping 
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Figure 2.13: Prediction for the WALLABY survey. We have assumed a An steradian 
survey with 602 000 galaxies, b = 0.7, (3 = 0.7 and K = 0.17h Mpc" 1 . 



term using = 0.17/i Mpc 1 (see eq. 2.17). We used a bias parameter of 1.6 



for TAIPAN and following our fiducial model we get /3 = 0.3, resulting in A = 
b 2 (l + 2/3/3 + /3 2 /5) = 3.1. For WALLABY we used a bias of 0.7 [based on the 
results found in the HIPASS survey; Basilakos et al. , 2007 . This results in (3 = 0.7 
and A = 0.76. To calculate the correlation function we used = 40 000/i 3 Mpc 3 
for TAIPAN and P = 5 000/i 3 Mpc 3 for WALLABY. 

The error bar for TAIPAN1 is smaller by roughly a factor of 1.7 relative to 
6dFGS, which is consistent with scaling by a/V^t and is comparable to the SDSS- 
LRG sample. We calculate the significance of the BAO detection for each log-normal 
realisation by performing fits to the correlation function using ACDM parameters 
and Qf, = 0, in exactly the same manner as the 6dFGS analysis described earlier. 
We find a 3.5 ± 0.8cx significance for the BAO detection for TAIPAN 1, 2.1 ± 0.7a 
for TAIPAN2 and 2.1 ± 0.7cr for WALLABY, where the error again describes the 
variance around the mean. 

We then fit a correlation function model to the mean values of the log-normal 
realisations for each survey, using the covariance matrix derived from these log- 
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normal realisations. We evaluated the correlation function of WALLABY, TAIPAN2 
and TAIPAN1 at the effective redshifts of 0.1, 0.12 and 0.14, respectively. With these 
in hand, we are able to derive distance constraints to respective precisions of 7%, 6% 
and 3%. The predicted value for WALLABY is not significantly better than that 
from 6dFGS. This is due to the significance of the 6dFGS BAO peak in the data, 
allowing us to place tight constraints on the distance. As an alternative figure- 
of- merit, we derive the constraints on the Hubble constant. All surveys recover 
the input parameter of Hq = 70 km s _1 Mpc _1 , with absolute uncertainties of 3.7, 
3 and 2.2 kms^Mpc" 1 for WALLABY, TAIPAN2 and TAIPAN1, respectively. 
Hence, TAIPAN1 is able to constrain the Hubble constant to 3% precision. These 
constraints might improve when combined with Planck constraints on Q^h 2 and 
Q m h 2 which will be available when these surveys come along. 

Since there is significant overlap between the survey volume of 6dFGS, TAIPAN 
and WALLABY, it might be interesting to test whether the BAO analysis of the 
local Universe can make use of a multiple tracer analysis, as suggested recently 



by Arnalte-Mur et al. 2011 . These authors claim that by employing two different 
tracers of the matter density field - one with high bias to trace the central over- 
densities, and one with low bias to trace the small density fluctuations - one can 
improve the detection and measurement of the BAO signal. Arnalte-Mur et aL 



2011 test this approach using the SDSS-LRG sample (with a very large bias) and 



the SDSS-main sample (with a low bias). Although the volume is limited by the 
amount of sample overlap, they detect the BAO peak at 4. la. Likewise, we expect 
that the contrasting high bias of 6dFGS and TAIPAN, when used in conjunction 
with the low bias of WALLABY, would furnish a combined sample that would be 
ideal for such an analysis. 

Neither TAIPAN nor WALLABY are designed as BAO surveys, with their pri- 
mary goals relating to galaxy formation and the local Universe. However, we have 
found that TAIPAN1 would be able to improve the measurement of the local Hubble 
constant by about 30% compared to 6dFGS going to only slightly higher redshift. 
WALLABY could make some interesting contributions in the form of a multiple 
tracer analysis. 
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2.9 Conclusion 

We have calculated the large-scale correlation function of the 6dF Galaxy Survey 
and detected a BAO peak with a significance of 2.4a. Although 6dFGS was never 
designed as a BAO survey, the peak is detectable because the survey contains a large 
number of very bright, highly biased galaxies, within a sufficiently large effective 
volume of 0.08/?T 3 Gpc 3 . We draw the following conclusions from our work: 

• The 6dFGS BAO detection confirms the finding by SDSS and 2dFGRS of 
a peak in the correlation function at around 105/i -1 Mpc, consistent with 
ACDM. This is important because 6dFGS is an independent sample, with 
a different target selection, redshift distribution, and bias compared to pre- 
vious studies. The 6dFGS BAO measurement is the lowest redshift BAO 
measurement ever made. 

• We do not see any excess correlation at large scales as seen in the SDSS-LRG 
sample. Our correlation function is consistent with a crossover to negative 
values at 140/i -1 Mpc, as expected from ACDM models. 

• We derive the distance to the effective redshift as Dy(z c s) = 456 ± 27 Mpc 
(5.9% precision). Alternatively, we can derive r s (z d )/Dy(z e fr) = 0.336 ±0.015 



(4.5% precision). All parameter constraints are summarised in Table 2.1 



Using a prior on Q m h 2 from WMAP-7, we find Q m = 0.296 ± 0.028. Inde- 
pendent of WMAP-7, and taking into account curvature and the dark energy 
equation of state, we derive Q m = 0.287 + 0.039(l + w) + 0.039fi fc ±0.027. This 
agrees very well with the first value, and shows the very small dependence on 
cosmology for parameter derivations from 6dFGS given its low redshift. 

We are able to measure the Hubble constant, H = 67.0 ± 3.2 km s^Mpc" 1 , 
to 4.8% precision, using only the standard ruler calibration by the CMB (in 
form of fl m h 2 and Vt^h 2 ). Compared to previous BAO measurements, 6dFGS 
is almost completely independent of cosmological parameters (e.g. f2 fc and w), 
similar to Cepheid and low- z supernovae methods. However, in contrast to 
these methods, the BAO derivation of the Hubble constant depends on very 
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basic early Universe physics and avoids possible systematic errors coming from 
the build up of a distance ladder. 

By combining the 6dFGS BAO measurement with those of WMAP-7 and 



previous redshift samples Percival et al. [from SDSS-DR7, SDDS-LRG and 



2dFGRS; 2010 , we can further improve the constraints on the dark energy 
equation of state, w, by breaking the H Q — w degeneracy in the CMB data. Do- 
ing this, we find w = — 0.97±0.13, which is an improvement of 24% compared 
to previous combinations of BAO and WMAP-7 data. 

We have made detailed predictions for two next-generation low redshift sur- 
veys, WALLABY and TAIPAN. Using our 6dFGS result, we predict that both 
surveys will detect the BAO signal, and that WALLABY may be the first 
radio galaxy survey to do so. Furthermore, we predict that TAIPAN has the 
potential to constrain the Hubble constant to a precision of 3% improving the 
6dFGS measurement by 30%. 
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Chapter 3 

z « measurements of the growth 
rate and erg 



Beutle r et al.\ 
MNRAS 423, 3430B (2012) 



Abstract 

We present a detailed analysis of redshift-space distortions in the two-point correlation 
function of the 6dF Galaxy Survey (6dFGS). The .fT-band selected sub-sample which we 
employ in this study contains 81 971 galaxies distributed over 17 000deg 2 with an effective 
redshift z e g = 0.067. By modelling the 2D galaxy correlation function, £(r p ,7r), we mea- 
sure the parameter combination f(z e s)as(z e s) = 0.423 ± 0.055, where / ~ SXn(z) is the 
growth rate of cosmic structure and a§ is the r.m.s. of matter fluctuations in 8/i _1 Mpc 
spheres. 

Alternatively, by assuming standard gravity we can break the degeneracy between as 
and the galaxy bias parameter, b. Combining our data with the Hubble constant prior 
from Riess et al. 2011] , we measure erg = 0.76 ± 0.11 and O m = 0.250 ± 0.022, consistent 



with constraints from other galaxy surveys and the Cosmic Microwave Background data 
from WMAP7. 

Combining our measurement of fas with WMAP7 allows us to test the cosmic growth 
history and the relationship between matter and gravity on cosmic scales by constraining 
the growth index of density fluctuations, 7. Using only 6dFGS and WMAP7 data we 
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find 7 = 0.547 ± 0.088, consistent with the prediction of General Relativity. We note 
that because of the low effective redshift of 6dFGS our measurement of the growth rate is 
independent of the fiducial cosmological model (Alcock-Paczynski effect). We also show 
that our conclusions are not sensitive to the model adopted for non-linear redshift-space 
distortions. 

Using a Fisher matrix analysis we report predictions for constraints on fag for the WAL- 
LABY survey and the proposed TAIPAN survey. The WALLABY survey will be able to 
measure fag with a precision of 4 — 10%, depending on the modelling of non-linear struc- 
ture formation. This is comparable to the predicted precision for the best redshift bins 
of the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that low-redshift 
surveys have a significant role to play in future tests of dark energy and modified gravity. 

3.1 Introduction 

The distribution of matter in the cosmos depends on the gravitational interaction 
and the expansion history of the Universe. Assuming that galaxies trace the mass 
distribution, a measurement of galaxy clustering can be used to derive fundamental 
properties of the Universe. 

On large scales the movement of galaxies is dominated by the Hubble recession, 
while on small scales the gravitational field introduces so-called peculiar velocities. 
Individual galaxy redshifts combine both Hubble recession and peculiar velocities 
indistinguishably. However, these effects can be statistically distinguished in a large 
sample of galaxy redshifts. This is the purpose of this paper. The difference between 
the redshift-inferred distance and the true distance is known as redshift-space dis- 
tortion. Redshift-space distortions effectively couple the density and velocity fields, 
complicating the models needed to accurately describe observed galaxy samples. On 
the other hand they permit measurements of the properties of the galaxy velocity 
field, which are difficult to access otherwise. In standard gravity we can use the 
amplitude of peculiar velocities to measure parameters that describe the matter 
content of the Universe such as Q m and a$- 

In this paper we report measurements of the parameter combination f(z e s)as(z e g), 
where / = d\n(D) / d\n(a) is the growth rate of cosmic structure (in terms of the 
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linear growth factor D and cosmic scale factor a) and as is the r.m.s. of the matter 
fluctuations in spheres of 8/i _1 Mpc. The measurement of f(z e s)as{z e s) can be used 
to test theories of dark energy and modified gravity, since a stronger gravitational 
interaction causes a larger growth rate /. It is interesting to note in this context the 
fact that a different form of gravitational interaction on large scales could be respon- 



sible for the accelerating expansion of the Universe (e.g. Dvali, Gabadadze & Porrati 



2000| |Wang|2008[ ). Probes such as type la supernovae, the Cosmic Microwave Back- 



ground (CMB) or Baryon Acoustic Oscillations, which have proven the existence of 
the current acceleration of the expansion of the Universe (e.g. Blake et al. 2011c), 
cannot distinguish between acceleration due to a dark energy component with nega- 
tive pressure or due to a modification of General Relativity. However, measurements 
of the growth of structure are able to distinguish between these models. 

In order to measure f(z e ff)a$(z e g) we have to model the effect of redshift-space 
distortions on the correlation function. While on large scales linear theory can be 
used to model these effects, on smaller scales non-linear contributions complicate 
the process. Several new approaches have been suggested in recent years to extend 
linear theory. We will discuss some of these models and apply them to our dataset. 

Redshift-space distortions have previously been analysed using both the correla- 
tion function and power spectrum using data from the 2dF Galaxy Redshift Survey 



(2dFGRS; |Peacock et al 2001 Hawkins et aL]|2003 Cole et al. 2005) and the Sloan 



Digital Sky Survey (SDSS; Tegmark et al. 2004; Zehavi et al. 2005 Tegmark et 



al. 2006 Cabre & Gaztanaga 2009 Song et al. 2010 Samushia et al. 2011). More 



recently it has become possible to do similar studies at higher redshift using the 



VVDS Guzzo et al., 2008 , WiggleZ Blake et al., 2011a and VLT VIMOS sur- 



veys Bielby et al, 2010 



Redshift-space distortion measurements are also sensitive to the overall ampli- 



tude of the clustering pattern of matter, commonly parameterised by as Lahav et 



al. , 2002 . This parameter is used to normalise the amplitude of clustering statistics 
such as the correlation function, £ oc er|. From the CMB we have a very accurate 
measurement of the matter fluctuations in the early universe (the scalar amplitude 
A s ) at the time of decoupling, z*. In order to derive erg (-2=0) we have to extrapo- 
late this measurement to redshift zero, involving assumptions about the expansion 
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history of the Universe. The CMB constraint on as heavily depends on these as- 
sumptions. Hence there is a clear advantage in obtaining low-redshift measurements 
of this parameter. The 6dF Galaxy Survey, which we analyse in this study, is one 
of the largest galaxy redshift surveys available. Its very small effective redshift and 
wide areal coverage (41% of the sky) make it a powerful sample for the study of the 
local galaxy distribution. 

Galaxy surveys usually have to consider degeneracies between redshift-space 
distortions and the Alcock-Paczynski effect, which arises from the need to assume a 
cosmological model to transform redshifts into distances. At low redshift this effect 
is very small, meaning that our measurement is fairly independent of the choice of 
the fiducial cosmological model. 

While at high redshift {z > 1) the matter density dominates both the expan- 
sion of the Universe and the growth of perturbations, at low redshift these two are 
partially decoupled, with dark energy mostly dominating the background expansion 
and the matter density dominating the growth of perturbations. As a result in 
ACDM and most proposed modified gravity models, low redshift measurements of 
the growth rate have a better constraining power than high redshift measurements. 
The measurement of the growth rate in 6dFGS therefore not only provides a new 
independent data point at very low redshift, but also promises to make a valuable 
contribution to tests of General Relativity on cosmic scales. 



The outline of this paper is as follows: In section 3.2 we introduce the 6dF 



Galaxy Survey. In section 3.3 we describe the details of the correlation function 
estimate and introduce the 2D correlation function of 6dFGS. In section 13.41 we 
derive the covariance matrix for the 2D correlation function. In section l3~5l we sum- 
marise the theory of redshift-space distortions, including extensions to the standard 
linear approach. We also discuss wide-angle effects and other systematics, such as 



the Alcock-Paczynski effect. In section |3.6| we fit the 2D correlation function to 
derive ge(z e g) = f(z e s)as{z e s) and as- Cosmological implications are investigated in 



section 3.7 In section 3.8 we make Fisher matrix predictions for two future low red- 
shift galaxy surveys, WALLABY and the proposed TAIPAN survey. We conclude 
in section I3~9l 

Throughout the paper we use r to denote real space separations and s to denote 
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separations in redshift-space. Our fiducial model assumes a flat universe with = 
0.27, w Rd = -1 and fif d = 0. The Hubble constant is set to H = 100/ikm s^Mpc -1 . 



3.2 The 6dF Galaxy Survey 



The 6dF Galaxy Survey [6dFGS; [Jones et al.[ |2004] [2006} [2009] is a near-infrared 
selected (JHK) redshift survey of 125 000 galaxies across four-fifths of the southern 
sky, with secondary samples selected in bj and rp. The |6| < 10° region around the 
Galactic Plane is avoided by the JHK surveys to minimise Galactic extinction and 
foreground source confusion in the Plane (as is |6| < 20° for 6j and rp). The near- 
infrared photometric selection was based on total magnitudes from the Two-Micron 



All-Sky Survey - Extended Source Catalog [2MASS XSC; |Jarrett et al.] |2000) . The 
spectroscopic redshifts of 6dFGS were obtained with the Six-Degree Field (6dF) 
multi-object spectrograph of the UK Schmidt Telescope (UKST) between 2001 and 
2006. The effective volume of 6dFGS is about the same as the 2dF Galaxy Redshift 



Survey [2dFGRS; Colless et al. 2001 and is a little under a third that of the Sloan 
Digital Sky Survey main spectroscopic sample at its Seventh Data Release [SDSS 
DR7; |Abazajian et al. 2009 . A subset of early- type 6dFGS galaxies (approximately 
10 000) have measured line-widths that will be used to derive Fundamental Plane 
distances and peculiar motions. 

The 6dFGS fC-selected sample used in this paper contains 81 971 galaxies se- 
lected to a faint limit of K = 12.75. The 2MASS magnitudes are on the Vega 
system. The mean completeness of the 6dFGS is 92 percent and median redshift 
is z = 0.05. Completeness corrections are derived by normalising completeness- 
apparent magnitude functions so that, when integrated over all magnitudes, they 
equal the measured total completeness on a particular patch of sky. This procedure 



is outlined in the luminosity function evaluation of Jones et al. 2006 and also in 



Jones et al., (in prep). The original survey papers Jones et al. 2004, 2009 describe 



in full detail the implementation of the survey and its associated online database. 

The clustering in a galaxy survey is estimated relative to a random (unclustered) 
distribution which follows the same angular and redshift selection function as the 
galaxy sample itself. We base our random mock catalogue generation on the 6dFGS 
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Figure 3.1: The solid black line shows the 6dFGS redshift distribution, while the 
dashed black line shows one of the random mock catalogues containing the same 
number of galaxies. The blue solid and dashed lines show the distribution after 
weighting with Pq = 1600/i 3 Mpc~ 3 (see section 
ployed weighting scheme). 



3.3.1 for more details on the em- 



luminosity function, where we use random numbers to pick volume-weighted red- 
shifts and luminosity function-weighted absolute magnitudes. We then test whether 
the redshift-magnitude combination falls within the 6dFGS K-band faint and bright 
apparent magnitude limits (8.75 < K < 12.75). 



Figure 3.1 shows the redshift distribution of the 6dFGS if-selected sample (black 
solid line) compared to a mock catalogue with the same number of galaxies (black 
dashed line). 



3.3 Correlation function measurement 



We calculate the co-moving distances for each galaxy using the measured redshift 



c f z dz' 

D ^wJ a W) (3 ' i; 
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with 



E(z) 



[n 



™(i + ^) 3 



fid! 1/2 



(3.2) 

1 — and describe dark 

energy as a cosmo logical constant (w = — 1). Given the low redshift of our dataset, 
these assumptions have a very small impact on our final results (see section 3.5.5). 



where we assume a flat universe with Q^ d = and Qjf 

fid 



We define the positions of two galaxies as si and s~2. The redshift-space separa- 
tion is then given by h = si — s*2, while s = (s*i + s*2)/2 is the mean distance to the 
pair of galaxies. Now we can calculate the separation along the line-of-sight 7r and 
the separation across the line-of-sight r p 



7T 



\s- h\ 



h\ 



IT' 



(3.3) 
(3.4) 



The absolute separation is then given by s = */7r 2 + r%. 

We measure the separation between all galaxy pairs in our survey and count 
the number of such pairs in each separation bin. We do this for the 6dFGS data 
catalogue, a random catalogue with the same selection function, and a combination 
of data-random pairs. We call the pair-separation distributions obtained from this 
analysis step DD, RR and DR, respectively. In the analysis we used 30 random 
catalogues with the same size as the real data catalogue and average DR and RR. 



The redshift-space correlation function itself is then given by Landy & Szalay 1993 



P' =14- 

Sdata T RR \ nd 



DD fn r \ 2 DR fur 
~ ' ~ 2 RR Vr7 d 



(3.5) 



where the ratio nr/n^ is given by 



n, 



E 



<1 W l 



(3.6) 



and the sums go over all random (N r ) and data (N d ) galaxies. The galaxies are 
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weighted by the inverse completeness Cj of their area of the sky 



Wi(z) = Ci 



(3.7) 



We will discuss further weighting techniques in the next section. 

There is a possible bias in the estimation of the correlation function due to the 
fact that we estimate both the mean density and the pair counts from the same 
survey. This leads to a non-zero difference between the true correlation function 
estimate of an ensemble of surveys and the ensemble average of £(s) from each 
survey. This is commonly known as the integral constraint (e.g. |Peebles 1980), 



which can be calculated as (see e.g. Roche et al. 2002) 



ic 



£6 



model 



RR 



ERR 

and enters our correlation function estimate as 



(3.8) 



6 



data 



£data + * C > 



(3.9) 



where £^ ata is the redshift-space correlation function from eq. 3.5 and £ mo dei is the 



model for the correlation function. In 6dFGS ic is typically around 6 x 10 4 and so 
has no significant impact on the final result. 



In Figure 3.2 we show the 2D correlation function calculated from the 6dFGS 
dataset. In this Figure we use bins of 0.5/i _1 Mpc, while for the analysis later on 



we use larger bins of 2h 1 Mpc (see Figure 3.6). The figure shows clearly the two 



effects of redshift-space distortions which we will discuss later in section 3.5, the 
"finger-of-God" effect at small r p , and the linear infall effect at larger r p which gives 
the correlation function a non-circular shape. 
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r p [h^Mpc] 



Figure 3.2: The 2D correlation function of 6dFGS using a density weighting with 
Po = 1600/i 3 Mpc~ 3 . For reasons of presentation we binned the correlation function 
in 0.5/i _1 Mpc bins, while in the analysis we use larger bins of 2h~ l Mpc. Both 
redshift-space distortion effects are visible: the "finger-of-God" effect at small angu- 
lar separation r p , and the anisotropic (non-circular) shape of the correlation function 
at large angular separations. 



3.3.1 Density weighting 

In Fourier space the error in measuring the amplitude of a mode of the linear power 
spectrum^] is given by 

a P{k) ^{b + f^fp{k) + {N), (3-10) 

where b is the linear bias, / is the growth rate, /i is the cosine of the angle to the 
line of sight and P(k) is the matter power spectrum. The first term on the right 
hand side of this equation represents the sample- variance error, while the second 

■"■As the correlation function and power spectrum are related by a Fourier transform, the fol- 
lowing discussion also holds true for a correlation function measurement. 
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term ((N)) represents the Poisson error. 

If the sample- variance error is dominant we can reduce the power spectrum error 
by employing a weighting scheme which depends upon the galaxy density n(z), such 



as the one suggested by Feldman, Kaiser & Peacock 1994 



where Pq describes the amplitude of the weighting. A stronger weighting (larger 
value of Pq) yields a smaller sample- variance error since it increases the survey 
volume by up-weighting sparsely sampled regions. However, such a weighting scheme 
also increases the Poisson error because it shifts the effective redshift to larger values 



with a smaller galaxy number density. This is illustrated in Figure 3.3(a) and 3.3(b) 
Such a weighting scheme is standard for large scale structure analyses. 

In a magnitude-limited sample such as 6dFGS, up- weighting higher redshift 
galaxies also has the effect of shifting the galaxy bias to larger values. The sample- 
variance error is proportional to the clustering amplitude, and so a larger bias results 
in a larger error. However, the weighting will still ensure that the relative error of 
the power spectrum, ap^)/P{k), is minimised. The redshift-space distortion signal 
is inversely proportional to the galaxy bias, ~ Q'J n (z)/b. If weighting increases the 
bias b, it also reduces the signal we are trying to measure. We therefore must inves- 
tigate whether the advantage of the weighting (the reduced relative error) outweighs 
the disadvantage (increasing galaxy bias). 

The situation is very different for measuring a signal that is proportional to the 
clustering amplitude, such as the baryon acoustic peak. In this case the error and the 
signal are proportional to the bias, and so weighting will always be beneficial. We 
stress that an increasing bias with redshift is expected in almost all galaxy redshift 
surveys. Therefore redshift-space distortion studies should first test whether galaxy 
weighting improves the measurement. The 6dF Galaxy Survey is quite sensitive to 
the weighting scheme employed because it has a high galaxy density, making the 
sample- variance error by far the dominant source of error. 

Finally, we have to consider the correlation between the bins in the measured 
power spectrum or correlation function. If the error is sample-variance dominated, 
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the bins will show large correlation (especially in the correlation function), while in 
the case of Poisson-noise dominated errors, the correlation is much smaller. Weight- 
ing will always increase the Poisson noise and hence reduce the correlation between 
bins. 

All of the density weighting effects discussed above need to be considered, before 
deciding which weighting is best suited to the specific analysis. We summarise as 
follows: 

• If a galaxy sample is sample-variance limited, density weighting will reduce 
the (relative) error of the clustering measurement. 

• In most galaxy redshift surveys, the density weighting increases the galaxy 
bias, which reduces the redshift-space distortion signal, by flattening the clus- 
tering anisotropy. 

• Galaxy weighting also reduces the correlation between bins in the power spec- 
trum and correlation function. 



Whether the signal-to-noise is actually improved by density weighting depends on 
the specific sample. For the 6dFGS redshift-space distortion analysis we found 
Po ~ 1600/i -3 Mpc 3 leads to the most accurate constraint on the growth rate. 

This discussion also indicates that low-biased galaxy samples have an advantage 
over a highly biased galaxy sample in measuring redshift-space distortions. We will 



discuss this point further in section 3.8 



3.4 Error estimate 

In this section we will derive a covariance matrix for the 2D correlation function 
using jack-knife re-sampling. We also use log-normal realisations to test the jack- 
knife covariance matrix. 
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Figure 3.3: (a) The relative error in the 2D correlation function as a function of 
line-of-sight separation 7r at a fixed r p = ll/i _1 Mpc. Other regions of the 2D 
correlation function behave in a similar manner. The solid lines show the Poisson 
error for different values of Po, while the data points show the total (Poisson + 
sample variance) error obtained as the diagonal of the covariance matrix derived 
using jack-knife re-sampling. The purpose of the weighting (P ) is to minimise the 
total error, which is achieved for a value of Po ~ 1600/i~ 3 Mpc 3 . The weighting 
reduces the error by almost a factor of two on most scales. The dashed line shows 
the error derived from log-normal realisations using Po — 1600/i~ 3 Mpc 3 and is in 
very good agreement with the jack-knife error, (b) same as (a) for a fixed line-of- 
sight separation 7r = 1 1/z 1 Mpc. 
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3.4.1 Jack-knife re-sampling 

We divide the dataset into N = 480 subsets, selected in R.A. and Dec. Each re- 
sampling step excludes one subset before calculating the correlation function. The 
N — 1 remaining subsets have a volume which is (N — 1)/N times the volume of the 
original data. The covariance matrix is then given by 

Cij = £ [**(**) - ec«)] - £(%)] , (3.12) 
fc=i 

where £ fc (<Si) is the correlation function estimate at separation s< with the exclusion 
of subset fc. The mean value is defined as 

1 - 

k=l 

The case z = j gives the error ignoring correlations between bins of = 



3.4.2 Log-normal realisations 



We can create a log-normal realisation Coles & Jones, 1991 Cole et al. 2005 



Percival et al. , 2007; Blake et al. 2011a Beutler et al. 2011] of a galaxy survey 
by deriving a density field from a model power spectrum, P(k), assuming Gaussian 
fluctuations. This density field is then Poisson sampled, taking into account the 
window function and the total number of galaxies. The assumption that the input 
power spectrum has Gaussian fluctuations can only be used if the fluctuations are 
much smaller than the mean density, otherwise the Gaussian model assigns a non- 
zero probability to regions of negative density. A log- normal random field, LN(x), 
can avoid this unphysical behaviour. It is obtained from a Gaussian field G(x) by 



LN(x) = exp[G(f) - a 2 G /2] - 1, 



(3.14) 



where <jq is the variance of the field. LN(x) is positive-definite but approaches 
G(x) whenever the perturbations are small. As an input power spectrum for the 
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log-normal field we use the linear model Kaiser, 1987 , given by 



(3.15) 



with b = 1.47 and = 0.35, which is consistent with the parameters we measure 



for 6dFGS (see section 3.6). Pss(k) is a linear density power spectrum in real space 



obtained from CAMB Lewis, Challinor & Lasenby, 2000 and P g (k,fi) is the galaxy 



power spectrum in redshift-space. Our method is explained in more detail in Beutler 



et al. 2011 , appendix A. 



We produce N = 1500 such realisations and calculate the 2D correlation func- 
tion for each of them, deriving a covariance matrix 



1 



N 



n=l 



(3.16) 



where the mean value £ l (r p , it) is defined as 

1 N 



(3.17) 



n=l 



and £n(r p , 7r) is the 2D correlation function estimate of realisation n at a specific 
separation (r p , ir). 



3.4.3 Discussion: Error analysis 



Figure 3.4 shows the correlation matrix 



= ^S=^> ( 3 - 18 ) 

derived from the two different covariance matrices, Cij. We plot the jack-knife result 
in the upper-left corner and the log-normal result in the lower-right corner. Both the 
correlation and covariance matrix are symmetric. While the log-normal correlation 
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Figure 3.4: The correlation matrix for the 2D correlation function 7r) with a 
bin size of 2 x 2h~ 1 Mpc. The upper-left corner shows the jack-knife estimate, while 
the lower-right corner shows the result of using 1500 log-normal realisations. Since 
this plot shows the correlation of all 15 x 15 bins it contains 225 x 225 entries. 



matrix indicates somewhat more correlation between bins, overall the correlation 
matrices are in rough agreement. The diagonal errors for both the log- normal and 



jack-knife covariance matrices are plotted in Figures 3.3(a) and 3.3(b) and show 
very good agreement at large scales. 



Every row/column in Figure 3.4 shows the correlation of one bin with all other 
bins in £(r p ,7r). Figure 3.5 shows an example of such a row/column obtained from 
jack-knife re-sampling as a 15 x 15 matrix, in this case for bin 127 (r p = 13/iT 1 Mpc, 
7T = m^Mpc). 

Log-normal realisations do not account for non-linear mode coupling and are 
very model-dependent in the quasi-linear and non-linear regime. Since our analysis 
relies on fits to fairly small scales, we decided to use the jack-knife covariance matrix 
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Figure 3.5: This plot shows the correlation of bin 127 



13/i _1 MpC, 7T 



17h 1 Mpc) with all other bins in the 2D correlation function, derived using jack- 
knife re-sampling. It corresponds to row/column 127 of the jack-knife correlation 



matrix which is shown in Figure 3.4| (upper- left corner). 



in our analysis. However, we find that none of the results reported in this paper 
depend significantly on which of the two covariance matrices is used. 



3.5 Modelling the 2D correlation function 

In this section we discuss the theory of redshift-space distortions, starting with 
the standard linear perturbation theory. We then describe different approaches 
for extending the linear model to include non-linear structure formation. We also 
discuss deviations from the plane-parallel approximation. 
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3.5.1 Linear redshift-space distortions 



The position of a galaxy in real space r 
redshift-space, s, via 

s = r + 



(x, y, z) is mapped to the position in 



Vz{r) 
aH(z) 



(3.19) 



where the unit vector z indicates the line-of-sight direction, and the quantity v z is the 
line-of-sight component of the velocity field, namely v z = v ■ z. The scale factor a is 
defined as 1/(1+2;) and H(z) is the Hubble constant at redshift z. The second term in 
the equation above represents peculiar velocities caused by gravitational interaction. 
On small scales this elongates structures along the line-of-sight and leads to the so- 



called "finger-of-God" effect Jackson, 1972 . On large scales matter falls in towards 



over-dense regions and systematically influences our distance measurement, making 
the over-densities appear more over-dense. The latter effect can be described by 
linear theory, while the "finger-of-God" effect is a non-linear phenomenon. 



The model of linear redshift-space distortions has been developed by Kaiser 



1987 and Hamilton 1992 assuming a plane-parallel approximation. In Fourier 



space the linear model can be written as 



(3.20) 



where Pss(k) is the matter density power spectrum and P g (k, /x) is the galaxy density 
power spectrum. In this model, linear redshift-space distortions are quantified by 
the parameter /3, which is defined as 



ldlnD(z) _ W m {z) 
b dln(a) b 



(3.21) 



where b is the linear galaxy bias factor, D(z) is the growth factor and 7 is the 



gravitational growth index, which takes the value 7 = 0.55 in ACDM Linder, 2005 
f2 m (z) is the matter density at redshift z, and is defined as 



H 



^fi m (, = 0)(l + z) 3 



(3.22) 
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with 

Jk-=[n m (i + z ) s + n A ]- 1 , (3.23) 

where we again follow our fiducial model of f2| d = and w m — — 1. Modifications 
of the gravitational force mainly affect 7, while changes in the expansion history of 
the Universe affect Q m (z). 



3.5.2 Parameterisation 

From now on we will formulate our equations in terms of ge(z) = f(z)a 8 (z) and 
gb(z) = ba$(z). We choose the parameter set \g$, g^\ instead of [(3, cr 8 , b], because as 



and the linear bias, b, are degenerate and difficult to disentangle White, Song & Per 



rival, 2008 Song & Percival, 2009) Song et al. , 2010 . Within this parameterisation 



the power spectrum is expressed as 



P g (k) = b 2 P 5S (k)=g 2 b Q S5 (k), 



(3.24) 



where Qss(k) is the unnormalised matter density power spectrum (see eq. 3.30). In 
terms of our new parameters, we can write 



9b 



(3.25) 



3.5.3 Extensions to linear theory 



For the remainder of this section we will discuss possible extensions of the linear 
model to include non-linear structure formation and wide-angle effects. The simplest 



model of non-linearities is the so-called streaming model Peebles, 1980 Hatton & 



Cole, 1998 



Here, the redshift-space correlation function is just a convolution of 
the linear correlation function in redshift-space with a pairwise velocity probability 
density function, F(v), 



H(z eS )a eS J 



F(v)dv, 



(3.26) 
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where a e g = 1/(1 + z e ff) is the scale factor and H(z c g) is the Hubble constant at 



the effective redshift. We chose F(v) to be an exponential distribution Peacock & 

V2\v 



DoddsJ|1996] 



F(v) 



a 



,y/2 



cxp 



On 



(3.27) 



which has been shown to successfully describe observations Davis & Peebles, 1982 



Fisher et al. , 1994; Marzke et al. , 1995; Landy, 2002 . Although the pairwise velocity 



dispersion within a halo is expected to follow a Gaussian distribution instead of an 
exponential, galaxies populate halos of a wide range of masses and velocity disper- 



sions, which combine to approximately form an exponential function Sheth, 1996 



Diaferio & Geller 


1996 


Seto & Yokoyama, 


1998 . 


galaxy type (e.g. 


Madgwick et al. 2003 


) and hence its 



The parameter a p depends on 



is limited. 



In recent years many improvements to the model discussed above have been 
suggested. We will initially discuss these models in Fourier space because the theo- 
retical motivation is clearer and many expressions needed for these models simplify 
considerably. However, since our data is in real space we will also give real-space 



expressions later on. We start with the streaming model of eq. 3.26 which in Fourier 
space is given by 



f,(M) = f(i + ^W 1 + t ,^ . 



(3.28) 



Hence in Fourier space, the convolution with an exponential function becomes a 
multiplication by a Lorentzian distribution. 

The model above assumes that there is a perfect correlation between the velocity 
field and the density field, which is given by Pss(k) = Pse(k) = Pee(k) where Pss(k) 
is the matter density power spectrum as before. Pee(k) = (l^l 2 ) is the velocity 
divergence power spectrum where 9 = V • v is the velocity divergence, and Pse(k) is 
the cross power spectrum. Non-linear effects will violate these assumptions, since the 
density power spectrum is expected to increase in amplitude at small scales because 
of non-linear effects, while the velocity field becomes randomised at small scales 
(e.g. within virialized galaxy clusters) and hence Pee{k) will decrease in amplitude 
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e.g. Carlson|2009 ) . Scoccimarro 2004 suggested expressing the 2D power spectrum 



without the assumption of linear relations between the density field and velocity 
field, by 



P g (k,fi) = F q {k,iJL,a v ) 

x [b 2 P 55 (k) + 2 l x 2 bfP 5e (k) + fx 4 f 2 P ee (k)] 
= F q (k,ii,a v ) 

x [dQss(k) + 2[i 2 g b g e Q S o(k) + ^g 2 e Q ee {k)\ 



(3.29) 



where the different Q xy are defined as 

Qss(k) = P S s(k)/a 8 (z cS ) 2 , 

Qse(k) = P se (k)/a s (z cS ) 2 , (3.30) 
Qee(k) = P ee (k)/cr 8 (z cS ) 2 . 

and the damping function F q (k, fi, a v ) is usually chosen to be a Gaussian of the form 

F g (k,fi,a v ) = e-^ 2 . (3.31) 

The parameter a v quantifies the non-linear dispersion in the bulk motion of halos. 
It is different to the a p parameter we introduced earlier that describes small-scale 
randomised motion (e.g. of galaxies within a halo). We can derive a v from the 
velocity power spectrum as 



6tt 2 



<r v (z) = -9^r / Qee(k)dk. (3.32) 



o 



Jennings, Baugh & Pascoli 2011a provide fitting formulae for Pse(k) and Pee(k) 
derived from N-body simulations. They find the following relation between the 
different power spectra 



a 2 + a 3 P S s{k) 



P,y{k) = " \T n"? ■ ( 3 - 33 ) 



where Pss(k) can be obtained from CAMB by including halofit Smith et al. , 2003 . 
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For the cross power spectrum P xy {k) = Pse{k) we use the (updated) parameters (a>o, 
a Xl a 2 , a 3 ) = (-12 483.8, 2.55430, 1381.29, 2.54069) and for P xy {k) = P ee (k) we 
use (-12480.5, 1.52404, 2165.87, 1.79640). The fitting formula reproduces P ge to 
better than 1% for k < 0.4/iMpc^ 1 , to 10% for 0.4 < k < 0.7/iMpc -1 , and to 15% 
for 0.7 < k < 1/iMpc -1 . It also reproduces P$o to less than 4% over the whole range 
k < 1/iMpc -1 (Jennings, private communication). We cut off the integral when the 
Jennings formula predicts negative values, although because of the high precision of 
the fitting formula up to large k, such a cut-off will not affect our measurement. 



We can express eq. 3.29 in real-space as 



£sc(r P ,7r) 



+ 



-,9b9e^2,5 



4 
7 l 



(r) 



5* 



(3.34) 



where Ve(fi) are the Legendre polynomials and the spherical harmonic moments 
£ i>xy {r) are given by 



xy\ 



OG 



1 k 2 dkdjjL 



-(fc/i(T„) 2 



/o J-i (2vr) 2 
x cos(krii)Q xy (k)Ve(n). 



(3.35) 



Appendix [D] shows a partial analytic solution for the double integral above. 

We therefore have two different models which we will apply to our data: The 
simple streaming model £, s t( r P , 7! ') and the Scoccimarro models £sc( r p, 7 0- Note that 
Csc(^p,tt) does not include the parameter a p and hence has one less free parameter 
than the streaming model. 

All equations above are based on the plane-parallel approximation. In the case 
of 6dFGS we also need to account for wide-angle effects. This means we will replace 
the equations above with more general descriptions, which we discuss in the next 
section. 
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3.5.4 Wide-angle formalism 



So far we have assumed that the separation between galaxy pairs is much smaller 
than the distance of the galaxies from the observer. The 6dF Galaxy Survey 
has a maximum opening angle of 180° and the (effective) redshift is fairly low at 
z cS — 0.067 (see section 3.3). We therefore include wide-angle correction terms. 



The wide-angle description of redshift-space distortions has been laid out in sev- 



eral papers Szalay et al. 1997 Szapudi 2004 Matsubara 2004 Papai & Szapudi 



2008 Raccanelli et al. , 2010 . Here we will expand on this work by formulating the 



equations in terms of gg and and by distinguishing between the density-density, 
velocity- velocity and density- velocity contributions £ee and £$g). The model will 



then correspond to the Scoccimarro model (eq. 3.34) in the last section. 

The general redshift-space correlation function (ignoring the plane-parallel ap- 
proximation) depends on 0, 9 and s. Here, s is the separation between the galaxy 
pair, 9 is the half opening angle, and is the angle of s to the line-of-sight (see 
Figure 1 in Raccanelli et al. 2010). The angles and 9 are not independent, but 



the relation between them is usually expressed through the two angles 0i and 2 
given by = ^(0i + 2 ) and 9 = \{4>\ — 2 ). The total correlation function model, 



including 0(9 ) correction terms, is then given by Papai & Szapudi, 2008 



£(0, 9, s) = a 00 + 2a 2 cos(20) + a 2 2 cos(20) + b 2 2 sin 2 (20) 
+ - 4a 02 cos(20) - 4a 22 - 46 22 - Aa 10 cot 2 (0) 
+ 4a n cot 2 (0) - 4a 12 cot 2 (0) cos(20) + Ab u 
-8&i 2 cos 2 (0)l# 2 + O(# 4 ). 



(3.36) 



This equation reduces to the plane-parallel approximation if 9 = 0. The factors a 
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and b xy in this equation are given by 



«02 
«22 
&22 

Qio 
a n 
an 
bn 



9b9e c2 

o ^2,S 



*9? 



21 



8 f:2 
S2,( 



15 
140 



(r) + 7^4, 



" — ^~ ^2,<5e( r ) + "^"£2,00 ( r ) + ^|^4,ee( r ) 



15^ e 
9e,2 



15 



o.e 



21 ^ 
21^' 



r + 



140 

35 



£ 2 

S4,< 



e ^2 

S4,< 



3 r 2 [^o,6»e( r ) 2^2,e 
|[2^(r)-3& 

5r 



in 

,,2 



ye t i 

c33 



5r 



where the spherical harmonic moments, 67a™ ( r )> are 



r 



xy v 



(k^ia v ) 2 



>o J-i (2^) 2 
x cos(krn)Q xy (k)Ve(n), 



(3.37) 



(3.38) 



with Qss(k), Qse(k) and Qee(k) as defined in equation 3.30| (see appendix D for an 
analytic solution). 

In order to obtain a model for the 2D correlation function (including wide- 



angle effects), we can simply integrate eq. 3.36 over 6. In case of the plane-parallel 



approximation this would correspond to eq. |3.34| To reduce the equation to the 
simple streaming model we have to set a v = 0, Qss(k) = Qse(k) = Qee{k) and 
convolve with F(v). 



Samushia et al. 2011 studied wide angle effects in the SDSS-LRG sample and 
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found that such effects are very small. The lower redshift and larger sky coverage 
of the 6dF Galaxy Survey mean that wide-angle effects are certainly larger in our 
data set than SDSS-LRG. However, we found no significant impact of such effects 
on our results. This is mainly because our analysis includes only rather small scales 
(n < 30/i -1 Mpc, r p < 30/iT 1 Mpc, see section 3.6). This agrees with the findings 



of Beutler et al. 2011 



We stress that the correction terms discussed above only capture first-order 
effects and so we additionally restrict our analysis to 9 < 50°. 

Many authors have found, using N-body simulations, that both the simple 
streaming and Scoccimarro models do not capture all non-linear and quasi-linear 



effects present in redshift surveys (e.g. Jennings, Baugh & Pascoli 2011b Kwan, 



Lewis & Linder 2011 de la Torre & Guzzo 2012 and references therein). For exam- 



ple the linear bias model seems to be too simplistic at small scales. In this study 
we restrict ourselves to the models discussed above but would like to point out that 
many other models have been suggested Matsubara, 2008[b Taruya, Nishimichi & 



Saito 2010 Reid & White 2011 Seljak & McDonald 2011 which could be com- 



pared to our data in future work. 



3.5.5 Systematics and the Alcock-Paczynski effect 

A dataset such as 6dFGS, containing galaxies with a high linear bias factor, may 
be prone to scale-dependent galaxy bias on small scales. We used the GiggleZ 
simulations (Poole et al., in preparation) to derive the form of the scale-dependent 
bias. We find a 1.7% correction to a linear bias at scales of s = lOh^ 1 Mpc, and note 
that this correction has a small but non-negligible impact on our results, depending 
on the smallest scales included in the fit. 

A further systematic uncertainty comes from the so-called Alcock-Paczynski 



(AP) effect Alcock & Paczynski, 1979 . By assuming a cosmological model for the 



conversion of redshifts and angles to distances (see section 3.3) we could introduce an 
additional anisotropic signal in the correlation function which depends on the angle 
to the line-of-sight, if the true cosmology differs from our fiducial cosmological model. 
These distortions are degenerate with the linear redshift-space distortion signal and 
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hence both effects need to be modelled to avoid systematic bias Ballinger, Peacock 



& Heavens, 1996 Simpson & Peacock, 2010 . Using the two scaling factors 



h 



H hd (z) 
~Hjz) 
D A (z) 
D 



Mi 



(3.39) 
(3.40) 



we can account for the AP effect by rescaling the corresponding axis 



7T 



7T 
J\ 

f± 



(3.41) 
(3.42) 



While this is a matter of concern for high-redshift surveys, the 6dF Galaxy Survey is 
very nearly independent of the Alcock-Paczynski distortions because the distances 
7r and r p are almost independent of the fiducial cosmological model when expressed 
in units of h~ 1 Mpc. For example the value of f± at z = 0.067 changes by only 
(0.3%, 1.5%) for 10% changes in (O m , w). The corresponding values for f\\ are 
(0.7%, 1.6%). The relative tangential/ radial distortion depends on the value of 
(1 + z)Da(z)H(z)/c relative to the fiducial model: 



Da{z)H{z) _ h 
Df(z)H^(z) /|| 



(3.43) 



This parameter combination changes by only 0.3 (0.1)% for 10% changes in O m (w). 
Finally, it has been shown that 7 has a degeneracy with the equation of state 



parameter for dark energy, w Simpson & Peacock, 2010 . High redshift measure- 



ments usually assume w — — 1, which could introduce a bias in the measured value 
of 7 if dark energy is not exactly a cosmological constant. Again 6dFGS is robust 
against such effects. 
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Table 3.1: Cosmological parameters derived from the 6dFGS 2D correlation func- 
tion. The effective redshift is z e g = 0.067. The last column indicates the pri- 
ors/assumptions which go into each individual parameter measurement. The prior 
on the Hubble constant comes from Riess et al. 2011 and the WMAP7 prior 



from Komatsu et al. 2011 . The asterisks denote parameters which are derived 



from fitting parameters. 



Summary of parameter 


measurements from 6dFGS 




0.423 ± 0.055 




9b(z e s) 


1.134 ±0.073 




(3* 


0.373 ± 0.054 






0.76 ±0.11 


[H = 73.8 ±2.4, 7 = 0.55] 




0.250 ±0.022 


[H = 73.8 ±2.4, 7 = 0.55] 


b 


1.48 ±0.27 


[#o = 73.8 ±2.4, 7 = 0.55] 




0.58 ±0.11 


[go + erg from 6dFGS] 


7 


0.547 ±0.088 


[WMAP7] 




0.271 ±0.027 


[WMAP7] 



3.6 Fitting the 2D correlation function 

We now fit the two models for the 2D correlation function we developed earlier 
(£st(?>, and ^sc(^p,tt)) to our data. For the final result we use the (,sc(r p ,7i) 
model since it gives similar measurements to the streaming model with one less free 
parameter. 



Other studies (e.g. Samushia et al. 2011) prefer to analyse the correlation func 



tion moments £o>£2 (and if possible £4), which carry the same information as the 
2D correlation function. The correlation function moments have the advantage that 
the number of bins grows linearly with the highest scales analysed, while for £(r p , n) 
the number of bins grows quadratically. This makes it easier to get reliable covari- 



ance matrices. Samushia et al. 2011 also show that the measurement errors of £^ 
are more Gaussian. However, the correlation function moments are integrals over 
fi and hence carry information from all directions, including /x = 0. Finger-of-God 
distortions can influence the correlation function moments up to large scales (20 - 
30/i -1 Mpc), while in the 2D correlation function they can be excluded via a cut in 
r p . In 6dFGS we have found that these non-linear effects have a strong impact on 
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Figure 3.6: The 2D correlation function in 2h~ 1 Mpc bins. The fitting area is 
indicated by the dashed lines, where black corresponds to the streaming model, 
(£ st (r p ,7r)) and red corresponds to the Scoccimarro model (£sc( r P , 7r)). The black 
and red contours show the best fitting models for £, 



sty pi 



(black) and £s c (?>,7r) 
(red). The deviations seen in the two contours at large scales are well within the 
error bars of the two models, which can be seen in Figure |3.8[ At small scales 



(< lAh~ l Mpc) the Scoccimarro model predicts much more clustering, while in the 
real data this clustering is smeared out along the line of sight because of the finger- 
of-God effect. 



the correlation function moments up to 30h 1 Mpc. We have therefore decided to 
focus on the 2D correlation function instead of the correlation function moments. 



3.6.1 Derivation of the growth rate, go = /cr 8 



In Figure 3.6 we show the 6dFGS 2D correlation function. For our analysis we 
bin the data in 2 x 2h~ 1 Mpc bins from to 30/i^ 1 Mpc in r p and n. Including 
larger scales does not add further information. At small r p , the finger-of-God effect 
becomes dominant and we expect any linear model to fail. Since our description 
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Figure 3.7: The value of ge{z e s) — f( z es) cr 8(z e s) as a function of the cut-off scale 
obtained by fitting the 6dFGS 2D correlation function with two different mod- 



„cut 



els (as described in section 3.5.3 and 3.5.4). At large scales the two models converge 



to similar values, while on small scales the models deviate from each other because 
of the different descriptions of non-linear evolution. For the final parameter mea- 

we chose model 2, £sc(r p , ^ 



surements in Table 3.1 
I6/1 



of iv«cut 
p 



x Mpc 
the quality of the fit 



with a conservative cut-off scale 
In the lower panel we plot the reduced x 2 as an indicator of 



of non-linearities, in both of our models, is limited in its capability to capture all 
non-linear effects, it is necessary to include a cut-off scale r^ ut marking a lower limit 
of the fitting range in r p . 



Figure 3.7 shows the measured value of gg as a function of the cut-off scale 
r cut f Qr Qur ^ WQ differed models. Above Vp Ut w Mpc the streaming model, 
6rt(^pj7r), approaches a constant value of gg. Our second model, ^sctjp,^), contains 
a systematic error up to much larger scales, before it comes into agreement with the 
streaming model at about r£ ut = 16/i -1 Mpc. This is expected since this model does 
not include a description of effects in the non-linear regime. For the final constraints 
we choose r™* = 10/z 1 Mpc for the streaming model and r^ ut = 16/i -1 Mpc for the 
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Scoccimarro model. We also note that since the Scoccimarro model is based on 
only two free parameters (gg and the error is generally smaller compared to the 
streaming model, which has three free parameters (gg, gb and a p ). Other studies fit 
for the parameter a v (e.g. de la Torre &: Guzzo||2012~ ) in the Scoccimarro model, but 



we derive it using eq. 3.32 



For £sc(?>,7r) we use the fitting range < n < 30/i _1 Mpc and 16 < r p < 
30/i -1 Mpc, which results in a total of 105 bins. The best-fitting results are gg = 
0.423 ±0.055 and gt> = 1.134 ±0.073, where the errors for each parameter are derived 
by marginalising over all other parameters. The x 2 of this fit is 115 with 103 degrees 
of freedom (d.o.f.), indicating a good fit to the data. 

For £,st{ r p, 7t '), we have the fitting range < n < 30/i _1 Mpc and 10 < r p < 
30/i -1 Mpc, which results in a total of 150 bins. The best fitting parameters are 
g = 0.389 ± 0.067, g b = 1.084 ± 0.036 and a v = 198 ± 81km/s. The reduced x 2 of 
this fit is given by x 2 /d.o.f. = 202/147 = 1.37. We compare the constraints on gg 
and gb from both models in Figure |3.8 



In the Scoccimarro model we could use the parameter a v oc gg instead of gg 



to test cosmology, as suggested by Song et al. 2010 . Our best fit gave a v 



2.59 ± 0.34/i -1 Mpc. However, this parameter depends on an additional integral 
over the velocity power spectrum, which adds a theoretical uncertainty. We therefore 
prefer to use gg in the following discussions. 

We can also express our results in terms of (5 which is given by /3 = gg/gb = 
0.373 ± 0.054. We summarise all measured and derived parameters in Table [3~T| 



3.6.2 Derivation of and Q TO 

In this section we use redshift-space distortions to directly measure as- The angular 
dependence of the redshift-space distortion signal in the 2D correlation function 
allows us to measure (3, which quantifies the amplitude of redshift-space distortions. 
Together with Q m (z) and 7 = 0.55, this constrains the linear bias b through the 
equation 




(3.44) 
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Figure 3.8: Likelihood distribution of gg and derived from the fit to the 2D 
correlation function. The solid black contours show model ^(V^ 71 "^ while the 



dashed contours show the streaming model (see section 3.5.3 and 3.5.4 for details of 
the modelling). The fitting range is < 7r < 30/i _1 Mpc and 10 < r v < 30/iT 1 Mpc 



for £, 



st \i pi 



and < 7T < 30/i- 1 Mpc and 16 < r p < 30/i" 1 Mpc for £ Sc (r p , vr). The 



black cross indicates the best-fitting value for the solid black contours. 



Knowing b we can use the absolute amplitude of the correlation function, [&cx 8 (z)] 2 , 
to constrain <r 8 (z=0) = [D (z=0) / D (z e g)] x cr 8 (,2; e fr). 

For computational reasons we use our first model, £ s t(r p , it), in this sub-section 
and fit the five parameters a 8 , fi m , b, Hq and a p using an MCMC approach. Since the 
shape of the correlation function is only sensitive to T = Q m h, we cannot constrain 
Q m and H at the same time. For the final results we include a prior on the Hubble 
constant (H = 73.8 ± 2.4 km s" 1 Mpc" 1 , Riess et al. 2011 , from now on referred to 
as HST prior) and marginalise over it. We use the same binning and fitting ranges 
as in the previous section. 

The best-fitting model results in x 2 /d.o.f = 1.35. We find a% = 0.76 ± 0.11, 
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Figure 3.9: This plot shows the likelihood distribution of the galaxy bias b and as, 
which we obtained by fitting the 6dFGS 2D correlation function assuming 7 = 0.55. 
The solid black line shows the result using a prior on the Hubble constant of H Q = 
x Mpc 



from Riess et al. 2011 



73.8 ± 2.4 km s" 

a prior of H = 67.0 ± S^kms^Mpc" 1 from iBeutler et all 12011 



while the dashed black line uses 
Although the 



detection of redshift-space distortions can partially break the degeneracy between b 
and o"s which exists in the ID correlation function, there is still a significant residual 
degeneracy. The black cross marks the maximum likelihood value for the solid black 
lines. 



r2. 



0.250 ± 0.022, b = 1.48 ± 0.27 and o v = 174 ± 73km/s. The remaining 



degeneracy between the bias b and <t 8 is illustrated in Figure 3J} We include all 
these results in Table 13.11 



Figure |3.10| compares the 6dFGS Q m — a 8 probability distribution to mea- 
surements from several other datasets: The CFHT wide synoptic Legacy Survey 
(CFHTLS) |Fu et al.[ [2008] , the SFI++ peculiar velocity survey iNusser & Davis 



2011], cluster abundance from X-ray surveys Mantz et al. 2010 and WMAP7 Ko- 



matsu et al. , 2011 . Many of the experiments shown in this figure have systematic 
modelling uncertainties when extrapolating to z = arising from assumptions about 
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the expansion history of the Universe. Only 6dFGS and SFI++ are at sufficiently 
low redshift to be independent of such effects. To illustrate the impact of these 
effects we plot in Figure [3 the probability distribution Q m — as from WMAP7 for 
different cosmological models. The CMB measures the scalar amplitude A s , which 
needs to be extrapolated from redshift z* ~ 1100 to redshift zero to obtain og. Ev- 
ery parameter that influences the expansion history of the Universe in this period 
affects the value of erg derived from the CMB alone. The relation between A s and 



ag is given by (e.g. Takada, Komatsu & Futamase 2004) 



a 2 8 (z) = A s 



2c 2 



00 I h 



k s dkD 2 (k,z)T 2 (k) 



k* 



n s — 1 



3 sin(A;_R) 3 cos(kR) 



(kRf 



(kRy 



(3.45) 



where R = 8/t" 1 Mpc, K = 0.02 Mpc" 1 and A s = (2.21 ± 0.09) x 10~ 9 Komatsu 



et al., 2009 



D(k,z) is the growth factor at redshift z and T{k) is the transfer 
function. If higher-order cosmological parameters such as the dark energy equation 
of state parameter w are marginalised over, then the measurements of fl m and erg 



weaken considerably (see Figure 3.11). 



We now assess the influence of the Hq prior. We replace the result of |Riess 
et al. [2011 with a measurement derived from the 6dFGS dataset using Baryon 
Acoustic Oscillations Beutler et al. 2011] . The prior from this study is lower than 
the former value and is given by Hq = 67.0 ± 3.2 km s -1 Mpc^ 1 . Using the 6dFGS 
value of # results in a 8 = 0.75 ± 0.13, Q m = 0.279 ± 0.028, b = 1.52 ± 0.29 and 
a p = 174 ± 106km/s. The quality of the fit is x 2 /d.o.f. = 1.35, very similar to 
the value obtained with the HST prior. Comparing the two results shows that a 
different prior in Hq shifts the constraint in as and b along the degeneracy shown 



in Figure |3.9| However, we note that the 6dFGS measurement of -f^o is derived 
from the same dataset as our present study and hence could be correlated with our 
measured growth rate. We use these results only for comparison, and include the 



values obtained using the HST prior in Table |3.1| as our final results of this section. 
Alternative methods for deriving as or fas at low redshift are provided by 



Gordon, Land & Slosar 


2007 


Abate & Erdogu 


2009 
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Figure 3.10: This plot shows the likelihood distribution in erg and Vt m for different 
cosmological probes. The solid blue contours show the 6dFGS result, the magenta 
dotted dashed contours shows the currently best result of weak lensing from the 
CFHT wide synoptic Legacy Survey Fu et al. , 2008 , the red solid contours show 
the result of the SFI++ peculiar velocity survey Nusser & Davis, 2011 , the green 



dashed contours show the result of Mantz et al. 2010 using cluster abundances and 



the black dotted contours are from WMAP7 Komatsu et al., 2011 



Nusser & Davis 


2011 


|Turnbull et al. 


2011 |Davis et al. 2011 


Hudson & Turnbull 



2012). 6dFGS will soon provide its own peculiar velocity survey of around 10 000 
galaxies. Velocity surveys have the advantage of tracing the matter density field 
directly, without the complication of a galaxy bias. However, they are much harder 
to obtain and current velocity surveys are 1-2 orders of magnitude smaller than 
galaxy redshift surveys. 
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6dFGS 

WMAP7.2 (ACDM) 
WMAP7.2 (wCDM) 
WMAP7.2 (oCDM) 
WMAP7.2 (owCDM) 




Figure 3.11: Likelihood distribution of <jg — Q m from WMAP7 for a ACDM model 
(solid black line), wCDM model (dashed red line), oCDM model (dotted dashed 
green line) and owCDM model (dotted magenta line). In blue we show the 6dFGS 
result. Both parameters are defined at redshift zero and WMAP constraints on these 
parameters depend on assumptions about the expansion history of the Universe. We 



use CosmoMC Lewis & Bridle, 2002 , together with the WMAP7.2 Komatsu et al. 



2011 dataset to produce these likelihood distributions. 



3.7 Cosmo logical implications 



In this section we test General Relativity by measuring the growth index 7. We 
would like to stress that the 7-parameterisation of modified gravity has its limitations 



and other more general parameterisations have been proposed (see e.g. Silvestri & 



Trodden|[2009| |Bean fc Tangmatitham||20T0l |Daniel fc Lmder]|20T0l |Clifton|[20lT 



Hojjati, Pogosian & Zhao 2011 Baker et al. 2011). However, this is going beyond 



the scope of this paper. 



We combine our result for ge(z c f() with the latest results from WMAP 7 Komatsu 
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Figure 3.12: Comparison of measurements of the growth of structure using galaxy 
surveys at different redshifts. The different data points belong to 6dFGS (solid blue 



circle, this paper), 2dFGRS (solid black circle; Hawkins et al. 2003[), SPSS (solid 
black boxes; Samushia et al.|20TT ), WiggleZ (solid black triangles; Blake et al.|2011a ) 
and VVDS (empty circle; Guzzo et al. 2008 ). We also included a WALLABY forecast 
with a 4% error-bar in red (see section 3.8). For the WiggleZ survey we also include 
the data points from Blake et al. 2011d| (empty triangles, shifted by Az = 0.005 



to the right for visibility), where the Alcock-Paczynski effect has been taken into 
account. We plot a ACDM model as well as a DGP model for comparison. 



et al. , 2011 , where we use the WMAP7.2 dataset provided on the NASA webpage^J 
While ACDM predicts 7 ~ 0.55, alternative theories of gravity deviate from this 
value. One example of such an alternative model is the DGP braneworld model 



of Dvali, Gabadadze & Porrati 2000 , in which our observable Universe is considered 



to be a brane embedded in a higher dimensional bulk space-time and the leakage of 
gravity force propagating into the bulk can lead to the current accelerated expansion 
of the Universe. Because of the missing dark energy component, this model predicts 
a larger growth index of 7 ~ 0.69 Linder, 2005 . 



In Figure 3.12 we compare measurements of the growth of structure gg from dif- 



http : //lambda. gsfc .nasa.gov/product/map/dr4/likelihood_get . cfm 
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ferent galaxy redshift surveys. For the WiggleZ survey we include data points which 
assume a correct fiducial cosmology (solid black triangles) as well as data points 
which account for the Alcock-Paczynski effect (empty black triangles). The degen- 
eracy between the Alcock-Paczynski effect and the linear redshift-space distortion 



signal increases the error by about a factor of two. As we showed in section 3.5.5 
the Alcock-Paczynski effect is very small in 6dFGS, which therefore yields a direct 
measurement of the redshift-space distortion signal. 

All data points seem to be in good agreement with the ACDM model (black 
solid line), while the DGP model generally predicts smaller values of gg. The value 
of <Tg for the two different models has been derived from the CMB scalar amplitude 



A s Komatsu et al. 2011 , where we use the corresponding Friedmann equation to 
calculate erg. 

The analysis method we apply in this section is summarised in the following 
four points: 



1. We produce a Monte Carlo Markov Chain (MCMC) with CosmoMC Lewis 



k Bridle| |2002] for a ACDM universe by fitting the WMAP7 dataset. The 
CMB depends on dark energy through the distance of last scattering and the 
late-time Integrated Sachs- Wolfe (ISW) effect. We avoid the contributions of 
the ISW effect by limiting the WMAP7 dataset to multipole moments £ > 100. 

2. Now we importance-sample the CosmoMC chain by randomly choosing a value 
of 7 in the range < 7 < 1 for each chain element. Since the value of cr 8 (z ef j) 
depends on 7 we have to recalculate this value for each chain element. First 
we derive the growth factor 



D(a eS ) = exp 



da' f(a')/a' 



"off 



(3.46) 



where a e g is the scale factor at the effective redshift a e g = l/(l+z e g). In order 
to derive o" 8i7 (-2 e fr) we have to extrapolate from the matter dominated region 
to the effective redshift, 

o"8 l7 (^eff) = D ^ ZeS } °s{zhi), (3.47) 
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where we use a%{zhi) from eq. 3.45 and zu = 50, well in the matter-dominated 
regime. 

3. We now calculate the growth rate using f 7 (z e g) ~ f2^(z cff ) and construct 



4. Finally we compare the model with gg = 0.423 ± 0.055 from Table 34^ and 
combine the likelihood from this comparison with the WMAP7 likelihood. 



The result is shown in Figure 3.13 Marginalising over the remaining parameters 
we get 7 = 0.547 ± 0.088 and Q m = 0.271 ± 0.027, which is in agreement with the 
prediction of a ACDM universe (7 rs 0.55). Our analysis depends only on the 
growth rate measured in 6dFGS and WMAP7. This makes our measurement of 
7 independent of systematic effects like the Alcock-Paczynski distortion which is a 
matter of concern for galaxy redshift surveys at higher redshift. 



3.8 Future low-redshift galaxy surveys: WALLABY 
and TAIPAN 



In this section we make predictions for the accuracy of fas measurements from 



future low-redshift galaxy surveys using a Fisher matrix analysis based on White, 



Song & Percival 2008 



The Wide-field ASKAP L-band Legacy All-sky Blind surveY (WALLABY)^ is 
an HI survey planned for the Australian SKA Pathfinder telescope (ASKAP), cur- 
rently under construction at the Murchison Radio-astronomy Observatory (MRO) 
in Western Australia. The survey will cover 75% of the sky and a proposal exists to 
fill up the remaining 25% using the Westerbork Radio Telescope. In this analysis we 



follow the survey parameters employed by Beutler et al. 2011] (see also Duffy et a!., 



in preparation), where for WALLABY we assume a An survey containing 600 000 
galaxies at a mean redshift of z = 0.04. The linear bias of a typical WALLABY 



galaxy is 0.7 Basilakos et al. , 2007 and the volume of the survey is 0.12/i Gpc 



3 http://www. atnf.csiro.au/research/ WALLABY 
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Figure 3.13: Likelihood distribution of 7 — Q m for our fit to g$ = 0.423 ± 0.055 from 



6dFGS and WMAP7 |Komatsu et^[20lT] . The 6dFGS contours (black) include 
a normalisation prior from WMAP7. Marginalising over the different parameters in 
the MCMC chain gives Vt m = 0.271 ±0.027 and 7 = 0.547 ±0.088. The low redshift 
of 6dFGS makes this measurement particularly sensitive to 7 and independent of 
systematic effects like the Alcock-Paczynski distortion. 



The TAIPAN survejj^] proposed for the UK Schmidt telescope at Siding Spring 
Observatory in New South Wales will cover a similar sky 6dFGS but will 

extend to a larger redshift such that z = 0.08. For our TAIPAN forecast we assumed 
the same sky-coverage as 6dFGS (/ s k y = 0.41), a bias of b = 1.4, a total of 400 000 
galaxies and a volume of 0.23/i~ 3 Gpc 3 . 

First we test the Fisher matrix prediction for fa 8 in the case of 6dFGS. We 
assume a survey volume of 0.08/i~ 3 Gpc 3 , with 81 971 galaxies. Using /c max = 
O.lh Mpc" 1 we forecast a measurement of /erg of 23%, while using /c max — 0.2/i Mpc -1 

4 TAIPAN: Transforming Astronomical Imaging surveys through Polychromatic Analysis of Neb- 
ulae 
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produces a 8.3% error. The actual error in fa 8 we found in this paper is 13%, 
somewhere between these two values. For WALLABY and TAIPAN we will report 
constraints for both A; max = 0.1/iMpc -1 and fc max = 0.2/tMpc -1 . 

With the specifications given above and using k m3iX = 0.1 (0.2)h Mpc -1 , the 
Fisher matrix forecast for WALLABY is a measurement of fa 8 with 10.5 (3.9)% er- 
ror. We included the WALLABY forecast with a 4% error-bar in Figure [3712] The 
model TAIPAN survey produces forecast errors of 13.2 (4.9)%, improving the results 
from 6dFGS by almost a factor of two. Although TAIPAN maps a larger volume of 
the Universe compared to WALLABY, it does not produce a better measurement of 
fas- WALLABY has a smaller galaxy bias, which increases the redshift-space dis- 
tortion signal by a factor of two compared to TAIPAN. This will also be very useful 



for breaking the degeneracy between bias and as using the technique of section 3.6.2 



The WALLABY survey will target galaxies rich in HI gas. Such galaxies will 
mostly populate under-dense regions of the Universe, because in groups and clus- 
ters galaxies are stripped of their gas by interactions with other galaxies and the 
intra-group and intra-cluster medium. This is the reason that Hl-selected galaxies 
possess a low bias (~ 0.7). However, this fact also implies that WALLABY galax- 
ies sample the density field in a manner that avoids high-density regions. These 
high-density regions are an important source of non-linear redshift-space distor- 
tions ( "finger-of-God" effect). We can hence suppose that non-linear effects will 
be smaller in amplitude in an HI survey compared to highly-biased surveys such as 
6dFGS (see e.g. |Simpson et al.||2011| or Figure 4 in |Reid fc White||2011[ ). This should 
allow the inclusion of much smaller scales in the analysis, producing more accurate 
measurements. A more detailed analysis using WALLABY mock catalogues is in 
preparation. 



McDonald & Seljak 2009 have suggested that multiple tracers within the same 



cosmic volume can be used to reduce the sampling variance and improve cosmological 



parameter constraints (see also Seljak 2009 Slosar 2009 Bernstein & Cai 2011). 
Using the ratio of the perturbation amplitudes of two surveys with different bias 
factors gives 

h + ffi 2 ab 2 + f[i 2 



b 2 + f^ 2 b 2 + ffi* ' 



(3.48) 
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where a = bi/b 2 . The angular dependence of this expression allows one to extract / 
without any dependence on the density field (see Figure 1 in Bernstein fc CaipOll ) . 
The density field is the source of the sampling variance error, since it will change, 
depending on the patch of the sky which is observed. Using the ratio of two tracers, 
the precision with which the growth rate / can be determined is (in principle) only 
limited by the shot noise, and not by the sampling variance. 

The three surveys discussed above, 6dFGS, WALLABY and TAIPAN, have a 
large overlapping volume which allows the potential application of this method. The 
technique works best for densely-sampled surveys with very different bias factors, b\ 
and b 2 . While TAIPAN and 6dFGS have very similar bias, the bias of WALLABY 
will be much smaller. 

We assume an overlap volume of 0.41 x 0.12/i -3 Gpc 3 = 0.049/i~ 3 Gpc 3 , where we 
multiply the sky coverage of 6dFGS and TAIPAN with the effective volume of WAL- 
LABY. For the different surveys we use the parameters as stated above. We forecast 
Fisher matrix constraints on fas of 10.3 (5)% using k mSuX = 0.1 (0.2)/iMpc _1 . Be- 
cause the overlap volume is only w 1/3 of the WALLABY volume, this result does 
not improve the measurement arising from WALLABY alone, especially considering 
that WALLABY itself may be able to include modes up to large fc max in the fitting 
process. 

Our results show that future surveys such as WALLABY and TAIPAN will pro- 
vide an accurate measurement of fas at low redshift, and will be able to complement 
future high-redshift surveys such as BOSS, which will have a similar accuracy for 
several data points over the higher redshift range 0.2 - 0.6 [Song fc Perched] [2009 



White, Song & Percival 2008 Reid & White 2011 



3.9 Conclusion 

In this paper we have measured the 2D correlation function of the 6dF Galaxy 
Survey. We derived a covariance matrix using jack-knife resampling as well as log- 
normal realisations and showed that both techniques give comparable results. We 
have modelled the 2D correlation function with a simple streaming model and a 
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more advanced approach suggested by Scoccimarro 2004 combined with the N- 



body calibrated results from Jennings, Baugh & Pascoli 2011a . We formulated 



these models in real-space including wide-angle corrections. For the final results on 



fas we chose the model by Scoccimarro 2004 , although we found that both models 



gave consistent results at sufficiently large scales. 

We analysed the measurement in two different ways. First we fitted for the two 
parameters ge(z e s) = f(z e s)a 8 (z e s) and gb(z e &) = ba 8 (z e g), where these constraints 
depend only on the 6dFGS data. Our second analysis method assumes a growth 
index from standard gravity (7 ~ 0.55) and fits for a 8 , b, Q m , H and a p , where we 
combine the 6dFGS measurement with a prior in the Hubble constant. All parameter 
measurements are summarised in TableHj] We can summarise the results as follows: 

• Our first analysis method found ge{z e s) = f{z e s)a 8 (z e s) = 0.423 ± 0.055 and 
gb{z e s) = ba 8 (z e ff) = 1.134 ± 0.073, at an effective redshift of z e s = 0.067. 
The 6dFGS measurement of gg, unlike high-redshift measurements, does not 
depend on assumptions about the expansion history of the Universe and the 
Alcock-Paczynski distortion. 

• In our second analysis method we used the angle dependence of redshift-space 
distortions in the 2D correlation function to break the degeneracy between 
the galaxy bias b and the normalisation of the matter clustering statistic a 8 , 
assuming standard gravity. We found a 8 = 0.76 ± 0.11, Q m = 0.250 ± 0.022, 
b = 1.48 ± 0.27 and a p = 174 ± 73km/s. This result uses a prior on H 



from Riess et al. 2011 . 



Combining our measurement of gg(z e s) with WMAP7 Komatsu et al. 2011 



allows us to measure the growth index 7, directly testing General Relativity. 
We found 7 = 0.547 ± 0.088 and Vt m = 0.271 ± 0.027, in agreement with the 
predictions of General Relativity (7 ~ 0.55). The 6dFGS measurement of this 
parameter is independent of possible degeneracies of 7 with other parameters 
which affect the correlation function at high redshift, such as the dark energy 
equation of state parameter w. 

We used a Fisher matrix analysis to forecast the constraints on fa 8 that would 
be obtained from two future low-redshift galaxy surveys, WALLABY and 
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TAIPAN. We found that WALLABY will be able to measure fa s to a forecast 
accuracy of 10.5% for k max = 0.1/iMpc -1 and 3.9% for /c max = 0.2/iMpc -1 . 
A combination of 6dFGS, TAIPAN and WALLABY, using the multiple-tracer 



method proposed by McDonald & Seljak 2009 , will be able to constrain /cr 8 
to 5 - 10.3%. These measurements would complement future large-volume 
surveys such as BOSS, which will measure the growth rate at much higher 
redshift (z > 0.2), and contribute to future precision tests of General Relativ- 
ity on cosmic scales. 
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Dependence of halo occupation on 
stellar mass 

[Beutler et al.\ 
MNRAS 429, 360 AB (2013) 

Abstract 

In this paper we study the stellar-mass dependence of galaxy clustering in the 6dF Galaxy 
Survey. The near-infrared selection of 6dFGS allows more reliable stellar mass estimates 
compared to optical bands used in other galaxy surveys. Using the Halo Occupation 
Distribution (HOD) model, we investigate the trend of dark matter halo mass and satellite 
fraction with stellar mass by measuring the projected correlation function, w p {r p ). We 
find that the typical halo mass (Mi) as well as the satellite power law index (a) increase 
with stellar mass. This indicates, (1) that galaxies with higher stellar mass sit in more 
massive dark matter halos and (2) that these more massive dark matter halos accumulate 
satellites faster with growing mass compared to halos occupied by low stellar mass galaxies. 
Furthermore we find a relation between M\ and the minimum dark matter halo mass 
(Mm in) of Mi ~ 22M m i n , in agreement with similar findings for SDSS galaxies. The 
satellite fraction of 6dFGS galaxies declines with increasing stellar mass from 21% at 
M stel i ar = 2.6 x 10 10 h- 2 M Q to 12% at M stellar = 5.4 x W 10 h~ 2 M indicating that high 
stellar mass galaxies are more likely to be central galaxies. We compare our results to 
two different semi-analytic models derived from the Millennium Simulation, finding some 
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disagreement. Our results can be used for placing new constraints on semi-analytic models 
in the future, particularly the behaviour of luminous red satellites. Finally we compare our 
results to studies of halo occupation using galaxy-galaxy weak lensing. We find good overall 
agreement, representing a valuable crosscheck for these two different tools of studying the 
matter distribution in the Universe. 



4.1 Introduction 



The first statistical studies of galaxy clustering Totsuji & Kihara, 1969 Peebles 



1973 Hauser & Peebles 1973 1974 Peebles 1974 found that the galaxy corre 



lation function behaves like a power law, which is difficult to explain from first 



principles Berlind & Weinberg, 2001 . More recent studies, however, found devia- 



tions from a power law. For example Zehavi et al. 2005 showed that the projected 



correlation function w p {r p ) of SDSS galaxies exhibits a statistically significant de- 
parture from a power law. They also showed that a 3-parameter Halo Occupation 



Distribution (HOD) model (e.g., |Jing, Mo fc Borner|[l998| |Ma fc Fry||2000[ [Pea- 
cock fc Smithl|2000l |Seljak||2000t |Scoccimarro et al.pOOTl |Berlind fc Weinberg||200T 



Cooray & Sheth 2002 ) together with a ACDM background cosmology, can account 



for this departure, reproducing the observed w p (r p ). 

Within the halo model the transition from the 1-halo term to the 2-halo term 
causes a "dip" in the correlation function at around 1 — 3/i _1 Mpc, corresponding 
to the exponential cutoff in the halo mass function. In case of a smooth transi- 
tion between the one- and two-halo terms, this can mimic a power-law correlation 
function. Studies with Luminous Red Galaxies (LRGs) found that the deviation 
from a power-law is larger for highly clustered bright galaxies Zehavi et al. , 2005[b 



Blake, Collister & Lahav, 2008 Zheng et al., 2009 Zehavi et al. 2010, and at 



high redshift Conroy, Wechsler & Kravtsov, 2006 , which agrees with theoretical 



predictions Watson et al., 2011 



While galaxy clustering is difficult to predict, dark matter clustering is domi- 
nated by gravity and can be predicted for a given cosmology using N-body simula- 
tions. Using models for how galaxies populate dark matter halos, which are usually 
motivated by N-body simulations, we can directly link galaxy clustering and matter 
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clustering. This can be modelled in terms of the probability distribution p(N\M) 
that a halo of virial mass M contains N galaxies of a given type. On strongly non- 
linear scales the dark matter distribution is given by the actual density distribution 
of the virialized halos, while on large and close to linear scales the dark matter 
distribution can be predicted from linear perturbation theory. 

HOD modelling has been applied to galaxy clustering data from the 2-degree 



Field Galaxy Redshift Survey (2dFGRS) [Porciani, Magliocchetti fe Norberg 2007 



Tinker et al. , 2006] and the Sloan Digital Sky Survey (SDSS) |van den Bosch, Yang 



Tinker et al. , 


& Mo, 


2003 



fc Mo[ [2003j IMagliocchetti fc Porciaru] [20031 |Zehavi et alj |2005|E1 |Tinker et al. 



2005 Yang et al. 2005, 2008 Zehavi et al., 2010 . More recently it also became 



possible to model the clustering of high-z galaxies using VVDS Abbas et al. , 2010 



Bootes Brown et al. , 2008 , DEEP2 Coil et al. , 2006 and Lyman-break galaxies at 



high redshift in the GOODS survey Lee et al. , 2006 . Such studies revealed, that the 
minimum mass, M m j n for a halo to host a central galaxy more luminous than some 
threshold, L is proportional to L at low luminosities, but steepens above L*. Massive 
halos have red central galaxies with predominantly red satellites, while the fraction of 



blue central galaxies increases with decreasing host halo mass. Furthermore Zehavi 



et al. 2005b found that there is a scaling relation between the minimum mass of 
the host halos, M min and the mass scale, Mi of halos that on average host one 
satellite galaxy in addition to the central galaxy, M\ w 23 M m ; n . Using a different 



HOD parameterization, Zheng, Coil & Zehavi 2007 found the relation to be M\ ~ 
18M min , very similar to Zehavi et al. 2010 who found M% ~ 17M min . 

The 6dF Galaxy Survey is one of the biggest galaxy surveys available today 
with a sky coverage of 42% and an average redshift of z — 0.05. The survey includes 



about 125 000 redshifts selected in the J, H, K, bj, rp-bands Jones et al. , 2004 , 2005 



2009 . The near infrared selection, the high completeness and the wide sky coverage 



make 6dFGS one of the best surveys in the local Universe to study galaxy formation. 
This dataset has been used to study the large scale galaxy clustering to measure the 



Hubble constant using Baryon Acoustic Oscillations Beutler et al. , 2011 , as well as 
the growth of structure at low redshift Beutler et al. , 2012 . While these previous 
studies used the i^-band selected sample, in this analysis we use the J-band. The 
J-band allows the most reliable stellar mass estimate of the five bands available in 
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6dFGS, because of its lower background noise. Together with the bj — r> colour we 



can derive stellar masses using the technique of Bell & De Jong 2001 , which leads 
to a dataset of 76 833 galaxies in total. The photometric near-infrared selection from 
2MASS makes the stellar mass estimates in 6dFGS more reliable than stellar mass 



estimates in other large galaxy surveys which rely on optical bands Drory et al. 



20041 |Kannappan fc Gawiser] [20071 |Longhetti fc Saraccoj [20091 |Grillo et al.[ [2008 



Gallazzi & Bell 2009 



Numerical N-body simulations are usually restricted to dark matter only. To 
understand galaxy formation, baryonic effects such as feedback and gas cooling, 
have to be included. Such simulations face severe theoretical and numerical chal- 
lenges. Semi-analytic models build upon pre-calculated dark matter merger trees 
from cosmological simulations and include simplified, physically and observationally 
motivated, analytic recipes for different baryonic effects. Semi-analytic models have 
been shown to successfully reproduce observed statistical properties of galaxies over 



a large range of galaxy masses and redshifts (e.g. Croton et al. 2006 Bower et al. 



20061 |De Lucia fc Blaizot|2007||Bertone et al.|2007l|Font et al.|2008~l |Guo et al.|2010[ ) 



and allow a level of understanding unavailable in N-body simulations. The underly- 
ing models are necessarily simplified and often use a large number of free parameters 
to fit different observations simultaneously. In this paper, we derive 6dFGS mock 



surveys from semi-analytic models based on the Millennium Simulation Springel et 



al, 2005 and compare the properties of these surveys with measurements in 6dFGS. 



Our results can be used to improve upon these semi-analytic models and further our 
understanding of baryonic feedback processes on galaxy clustering. 



Mandelbaum et al. 2006 studied halo occupation as a function of stellar mass 



and galaxy type using galaxy-galaxy weak lensing in SDSS. They found that for 
a given stellar mass, the halo mass is independent of morphology below M ste ii a r = 
10 11 Mq, indicating that stellar mass is a good proxy for halo mass at that range. 



We compare our results with Mandelbaum et al. 2006 which represents a valuable 
crosscheck of the HOD analysis using two very different techniques, weak lensing 
and galaxy clustering. 

This paper is organised as follows: First we introduce the 6dF Galaxy Survey in 



section 4.2 together with the technique to derive the stellar masses. We also explain 
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how we derive our four volume-limited sub-samples in stellar mass and redshift, 



which are then used for the further analysis. In section 4J3 we calculate the projected 
correlation function, w p (r p ), for each sub-sample and use jack-knife re-sampling to 
derive the covariance matrices. In section |4.4 we fit power laws to the projected 
correlation functions of the four sub-samples. In section 4J3 we introduce the HOD 
framework and in section |4~6 we apply the HOD model to the data. In section 4.7 



we derive 6dFGS mock samples from two different semi-analytic models, which we 



then compare to our results in section 4.8 together with a general discussion of our 



findings. We conclude in section 4.9 



Throughout the paper we use r to denote real space separations and s to denote 
separations in redshift space. Our fiducial model to convert redshifts into distances 
is a flat universe with Q m = 0.27, w = —1 and Qk — 0. The Hubble constant is set to 
Hq = 100/ikms~ 1 Mpc _1 which sets the unit of stellar masses to h~ 2 M & , while most 
other masses are given in h~ 1 M & . The HOD model uses cosmological parameters 



following WMAP7 Komatsu et al. 2011 



4.2 The 6dF Galaxy survey 



The 6dF Galaxy Survey [6dFGS; | Jones et al] |2004] |2005| [2009] is a near-infrared 
selected (J, H, K) redshift survey covering 17000 deg 2 of the southern sky. The J, H 
and K surveys avoids a ±10° region around the Galactic Plane to minimise Galactic 
extinction and foreground source confusion in the Plane. The near-infrared photo- 
metric selection was based on total magnitudes from the Two-Micron All-Sky Sur- 



vey Extended Source Catalog [2MASS XSC; Jarrett et al. , 2000] . The spectroscopic 
redshifts of 6dFGS were obtained with the Six-Degree Field (6dF) multi-object spec- 
trograph of the UK Schmidt Telescope (UKST) between 2001 and 2006. 

The 6dFGS J-selected sample used in this paper contains (after completeness 
cuts) 76 833 galaxies selected with 9.8 < J < 13.75. We chose the J-band because 
it has the highest signal-to-noise of the three 2MASS bands. While there is slightly 
less extinction in the iT-band compared to the J-band, for practical purposes, the 
J-band has better S/N because the night sky background glow is much less in the 
J- than in the i^-band. The near infrared selection makes 6dFGS very reliable for 
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stellar mass estimates. 

The mean completeness of 6dFGS is 92 percent and the median redshift is z = 
0.05. Completeness corrections are derived by normalising completeness-apparent 
magnitude functions so that, when integrated over all magnitudes, they equal the 
measured total completeness on a particular patch of sky. This procedure is outlined 



in the luminosity function evaluation of Jones et al. 2006 and also in Jones et al., 



(in prep). The original survey papers Jones et al. 2004, 2005, 2009 describe in full 
detail the implementation of the survey and its associated online database. 

The clustering in a galaxy survey is estimated relative to a random (unclus- 
tered) distribution which follows the same angular and redshift selection function 
as the galaxy sample itself. We base our random mock catalogue generation on the 
6dFGS luminosity function Jones et al. , 2006 , where we use random numbers to 
pick volume-weighted redshifts and luminosity function-weighted absolute magni- 
tudes. We then test whether the redshift-magnitude combination falls within the 
6dFGS J-band faint and bright apparent magnitude limits (9.8 < J < 13.75). We 
assigned a bj-rp colour to each random galaxy using the redshift- bj-rp colour re- 
lation measured in the data and used these to derive stellar masses for the random 



galaxies using the same technique as for the actual galaxies (see section 4.2.1). 



4.2.1 Stellar mass estimate and volume-limited sub-samples 



To calculate the stellar mass for our dataset we use the stellar population synthesis 



results from Bruzual & Chariot 1993 together with a scaled Salpeter initial mass 



function (IMF) as reported in Bell & De Jong 2001 



log 10 (M stellar /Lj) 
log 10 (Lj) 

log 10 (M ste llar) 



-0.57C, 



+ 0.48 



(M} un -Mj)/2.5 
log 10 (M ste ii ar /Lj) + log 10 (Lj) 



(4.1) 



with the 2MASS bj-r F colours, rp the J-band absolute magnitude, Mj and the 



J-band absolute magnitude of the sun, M} un = 3.70 Worthey, 1994 . The biggest 
uncertainty in stellar mass estimates of this type is the choice of the IMF. Assuming 
no trend in IMF with galaxy type, the range of IMFs presented in the literature cause 
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Figure 4.1: The distribution of 6dFGS galaxies in log stellar mass and redshift. The 
redshift and stellar mass cuts imposed to create the four volume-limited samples 
(S\ to 5*4 see Table 4.1) are shown by the coloured lines. All galaxies in the upper 
left quadrant created by the two correspondingly coloured lines are included in the 
volume-limited sub-samples. The plot shows a randomly chosen set of 20% of all 
galaxies. 



uncertainties in the absolute normalisation of the stellar M/L-ratio of a factor of 
2 in the near- infrared Bell & De Jong 2001 . The 6dFGS stellar mass function as 
well as a comparison of different stellar mass estimates is currently in preparation 
(Jones et al. in prep). 

We create four volume-limited sub-samples in redshift and stellar mass. This is 
done by choosing an upper limit in redshift (z max ) and then maximising the number 
of galaxies by choosing a lower limit in stellar mass (M^J^), meaning that every 
galaxy above that stellar mass will be detected in 6dFGS, if its redshift is below 
the redshift limit (see Figure 4.1). Because of the distribution in bj-rp colour, a 
clear cut in absolute magnitude does not correspond to a clear cut in stellar mass 



133 



4.2. THE 6DF GALAXY SURVEY 



0.007 
0.006 
0.005 
0.004 



n. 0.003 
c 



0.002 
0.001 

0, 



_ 












L/vv 






J 






2 





























0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 



Figure 4.2: Galaxy density as a function of redshift for the four different volume- 
limited sub-samples S% - S4. All samples follow an approximately constant number 



density indicated by the dashed lines and listed in the last column of Table 4.1 



and therefore our samples are not perfectly volume-limited. We use the absolute 
magnitude limit which corresponds to a chosen redshift limit and derive a stellar 
mass limit, using the bj-rp colour corresponding to the 50% height of the bj-rp 
distribution (bj-rp(50%) = 1.28). We create four sub-samples (S1-S4) with upper 
redshift cuts at z max = 0.05, 0.06, 0.07 and 0.08 and with the mean log stellar 
masses ranging from log 10 (M stollar //i~ 2 M ) = 10.41 to 10.73. We also include a low 
redshift cut-off at z min = 0.01. Since the stellar mass distributions for these samples 
overlap, especially for the higher stellar mass samples, the results are correlated 



to some extent (see Figure 4.1). All our sub-samples are summarised in Table 4.1 



Figure |472] shows the galaxy density as a function of redshift for the four sub-samples. 
The roughly constant number density with redshift shows that our sub-samples are 
close to volume-limited. 
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4.3 Data analysis 

In order to analyse the galaxy clustering in all three dimensions, we calculate the 
co-moving distances from us to each galaxy, 

°c = ± I' 4L> («) 



H„J E(z _ 
where z is the measured redshift, and 

E(z)=[n™(l + zT + nf] 1/2 . (4.3) 

We assume a flat universe with = and f2| d = 1 — and treat dark energy 
as a cosmological constant (w fid = — 1). Given the low redshift of our dataset, these 
assumptions have no impact on the final results. 

We measure the separation between galaxies in our survey along the line of sight 
(it) and perpendicular to the line of sight (r p ) and count the number of galaxy pairs 
on this two-dimensional grid. We do this for the 6dFGS data catalogue, a random 
catalogue with the same selection function, and a combination of data-random pairs. 
We call the pair-separation distributions obtained from this analysis step DD(r p , 7r), 
RR(r p , 7r) and DR(r p , n), respectively. In the analysis we used 30 random catalogues 
with the same size as the real data catalogue and average DR(r p , n) and RR(r p , it). 



The random mocks are sampled from the 6dFGS luminosity function (Jones et al. 



2006 and also Jones et al, in prep.), and hence they contain the same evolution of 



luminosity with redshift that we see in 6dFGS itself. The redshift-space correlation 



function is then given by the Landy & Szalay 1993 estimator 



i+ RR(r p ,n) [nj RR(r p ,n) \n d > ' 



where the ratio n r /rid is given by 



i W 



(4.5) 



n d J2f - Wj 

and the sums go over all random (N r ) and data (Nd) galaxies. Here we employ a 
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completeness weighting, Wi, where we weight each galaxy by the inverse sky- and 
magnitude completeness at its area of the sky (Jones et al., in prep.). 

From the two-dimensional correlation function, £(r p ,7r), we calculate the pro- 
jected correlation function 



w p (r p ) = 2 d7r£(r p ,7r), (4.6) 
Jo 

where we bin £(r p ,7r) in 30 logarithmic bins from 0.1 to 100/i _1 Mpc in r p and it. 



The upper integration limit in eq. 4.6 was chosen to be 7r max = 50 h 1 Mpc for all 
sub-samples. 

To derive a covariance matrix for the projected correlation function we use the 
method of jack-knife re-sampling our galaxy samples. First we divide the dataset 
into iV = 400 subsets, selected in R.A. and Dec. Each re-sampling step excludes one 
subset before calculating the correlation function. The covariance matrix is then 
given by 

c l3 = Yl K(^;) - K(^) - . ( 4 - 7 ) 

k=l 

where w h (r p ) is the projected correlation function estimate at separation r p with the 
exclusion of subset k. The mean value is defined as 



N 

w(r, 



v> N 

k=i 



1 N 



In Figure 4.3 we show the four correlation matrices for the different volume-limited 



sub-samples, which is defined as 

r« = (4-9) 
V 

with Cij being the covariance matrix. The highly correlated part above r p m 0.1 
- lh^ 1 Mpc is caused by a cosmic variance dominated error, while on small scales 
the Poisson dominated error leads to almost uncorrelated bins. With increasing 
redshift, the number density decreases, which shifts the transition from the Poisson 
noise dominated part to the sample variance limited part to larger scales. 
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Figure 4.3: The correlation matrix for the projected correlation functions, w p (r p ) 
derived from the covariance matrices for the four volume-limited 6dFGS sub-sample 
(Si - S4) used in this analysis. The colour indicates the correlation with blue indi- 
cating no correlation and red indicating high correlation. 



138 



CHAPTER 4. DEPENDENCE OF HALO OCCUPATION ON STELLAR MASS 




Figure 4.4: The angular correlation function of 6dFGS (blue data points) as well 
as the target catalogue (black solid line). The deviation between the two angular 
correlation functions is caused by fibre proximity limitations. The lower panel shows 
the ratio of the two correlation functions, which can be used as a weight to correct 
for the fibre proximity limitations (see section 4.3.1). The open data points show 
the corrected 6dFGS angular correlation function which is in very good agreement 
with the target catalogue. 



We also note that wide-angle effects can be neglected in this analysis, since we 
are interested in small-scale clustering (see Beutler et al. |2011| and Beutler et al. 



2012 for a detailed investigation of wide-angle effects in GdFGS). 



4.3.1 Fibre proximity limitations 



The design of the 6dF instrument does not allow fibres to be placed closer than 
5.7arcmin Jones et al. , 2004 , which corresponds to a distance of r„ pa 0.3/iT 1 Mpc 



at redshift z = 0.07. This limitation is relaxed in 6dFGS where about 70% of the 
survey area has been observed multiple times. However, for the remaining 30% we 
have to expect to miss galaxy pairs with a separation smaller than 5.7arcmin. 
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Figure 4A shows the angular correlation function for the 6dFGS redshift cata- 
logue (w z , blue data points) and the target catalogue (w p , solid black line). While the 
target catalogue contains 104 785 galaxies, the redshift catalogue of galaxies which 
have a J-band, 6j-band and rp-band magnitude contains 76 833. The dashed line 
indicates the angular scale of the fibres. The two angular correlation functions agree 
on scales 9 > 0.1°, with the redshift catalogue falling below the target catalogue at 
lower scales. 



In the lower panel of Figure 4.4 we show the ratio (1 + w p )/(l 



w 7 



The 



dashed line shows a spline fit to the blue data points, which than can be used to 
up-weight galaxy pairs with small angular separations and correct for the fibre prox- 
imity effect Hawkins et al. , 2003 . This weighting is additional to the completeness 
weighting Wi we introduced earlier. While the completeness weighting is applied to 
single data- as well as random galaxies, the fibre proximity weighting is only applied 
to data galaxy pairs. Applying this weight to the 6dFGS redshift catalogue results 



in the open data points in Figure 4.4 which now agree very well with the target 
catalogue. 



4.4 Power law fits 



The projected correlation function can be related to the real-space correlation func- 



tion £(r), using Davis & Peebles, 1982 



w p (r. 



POO 

p) = 2/ dy{[(rl+yr /2 } 
Jo 



r dr ^(r)(r 2 — r'i) 1 ^ 2 . 



(4.10) 



If the correlation function is assumed to follow a power law, £(r) = (r/ro) 7 , with 
the clustering amplitude r and the power law index 7, this can be written as 



w p {r p ) = r p 



r 



(4.11) 
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Table 4.2: Best fitting parameters, r and 7 for power law fits to the projected 
correlation functions w p (r p ) of the four volume-limited 6dFGS sub-samples (Si - 
£4). The fitting range in all cases is 0.1 < r p < 40 /i _1 Mpc with 24 bins and 2 free 
parameters. The last column shows the reduced \ 2 indicating the goodness of the 
fit. 



sample 


r [/i -1 Mpc] 




7 


X 2 /d.o.f. 




Si 


5.14 ±0.23 


1. 


849 ± 0.025 


23.5/(24- 2) = 


1.07 


So 


5.76 ±0.17 


1. 


826 ±0.019 


19.5/(24-2) = 


0.89 




5.76 ±0.16 


1. 


847 ±0.019 


32.6/(24- 2) = 


1.48 


Si 


6.21 ±0.17 


1. 


846 ± 0.019 


35.8/(24- 2) = 


1.63 



with T being the Gamma-function. Using this equation we can infer the best-fit 
power law for £(r) from w p (r p ). 



Table 4.2 summarises the results of the power law fits to the four 6dFGS sub- 



samples. The best fitting power laws are also included in Figure 4.5 together with 
the four projected correlation functions. We set the fitting range to be 0.1 < r p < 
40 /i -1 Mpc which includes 24 bins. Although 6dFGS has very good statistics at 
scales smaller than 0.1/i -1 Mpc we do not use them for our fits, since at such scales 
the fibre proximity correction becomes more than 30% (see section 4.3.1 and Fig- 
ure 



4.4). 



The value of the clustering amplitude, ro, increases from 5.14 to 6.21/i -1 Mpc 
with increasing stellar mass, while there doesn't seem to be a clear trend in 7, 
which varies around the value of 7 = 1.84. The reduced x 2 i n the last column of 
Table |4.2 indicates a good fit to the data for the first two sub-samples but grows to 
X 2 /d.o.f. = 1.48 and 1.63 for the high stellar mass sub-samples, indicating deviations 
from a power law. The lower panels of Figure 4.5| show the different projected 
correlation functions divided by the best fitting power law (blue data points). Here 
we can see that the deviations show systematic patterns. Such effects could be 



related to the strong correlations between bins as seen in Figure 4.3 Nevertheless, 
these patterns can be addressed with a full HOD analysis, which we will pursue in 
the next section. We will compare the result of our power law fits with other studies 
in a more detailed discussion in section 14.81 
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Figure 4.5: The projected correlation function for the four different volume-limited 
6dFGS sub-samples (Si - S , 4 ). The dashed black lines show the best fitting power 



laws (see Table 4.2) while the solid black lines show the best fitting HOD models, 



derived by fitting the data between 0.1 < r p < 40h 1 Mpc. The lower panel shows 
the data and HOD models divided by the best fitting power laws. 
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We remark at this point, that there is no satisfying theoretical model that 
predicts a power law behaviour of the correlation function, and hence the motivation 
of such a fit can only be empirical. While alternative approaches (e.g. the HOD 
model) rely on assumptions about the clustering behaviour of galaxies and dark 
matter, they are physically motivated and hence allow us to learn more about galaxy 
clustering than a pure empirical power law fit. 



4.5 Theory: Halo Occupation Distribution 



The Halo Occupation Distribution (HOD) model describes the relation between 
galaxies and mass in terms of the probability distribution p(N\M) that a halo of 
virial mass M contains N galaxies of a given type. Knowing how galaxies populate 
dark matter halos, we can use a dark matter correlation function and infer the galaxy 



correlation function. We use CAMB Lewis & Bridle, 2002 to derive a model matter 



power spectrum which we turn into a correlation function using a Hankel transform 



1 

2^2 



dk p(k)k 2S ^l. (4.12) 
kr 



The underlying cosmological model is fixed to (Qbh 2 , Q c h 2 , n s , as) = (0.02227, 
0.1116, 0.966, 0.8) as reported in Komatsu et al. |2011| (we set erg = 0.9 for one 
special case). We have to be aware of the fact that the fitted HOD parameters 
depend somewhat on the assumed values of Q m and as and hence the absolute 
values of the HOD parameters could be biased, if the assumed cosmology is wrong. 
However in this study we focus on the relative HOD parameters for different stellar 
mass selected sub-samples, which is fairly robust against such uncertainties. 



Here we employ an analytic HOD methodology that is similar to that of |Zheng7 

We 



Coil & Zehavi 2007 ; Blake, Collister & Lahav 2008 and Zehavi et al. 2010 



utilise analytic approximations for the halo mass function Tinker et al. , 2008 , the 



biased clustering of halos Tinker et al. , 2005 , the profile of dark matter within 



halos Navarro, Frenk & White, 1996 , the concentration of dark matter halos as 



a function of mass Duffy et al. 2008 and the nonlinear dark matter power spec- 



trum Smith et al. , 2003 . In the next section we will discuss the basic inputs for 
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our HOD model before discussing the HOD formalism itself. 



4.5.1 Halo mass function and the halo-mass relation 

The spatial density of halos as a function of their mass M is specified by the halo 
mass function 



dn 



f(°)Pr. 



dln(cr) 



(4.13) 



d\n(M) J v /r " L dM 

where a is the variance of the linear density field, when it is smoothed with a top- 
hat filter of mass scale M = ^7iR 3 p m with the matter density p m = Q m p c and 
p c = 3Hq/8tiG being the critical density of the Universe. G is the gravitational 
constant. 

In linear theory we can calculate a by 

D(z) 2 



a 2 (M,z) 



2tt 2 



P{k)W 2 M {k)k 2 dk, 



(4.14) 



where P(k) is a linear model power spectrum at redshift z = 0, D(z) is the linear 
growth factor normalised to unity at z = and Wj^ik) is the Fourier transform of 
the real-space top-hat window function at radius R 



W M {k) 



(kRf 



[sm(kR) - kRcos(kR)} 



(4.15) 



To a first approximation, the function f(a) is universal, and the effects of power 
spectrum shape, redshift evolution, and background cosmological model (e.g. Q m 
and enter only through determining \d\na/dM\ and p m . Using N-body simu- 



lations Tinker et al. 2008 found 



+ i 



-c/a 2 



(4.16) 



with the best fitting values A = 0.1086, a = 1.47, b = 2.57, c = 1.19 for z = halos, 
defined to be spherical regions centred on density peaks enclosing a mean interior 
over-density of 200 times the cosmic mean density p m . 

For a given peak height v = 8 c /a(M, z) (with 5 C = 1.686 being the linear over- 
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density at which a spherical perturbation collapses), we can also calculate the halo 



bias b(M), where we use the fitting function proposed by Tinker et al. 2010 

b(M) 



1- A- 



v 



Bv b + Cu c 



with A = 1.0000597, a = 0.13245, B = 0.183, b = 1.5, C = 0.26523 and c = 2.4 



(4.17) 



4.5.2 Halo density profile 

N-body simulations provide many predictions about structure formation, one of 
which is the principal prediction that dark matter halos should have a cuspy density 



profile, on average well represented by the two-parameter NFW profile Navarro, 

Ps 



Frenk & White, 1996 



p(r,M) = — 

{r/r s )(l + r/r s y 

where r s is the characteristic scale radius and p s is given by 



M 



Ps 



Anr 3 s [ln(l + c) - c/(l + c)] 



(4.18) 



(4.19) 



with the concentration parameter c = c(M, z). The profile is truncated at the virial 
radius, which for a given halo mass M is defined by 



/ 3M y /3 

\4ttAp J 



(4.20) 



with A = 200. 



The relation between the concentration parameter, the core radius r s and the 
virial radius r v j r is c = r v ; r /r s . In this analysis we use the halo concentration mass 



relation from Duffy et al. 2008 



c(M,z) = 6.71 



M 



M, 



pivot 



-0.091 



[l + Z 



,-0.44 



(4.21) 



with Mpi vot = 2 x 10 h M & . We tried replacing eq. 4.21 with the form suggested 



by Bullock et al. 1999 , and found that the best-fitting HOD parameters changed 
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by much less than the statistical errors. 



4.5.3 HOD framework and formalism 



In the HOD parametrisation it is common to separate the clustering contributions 
from the most massive galaxies, which are assumed to sit in the halo centre, from 
satellite galaxies. This picture of how galaxies populate halos is supported by hydro- 



dynamic simulations (e.g. Berlind et al. 2003 Simha et al. 2009) and semi-analytic 
models (e.g. White & Frenk 1991; Kauffmann, White & Guiderdoni 1993; Croton 



et al. 2006 Bower et al. 2006). The mean central and satellite number density of 



galaxies that populate dark matter halos of mass M is Zheng et al. , 2005 



(N C (M)) 
(N S (M)) 



if M < M n 



1 if M > M min 



(4.22) 



' M_ 
Mi 



where M min is the minimum dark matter halo mass which can host a central galaxy, 
Mi corresponds to the mass of halos that contain, on average, one additional satellite 
galaxy ((iVs(Mi)) = 1) and a sets the rate at which halos accumulate satellites when 
growing in mass. In very massive halos the number of satellites is proportional to 
halo mass M to the power of a. The total HOD number is given by 



(N t {M)) = (N c (M))[l + (N s (M))) 



(4.23) 



so that a dark matter halo can only host a satellite galaxy if it contains already a 
central galaxy. In our model we assume a step like transition from (N C (M)) = to 
(N C (M)) = 1. In reality this is more likely to be a gradual transition with a certain 
width 0"i O gM More et al. , 2009 . To account for this, other studies (e.g. Zehavi et al. 



2010) modify the HOD parametrisation of the central halo term to 
(N c (M)) = l I 



crf / log 10 (M)-log 10 (M mir 



ClogM 



(4.24) 
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which turns into eq. 14.221 for the case <j\ 



ogM 



0. For our data we found that the 



reduction in \ 2 obtained from fits including a\ og M as an additional parameter is not 
big enough to justify this parameterisation. For our largest sub-sample (S2) we 
found Ax 2 = —0.78, while a new parameter would be justified when A% 2 < —2 
following the Akaike information criterion Akaike et al. , 1974 



There are also higher parameter models for the satellite fraction such as 



(N S (M)) 



M -M 
Mi 



(4-25) 



where Mo is the minimum halo mass at which satellites can exist. Again we tested 
this model and found that the reduction in \ 2 does not justify this additional pa- 
rameter. Hence we chose the two parameter model (Mi, a) of eq 



4.22 



Within the halo model we can account separately for the clustering amplitude of 
galaxies which sit in the same dark matter halo (one halo term) and galaxies which 
sit in different dark matter halos (two halo term) . At small scales the clustering will 
be dominated by the one halo term and at large scales it will be dominated by the 
two halo term. For the correlation function this can be written as 



(4.26) 



where £ih( r ) and £ 2 /i ( r ) represent the one halo and two halo terms respectively. 



4.5.4 The 1-halo t 

We separate the 1-halo term into contributions from central-satellite galaxy pairs 
and satellite-satellite galaxy pairs. The central-satellite contribution is given by 

where n g is the galaxy number density and p(r, M) is the halo density profile (see 
eq. 4.18). The lower limit for the integral is the virial mass M v i r (r) corresponding 



1 -M m in is fixed by the number density as will be explained in section 4.5.5 
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to the virial separation r vir (see eq. 4.20). 

The satellite-satellite contribution is usually given by a convolution of the halo 
density profile with the halo mass function. Here we calculate this term in k-space 
since a convolution then turns into a simple multiplication 

P~(k) = ±J" dM?^N e (M)N?(M)\{k, M)\ (4.28) 

where X(k, M) is the normalised Fourier transform of the halo density profile p(r, M) 
which is given by 

X(k, M) = ^ f °° drp(r, M)r 2 ^^. (4.29) 
JV1 In kr 



The real-space expression of eq. (4.28) is 



«'' s<r) = h r* p ^ w * 2 ^' (43o) 

which can then be combined with the central-satellite contribution to obtain the 
1-halo term 

6ft(r) = *£*(r)+£r(r). (4-31) 



4.5.5 The 2-halo term, £ 2h (r) 

The 2-halo term, £2/1 ( r ); can be calculated from the dark matter correlation func- 
tion since on sufficiently large scales the galaxy and matter correlation function are 
related by a constant bias parameter. We calculate the two halo term in Fourier 
space as 



P2h(k,r) 



P m (k) x 

M Um (r) 



dM \{k,M) 



dM 



n'Jr) 



(4.32) 



where P m (k) is the non-linear model power spectrum from CAMB including halofit [Smith 



et al. , 2003| , Mu m (r) is the halo mass limit for which we can find galaxy pairs with a 
separation larger than r v i r and bh(M, r) is the scale dependent halo bias at separation 
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r for which we assume the following model Tinker et al. 2005 

b 2 h (M,r 



b\M) 



[l + 1.17£ ro (r)] 



1.49 



[l + 0.69£ m (r)] 



2.09 ' 



(4.33) 



where £m(r) is the non-linear matter correlation function with b(M) given by eq. 4.17 
We define the restricted galaxy number density n' g (r) at separation r 



n'(r) 



o 



M m m *m Ntm . 

dM K ' 



(4.34) 



For smaller separations r, only very small mass halos contribute to the 2-halo term 
because high mass halos have a virial radius r vir which is larger than the galaxy sep- 
aration, and hence the galaxy sits in the same halo. This mass limit is already given 
by Mu m (r) and can be derived using the "n^-matched" approximation described 
by Tinker et al. 2005 . First we calculate the restricted number density using 



n ,2 (r) 



dM 2 d ^^N t (M 2 )P(r, M u M 2 ) 



(4.35) 



where p(r, Mi, M2) quantifies the probability of non-overlapping halos of masses 
Mi and M2 with separation r. Defining x = r/(r v i r; i + r^p), where r v i rj i and 
r vir 2 are the virial radii corresponding to masses Mi and M 2 (see eq. 4.20). Using 
y — (x — 0.8)/0.29, p(y) is then given by 



p(y) 



1 


3y 2 



2y 3 



if y > 1 
if y < 
otherwise. 



(4.36) 



Given the value of n' (r), we increase My, m (r) in equation 4.34 to produce a number 



density which matches the measured n g . This value is then used in equation 4.32 
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We obtain the correlation function for the 2-halo term with 

^ (r) = h r dk P2h{K r)k2% -^ 1 - (4 - 37) 

Now we have to correct this term from the restricted galaxy density to the entire 
galaxy population which leads us to the final two halo term 

f n! (r) \ 2 

l + 6,(r)=(^J + (4-38) 

So far our HOD model has three free parameters: M m j n , Mi and a. However we use 
the galaxy number density to fix M min for each set of Mi and a via 



n, = / dM^-N t (M) (4.39) 
*M^ Nt( M, ,4,0, 

lv *xaxn 

We can now calculate the real-space correlation function, £(r) = [1 + Cift.( r )l +£,2h(r) 
for any set of Mi and a, which we then turn into the projected correlation function, 
w p (r p ), using eq. 4.10 



4.5.6 Derived quantities 

Since the HOD model is directly based on a description of dark matter clustering 
and its relation to galaxy clustering, an HOD model can tell us much more about a 
galaxy population, than just the two parameters Mi and a. 

For X(k,M) —7- 1 (which corresponds to large separations r in real-space) the 
two halo term simplifies to 

P 2h (k,r)^b 2 cS P m (k), (4.41) 



where the effective bias is the galaxy number weighted halo bias factor Tinker et 



al. , 2005 



b cS = — r dM d ?^b h (M)N t (M). (4.42) 
n g J dM 
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We can also ask what is the average group dark matter halo mass for a specific set 
of galaxies (often called host-halo mass). Such a quantity can be obtained as 

M cff = — f° dM^j^MN t (M), (4.43) 
n g J Q dM 

which represents a weighted sum over the halo mass function with the HOD number 
as a weight. 

The averaged ratio of satellite galaxies to the total number of galaxies is given 

by 

SZ°dM**®.N e (M)[l + N.(M)] 
and the central galaxy ratio is given by f c = 1 — / s . 



4.6 HOD parameter fits 

To compare the HOD models to the data we use the same fitting range as for the 
power law fits earlier (0.1 < r p < 40 h^ 1 Mpc). The free parameters of the fit are 
Mi and a as described in the last section and hence we have the same number of 
free parameters as for the power law fits. 



All the fitting results are summarised in Table 4.3 and Figure 4.5 The reduced 
X 2 in the last column of Table 4.3 indicates a good fit to the data in all cases and 



the ratio of data to best fitting power law in the lower panels of Figure 4J3 shows 
that the HOD model reproduces the double peak structure present in the data. We 
also note that the reduced \ 2 m case °f the HOD fits is uniformly lower than for 
the power law fits indicating a better fit to the data for all sub-samples. 



Figure 4.6 shows the 2D probability distributions in Mi and a for the four 
different volume-limited sub-samples. There is a strong trend of increasing Mi with 
stellar mass and a weaker but still significant trend of increasing a with increasing 
stellar mass. This indicates that dark matter halos that host a central galaxy with 
higher stellar mass have their first satellite on average at a larger dark matter halo 
mass. However the number of satellites increases more steeply with halo mass for 
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Figure 4.6: Two dimensional proba- 
bility distribution for log 10 (Mi) and 
a for the HOD fits to the four 6dFGS 
volume-limited sub-sample (Si - S4). 
We include two fits for the largest 
sample 5*2, where for S 1 ^ 8-0 ' 9 we 
change our standard assumption of 
a s = 0.8 to er 8 = 0.9. The parameters 
derived from the fits are summarised 
in Table |4.3| The best fitting values 



are marker with black crosses. 
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Figure 4.7: The effective galaxy bias (see eq. 4.42) as a function of log stellar mass 
for 6dFGS (blue data points). The increase in galaxy bias with stellar mass can be 
described by eq. |4.45| which corresponds to the black line. 



samples with higher stellar mass. This trend indicates that halos of higher mass have 
greater relative efficiency at producing multiple satellites. Similar trends were found 
in sub-samples of SDSS galaxies by Zehavi et al.||2005b As a test for the sensitivity 
of the fitting results to the upper fitting limit we performed a fit to sub-sample S 2 
with the fitting range 0.1 < r p < 20h~ 1 Mpc. All parameters agree within la with 



the fit to the larger fitting range. We called this fit S' 2 and included it in Table 4.3 



Table |4.3| also includes derived parameters like the minimum dark matter halo 
mass, M m j n , the satellite fraction, f s , the effective dark matter halo mass, M e g and 
the effective galaxy bias, 6 e fr- The errors on these parameters are calculated as 
the 68%-confidence level of their ID probability distribution. The effective dark 
matter halo mass seems to be almost constant for all samples while the minimum 
dark matter halo mass M min increases with stellar mass. The satellite fraction / s 
decreases with increasing stellar mass, indicating that galaxies with high stellar mass 
have a higher probability to be central galaxies. The increasing effective galaxy bias 
indicates that galaxies with higher stellar mass are more strongly clustered and 
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hence reside in high density regions of the Universe. The increasing effective galaxy 
bias is well described by a power law (in log stellar mass) of the form 

Hollar) = (1-05 + M stellar /M*) X (0.8/<r 8 ) (4.45) 

with M* = 1.18 x 10 ll h~ 2 M Q . We compare this function to the measurements in 



Figure 4.7 The absolute stellar masses are subject to significant uncertainties and 
different methods to derive stellar masses can come to very different conclusions. 
However, most stellar mass estimates are related by a simple constant offset, and 



relation 4.45 can be scaled accordingly. The equation above is only valid for the 



stellar mass range probed in this analysis (M ste n ar = 2.6 — 5.4 x 10 10 h 2 M ) since 



the underlying dynamics are most likely not captured in eq. |4.45 

We will discuss the implications of all these results in the next sections. First 
we will derive 6dFGS mock samples using different semi-analytic models. We will 
then compare the predictions from these semi-analytic models to our data followed 
by a comparison to other studies. 



4.7 semi-analytic mock catalogues 



To compare our result with theory, we derive 6dFGS mock catalogues from two dif- 



ferent semi-analytical models Croton et al. , 2006 Bower et al. , 2006 , both based 



on the Millennium Simulation jSpringel et al 2005 publicly available through the 
Millennium Simulation database^] Semi-analytic models are based on an under- 
lying N-body simulation together with theoretically and observationally motivated 
descriptions of gas cooling, star formation and feedback processes. 

The Millennium Simulation is a dark matter only N-body simulation which 
traces the hierarchical evolution of 2160 3 particles in a periodic box of 500 3 /i~ 3 Mpc 3 
from redshift z = 127 to z — 0. The underlying cosmological model follows WMAP1 
cosmology Spergel et al. , 2003 given by a matter density of Q m = Q dm + Q b = 0.25, 
a cosmological constant of Q\ = 0.75, a Hubble constant of H = 75km/s/Mpc, a 



2 http://gavo.mpa-garching.mpg.dc/MyMillennium/ 
http://galaxy-catalogue.dur. ac.uk:8080/Millennium/ 
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spectral index of n s = 1 and a r.m.s. of matter fluctuations in 8k' 1 Mpc spheres of 
<j$ = 0.9. The individual particle mass of the simulation is 8.6 x 10 8 h~ 1 M & and halos 
and sub-halos are identified from the spatial distribution of dark matter particles 



using a standard friends-of-friends algorithm and the SUBFIND algorithm Springel 



et al, 2001 . All sub-halos are then linked together to construct the halo merger 
trees which represent the basic input of the semi-analytic models. 

Here we are using the z = output of the Millennium Simulation. To ensure 
that the stellar masses are calculated in a consistent and comparable way, we follow 
the following procedure to derive the 6dFGS mock catalogues: 

1. We apply the 6dFGS i\~-band apparent magnitude limits of 8.85 < K < 12.75 
to the full 500/i~ 3 Mpc 3 simulation box. We have to use the i^-band instead 
of the J-band, which is actually used in this analysis, because none of the 
semi-analytic models provide J-band magnitudes. However, the 6dFGS J- 
band and K-bsnd samples have significant overlap, meaning that almost all 
galaxies which have a i^-band magnitude also have a J-band magnitude. We 
also account for sky- and magnitude incompleteness. 

2. We re-calculate stellar masses for each galaxy using the technique described 
in section 4.2.1 but with the corresponding JY-band relations Bell & De Jong 
2001 , instead of J-band. The re-calculation of the stellar masses ensures 
that potential disagreement with our measurement is not caused by a different 
technique of deriving stellar masses or a different assumptions about the IMF. 

3. We apply the same redshift and stellar mass limits to the semi-analytic cat- 
alogues, which we used to produce the four volume-limited samples in the 



6dFGS dataset (see Table 4.1) 



We found that all mock 6dFGS catalogues derived from these semi-analytic models 
contain fewer galaxies than the data sample (by about 40%). This could be related 
to the slightly different cosmology used in these simulations, which should have its 
largest impact on large clusters, which are sampled in 6dFGS. We are not attempting 
to correct for such differences in the cosmological model. The aim of this part of 
our analysis is to test the current predictive power of semi-analytic models. 
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Semi-analytic models are often grouped into "Durham models" and "Munich 



models". The Bower et al. 2006 model belongs in the group of "Durham models" 



while the Croton et al. 2006 model belongs in the group of "Munich models". 



4.7.1 Durham models 



In the Durham models, merger trees are produced following Helly et al. 2003 which 



are independent of those generated by Springel et al 2005 . When the satellite galaxy 
falls below a certain distance to the central galaxy given by -R m er g e = r c + r s , where 
r c and r s are the half mass radii of the central and satellite galaxy, respectively, 
the satellite and central galaxy are treated as one. The largest of the galaxies 
contained within this new combined dark matter halo is assumed to be the central 



galaxy Benson et al., 2002 , whilst all other galaxies within the halo are satellites. 



The dynamical friction and tidal stripping which are present in such a system are 
modelled analytically. These models are based on NFW density profiles for the 
central halo as well as the satellite halo, while galaxies are modelled as a disc plus 
spheroid. 

4.7.2 Munich models 



The Munich models are based on the original merger trees by Springel et al 2005 



One of the key differences between these merger trees and the ones used in the 
Durham models is that the Munich models explicitly follow dark matter halos even 
after they are accreted onto larger systems, allowing the dynamics of satellite galax- 
ies residing in the in-falling halos to be followed until the dark matter substructure 
is destroyed. The galaxy is than assigned to the most bound particle of the sub-halo 
at the last time the sub-halo could be identified. 

In the Munich models, a two-mode formalism is adopted for active galactic nuclei 
(AGNs), wherein a high-energy, or "quasar" mode occurs subsequent to mergers, and 
a constant low-energy "radio" mode suppresses cooling flows due to the interaction 
between the gas and the central black hole Croton et al. , 2006 . In the quasar 



model, accretion of gas onto the black hole peaks at z ~ 3, while the radio mode 
reaches a plateau at z ~ 2. AGN feedback is assumed to be efficient only in massive 
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halos, with supernova feedback being more dominant in lower-mass halos. 



The galaxies in the Croton et al. 2006 model can be of three different types: 
central galaxies (type 0), satellites of type 1 and satellites of type 2. Satellites of 
type 1 are associated with dark matter substructures, which usually refers to recently 
merged halos. Satellites of type 2 are instead galaxies whose dark matter halo has 
completely merged with a bigger halo and are not associated with a substructure. 
We treat both type 2 and type 1 as satellite galaxies. 



4.7.3 Testing semi-analytic models 



In their original paper Bower et al. 2006 compare the K-band luminosity function, 



galaxy stellar mass function, and cosmic star formation rate with high-redshift ob- 
servations. They find that their model matches the observed mass and luminosity 
functions reasonably well up to z ~ 1. Kitzbichler & White 2007 compare the 



magnitude counts in the bj, r>, I and .fT-bands, redshift distributions for i\~-band 
selected samples, bj- and i^-band luminosity functions, and galaxy stellar mass func- 



tion from Croton et al. 2006 and the very similar model by De Lucia & Blaizot 
[2007 with high redshift measurements. They find that the agreement of these mod- 
els with high-redshift observations is slightly worse than that found for the Durham 
models. In particular, they find that the Munich models tend to systematically 
overestimate the abundance of relatively massive galaxies at high redshift. 



Snaith et al. 2011 compared four different semi-analytic models De Lucia & 



Blaizot 


2007 


Bower et al. 


2006; 



vations and found that all models show a shallower, wider magnitude gap, between 
the brightest group galaxy and the second brightest, compared to observations. 



de la Torre et al. 2010 compared measurements of VVDS with the model by De 
Lucia & Blaizot 2007 . They found that the model reproduced the galaxy clustering 
at z > 0.8 as well as the magnitude counts in most bands. However the model 
failed in reproducing the clustering strength of red galaxies and the bj - I colour 
distribution. The model tends to produce too many relatively bright red satellites 
galaxies, a fact that has been reported in other studies as well (over-quenching 



problem: Weinmann et al. 2006 Kimm et al. 2009 Liu et al. 2010). 
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Figure 4.8: This plot shows the friends-of-friends halo mass log 10 (Mf o f) of the satel- 
lite and central galaxies from the Bower et al. 2006 semi-analytic catalogue with 
the 6dFGS selection. The bimodal distribution is caused by a dominant fraction 
of halos, in which the 6dFGS galaxies are the central galaxy and a smaller fraction 
with larger halo mass, in which the 6dFGS galaxies are satellite galaxies. 



4.8 Discussion 



4.8.1 Effective halo mass and satellite fraction 



As is evident from Table 4J3 the 6dFGS galaxies sit in massive central dark matter 
halos and most of our galaxies are central galaxies in these halos with only a small 
fraction being satellite galaxies. Because of the way 6dFGS galaxies are selected, 
the large majority are red elliptical galaxies and hence our findings agree very well 
with previous studies, which also found that such galaxies are strongly clustered and 
therefore must reside in high density regions. 

First we will discuss the effective dark matter halo mass, which represents the 
effective group or cluster mass for the distribution of galaxies (not the mass of the 
individual halos which host the galaxies). This quantity appears almost constant 
for our four different sub-samples. The increasing central galaxy mass (M min ) with 
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Table 4.4: Summary of parameters derived from different semi-analytic models. We 
impose the 6dFGS f^-band apparent magnitude selection (J-band is not available 
for these semi-analytic models) as well as correct for incompleteness. The stellar 



masses are re-calculated using the technique described in section |4.2.1[ To calculate 
the effective dark matter halo mass we used the friends-of-friends halo mass. The 
stellar mass and the effective halo mass are calculated as the mean of the distribution. 



sample 



Croton et al. 



SI 

S2 
S3 
S4 



20061 



fs log 

0.138 
0.099 
0.085 
0.090 



10 



(M s{ 



h-' 2 M e 



log 



10 



10.38 
10.48 
10.57 
10.67 



(Bower et al. , 2006 



Meg 
h- 1 M (3 



13.400 
13.389 
13.507 
14.013 



SI 
S2 
S3 
S4 



0.28 
0.25 
0.22 
0.18 



10.36 
10.47 
10.58 
10.68 



13.695 
13.661 
13.674 
13.777 



stellar mass is offset by the decreasing satellite fraction, resulting in a fairly constant 
M e ff. While the HOD model allows us to derive only averaged parameters for each 
sample, semi-analytic models directly connect dark matter halo masses with single 



galaxies. Figure 4.8 shows the distribution of central and satellite galaxies as a 



function of friends-of-friends (fof) halo mass derived from the Bower et al. 2006 



semi-analytic catalogue together with the 6dFGS selection criteria. The catalogue 
contains fewer satellites than central, but the satellites sit in very massive dark 
matter halos and hence have a significant impact on M e g. This plot shows that 
most 6dFGS galaxies sit in 10 11 - 10 12 h~ l M & halos, while the satellites sit in very 
massive groups and clusters of up to 10 15 /i~ 1 M o . While the median would be 
around 10 12 /i _1 M , the effective mass M e g from the HOD model is the mean of this 
distribution, which is pushed to very large values by the satellite fraction. 

We calculated the effective halo mass M e g = (Mf f) for the four volume-limited 
sub-samples derived from the two semi-analytic models and summarise these results 
in Table 4.4| While the effective halo mass appears constant in case of the Bower et 



al. |2006 model, the Croton et al. 2006] model shows an increase with increasing 
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stellar mass. These different behaviours for the two different semi-analytic models 
are most likely connected to the different trend in the satellite fraction shown in 
Figure 4.9 While the Bower et al. 2006 model shows a constant decrease in the 



satellite fraction, the Croton et al. 2006 model reaches a constant at large stellar 
mass which causes the effective halo mass to rise. A lower satellite fraction at 
a fixed stellar mass means that at fixed halo mass the satellite galaxies are less 
massive 



Weinmann et al. 2006 Kimm et al., 2009 Liu et al., 2010 



From Figure 4.9 we can see that both the Croton et al. 2006 and the Bower et 



al. 2006 model give slightly smaller stellar masses than measured in 6dFGS. This 



indicates that the semi-analytic catalogues contain too few bright galaxies at low 
redshift. 



Figure 4.9 also includes the fit to sample S% which assumes <j$ = 0.9 (red data 
point), which agrees with the assumptions of the Millennium Simulation, while the 
blue data points assume erg = 0.8. A larger <7g causes a larger satellite fraction which 
than causes a larger effective halo mass. This effect seem to bring our results closer 
to the Bower et al. 2006 model. 



4.8.2 The Mi — M min scaling relation 



The HOD analysis of early data from SDSS by Zehavi et al. 2005b showed the 
relation, M\ ~ 23 M min , between the mass of halos that on average host one ad- 
ditional satellite galaxy, M\ and the minimum dark matter halo mass to host a 
central galaxy, M min . This has been confirmed in subsequent studies (Mi 18 M min 
Zheng, Coil k Zehavi||2007| and M x « 17 M min iZehavi et al.||2010l). This relation 



in 



implies that on average a halo hosting two galaxies of the type studied in their anal- 
ysis has a mass ~ 23 times the mass of a halo hosting only one galaxy of the same 
type. 



Zehavi et al. 2010 also found that this scaling factor is somewhat smaller at 



the high luminosity end, corresponding to massive halos that host rich groups or 
clusters. This latter trend likely reflects the relatively late formation of these massive 
halos, which leaves less time for satellites to merge onto central galaxies and thus 
lowers the satellite threshold M\. 
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Figure 4.9: Satellite fraction of the 6dFGS stellar mass sub-samples (blue) and the 



corresponding semi-analytic sub-samples based on the Millennium simulation Cro- 



ton et al. , 2006 Bower et al. , 2006 . The red data point represents a fit to sample 
S2 which assumes as = 0.9, while the blue data points assume as = 0.8. 



Theoretical studies of HODs in dark matter simulations iKravtsov et al. . 2004 



and those predicted by Smoothed-particle hydrodynamics (SPH) and semi-analytic 
galaxy formation models Zheng et al. , 2005 reveal a similar relation with a scaling 
factor of ~ 20. The large gap between Mi and M min arises because in the low 
occupation regime, a more massive halo tends to host a more massive central galaxy, 



rather than multiple smaller galaxies Berlind et al., 2003 



Abbas et al. 2010 did a similar study using data from the /-band selected 
VIMOS-VLT Deep Survey (VVDS) and found the ratio Mi/M min to be « 40 - 
50. The galaxies in this study are at much higher redshift (z ~ 0.83) compared 
to the SDSS galaxies. This result means that in order to begin hosting satellite 
galaxies, halos sampled by the VVDS survey need to accrete a larger amount of 
mass compared to SDSS halos. 



Wake et al. 2011 also studied high redshift galaxies (1.1 < z < 1.9) in the 
NEWFIRM Medium Band Survey (NMBS) and found the ratio Mi/M min to be 
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Figure 4.10: The relation between Mi and M m i n for the four different volume- 
limited 6dFGS samples (blue data points), compared to |Zehavi et al. |2005b| (black 



2010 



et al. 



solid squares), iBlake, CoUister & Lahav |2008| (black open triangles), Zehavi et al. 



(black open squares), Abbas et al. 2010] (black solid triangles), Matsuoka 



2011] (black open points), Wake et al. |2011 



Leauthaud et al. |2012| (magenta solid triangles). Zehavi et al. 2005b 



(magenta open triangles) and 
don't report 



errors on their parameters and the 6dFGS errors are smaller than the data points. All 
studies which investigate HOD as a function of stellar mass are coloured (including 
our analysis), while luminosity defined studies are in black. 



« 4 — 10. This, together with the result by Abbas et al. 2010 shows, that this 
relation is not as fundamental as originally thought, but strongly dependent on the 
type of halos probed in each analysis. 



Leauthaud et al. 2012 studied galaxy clustering in the COSMOS survey using 



threshold stellar mass samples. Using the equations discussed in Leauthaud et al. 



2011 we can derive M 1 and M min from their HOD parametrisation, which is included 



in Figure 4.10 



Matsuoka et al. 2011 analysed ~ 60 000 massive (log lo (M stellar //i- 2 M ) > 10.7) 



galaxies from the UKIRT Infrared Deep Sky Survey (UKIDSS) and the SDSS II Su- 
pernova Survey. This analysis shows a very different clustering amplitude depending 
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on whether the observed and theoretical number densities are matched up or not. 
This dependency makes a comparison of our results with their derived parameters 
very difficult. Nevertheless, the ratio Mi/M min in their analysis does not depend 
significantly on their initial assumptions and hence we included their results in Fig- 
Their results indicate a lower Mi/M min ratio at large Mi consistent 



urc 



4.1C 



with Zehavi et al. 2010 and Blake, Collister & Lahav 2008 



In our analysis we found a scaling relation of Mi « 22 M min (the exact values 
are 22.29 ± 0.39, 21.81 ± 0.30, 22.08 ± 0.28 and 22.22 ± 0.25 for S^S^, respectively). 
In Figure |4.10 we compare our results with other stellar mass selected samples 
(coloured data points) and with luminosity threshold selected samples (black data 
points). The ratio Mi/M min found in 6dFGS is in agreement with Zehavi et al. 



2005b and slightly larger than Zehavi et al. 2010 . While Zehavi et al. 2005b 



and Zehavi et al. 2010 study the same type of galaxies, Zehavi et al. 2005b 



uses 



an HOD parameterisation very similar to ours, while Zehavi et al. 2010 



uses a 



parameterisation based on 5 free parameters, suggesting that the differences might 
be related to the parameterisation. 



4.8.3 Comparison to Mandelbaum et al. |2006] 



An alternative method for probing the connection between stellar mass and halo 



mass is galaxy-galaxy weak lensing. Here we will compare our findings to Mandel- 
baum et al. 2006 , who used weak lensing of ~ 350 000 galaxies from SDSS and 



looked at the dependence of the amplitude of the lensing signature as a function 
of galaxy type and stellar mass. They used stellar masses derived from the z-band 



magnitude and the ratio M/L z from Kauffmann et al. 2003 with the assumption 



of a Kroupa 2001 IMF 



The HOD model employed in their analysis is given by 



N t {M) = N c {M) + N s {M) 



(4.46) 



with the central galaxy number given by eq. 4.22 while the satellite galaxies number 



3 We use their results, in which they match the observed and theoretical number densities, since 
this agrees with our method. 
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Figure 4.11: The effective dark matter halo mass of the four 6dFGS sub-samples 
(blue data points) compared to early type (black triangles) and late type (black 
squares) galaxies from Mandelbaum et aL] [2006] . In red we also include the fit to 
sample £2 which assumes <jg = 0.9 in agreement with the assumptions of Mandel- 
baum et al. 2006 , while the blue data of assume as = 0.8. 



is modelled by a step-like function 



N S {M) 




if M > 3M min 

if M min < M < 3M n 

if M < M min . 



(4.47) 



The normalisation constant k can be determined by matching the measured satellite 
fractions /™ a and is given in appendix [El Here we have adjusted the nomenclature 



to the one used in our analysis. We also note that Mandelbaum et al. 2006 assumed 
erg = 0.9, while we assumed as = 0.8 for most of our fits. 



From the equation above we can derive the effective halo mass using eq. 4.43 
together with the dark matter halo mass function. 



We compare the reported satellite fraction, /™ a , of Mandelbaum et al. 2006 and 



the derived effective halo mass, M e R , with our results in Figure 4.11 and Figure 4.12 
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Figure 4.12: The satellite fraction of the four 6dFGS sub-samples (blue data points 
compared to early type (black triangles) and late type (black squares) galaxies 
from Mandelbaum et aL] [2006] . In red we also include the fit to sample S2 which 
assumes <7g = 0.9 in agreement with the assumptions of Mandelbaum et al. 2006 
while the blue data points assume erg = 0.8. 



Figure 4.11 shows that the galaxies sampled by Mandelbaum et al. 2006 have a very 



similar effective halo mass compared to 6dFGS galaxies. The 6dFGS measurements 
follow the trend of early-type galaxies in these measurements, while the late types 
have a slightly smaller effective halo mass. However, the errors in case of the results 



by Mandelbaum et al. 2006 don't allow one to distinguish between the early and 



late type galaxies. Since effective halo mass and satellite fraction are strongly linked, 



we also compare the satellite fraction of Mandelbaum et al. 2006 with our results 



in Figure 4.12 Both figures also include the fit to sample S2 where we assumed 
erg = 0.9, instead of the standard <jg = 0.8, since this agrees with the assumption 
m 



Mandelbaum et al. 2006 . A larger <jg increases the effective halo mass by a small 



amount and hence does not impact this comparison significantly. 



Overall we see very good agreement between our results and the results of Man- 



delbaum et al. 2006 , which is reassuring since the two techniques are subject to 



different systematic uncertainties. We can also emphasise that the 6dFGS results 
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are extremely precise compared to the lensing results. 



4.8.4 Comparison to other studies 



Meneux et al. 2007 study galaxy cluster mass dependence on stellar mass in the 



VIMOS-VLT Deep Survey (VVDS) at redshift 0.5 < z < 1.2. The stellar masses 
in this sample cover the range log lo (M ste n ar //i~ 2 M ) = 8.7 - 10.7. To quantify the 
clustering, they use power laws fits to the projected correlation function w p (r p ) and 
found an evolution in ro from 2.76 ± 0.17k- 1 Mpc at log lo (M steUar //i~ 2 M ) > 8.7 to 
r = 4.28 ± 0.45k- 1 Mpc at log lo (M stcUar //i- 2 M ) > 10.2. The slope changes over 
the same range from 7 = 1.67 ± 0.08 to 7 = 2.28 ± 0.28. Comparing the results 



by Meneux et al. 2007 with Li et al. 2006 who used the clustering of 200 000 SDSS 



galaxies at z = 0.15, showed that the evolution of the amplitude and shape of the 
correlation function w p (r p ) with redshift is faster for low stellar mass objects than 
for high stellar mass objects. At low stellar mass, the amplitude of w p {r p ) increases 
by a factor ~ 2 - 3 from high to low redshift, while at the high stellar mass range 
log 10 (M s teiiar/^~ 2 ^©) = 10.2 - 10.7, the amplitude at z ~ 0.85 and z ~ 0.15 are very 
similar, within the error bars. In 6dFGS we found a significantly larger clustering 
amplitude, ranging from r = 5.14 ± 0.23/i~ 1 Mpc at log lo (M ste ii ar //i~ 2 M ) = 10.41 
to r = 6.21 ± 0.17/i -1 Mpc at log 10 (M st eiiar/^~ 2 M3) — 10.73, while our power law 
index 7 stays constant at 7 ~ 1.84. This shows that the trend observed in Meneux 



et al. 2007 and Li et al. 2006 continues to the 6dFGS redshifts, even for high 



stellar mass galaxies. 



4.9 Conclusions 

We present in this paper an analysis of the clustering properties of four stellar mass 
selected volume-limited sub-samples of galaxies from the 6dF Galaxy Survey. The 
stellar masses are calculated using the J-band magnitude and the bj — r F colour 



following the technique by Bell & De Jongl [2001 1 . The average log-stellar mass for 



the different sub-samples ranges from log 10 (M ste n ar //i 2 M ) = 10.41 — 10.73. Our 
analysis has the following main results: 
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The projected correlation function, w p (r p ), for the two low stellar mass sub- 
samples (Si and S2) can be described by a power law with an acceptable x 2 > 
while the two high stellar mass sub-samples (S3 and S4) have a reduced \ 2 °f 
1.48 and 1.63, respectively. We also find patterns in the deviations between 
the best power law fit and the data, which can naturally be explained within 



the halo model. This is in agreement with theoretical studies Watson et al. 



2011 , predicting that the disagreement of w p {r p ) with a power law fit should 



grow with clustering amplitude. 

We used an HOD parameterisation with three free parameters (Mi, a, M min ), 
representing the typical halo mass, satellite power law index and minimum 
host halo mass, respectively. The minimum halo mass is fixed by the galaxy 
number density for all our parameter fits. We performed fits to the projected 
correlation functions of the four different volume-limited sub-samples. We 
tested alternative HOD parameterisations, but found that our data does not 
justify more free parameters. Mi and a show increasing trends for increasing 
stellar mass, with log 10 (Mi//i -1 M ) ranging from 13.4 — 14 and a from 1.21 — 
1.4. This means that galaxies with larger stellar mass populate larger dark 
matter halos, which accrete satellites faster with increasing mass, compared 
to dark matter halos populated by galaxies with lower stellar mass. 

From the halo model, we can derive averaged parameters for the four samples 
such as the satellite fraction, effective dark matter halo mass and the effective 
galaxy bias. We found that the satellite fraction decreases with stellar mass 
from 21% at log 10 (M stellar //i- 2 M ) = 10.41 to 12% at log 10 (M stcllar //^ 2 M ) = 
10.73. The effective dark matter halo mass stays constant at log 10 (M e g / h~ 1 M Q ) 
13.55 for all four sub-samples. The effective galaxy bias increases with stel- 
lar mass indicating that galaxies with higher stellar mass reside in denser 
environments. The increase in the effective galaxy bias can be described by 
(1.05 + M ste i lar /M*) x (o-g/0.8) with M* = 1.18 x lO n /i- 2 M . 



We use two semi-analytic models, based on the Millennium Simulation (Croton 



et al.||2006 and Bower et al.||2006 ) to derive 6dFGS mock surveys. We compare 



the results of these semi-analytic models with our measurements. The |Croton 
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et al. 2006 model under-predicts the satellite fraction, while the Bower et al. 



2006 model is in better agreement with our observations, although it slightly 



over-predicts the satellite fraction. Since the effective dark matter halo mass is 



strongly linked to the satellite fraction, the Bower et al. 2006 model prediction 



of M e ff is again in better agreement with our observations. 6dFGS allows a 
powerful test of semi-analytic models, because of (1) the robust stellar mass 
estimates and (2) the focus on "red satellites", which semi-analytic models 
struggled to reproduce in the past. Our results can be used as a new constraint 
on semi-analytic models in the future. 

For the four volume-limited samples we find a constant scaling relation be- 
tween Mi and M min of M\ ~ 22M min , in agreement with studies of SDSS 



galaxies Zehavi et al.j 2005b ; 2010 . This indicates that 6dFGS galaxies pop- 



ulate dark matter halos in a similar way to SDSS galaxies. However, we see a 



wide spectrum of this ratio ranging from > 40 in Abbas et al. 2010 to < 10 
in 



Matsuoka et al. 2011; Wake et al. , 2011 indicating that this ratio is not 



universal for all types of galaxies. 



We compare our results with the results of Mandelbaum et al. 2006 from 



galaxy-galaxy weak lensing. We find overall good agreement which represents 
a valuable crosscheck for these two different clustering measurements. Al- 
though our analysis depends on slightly more assumptions and covers a smaller 
range in stellar mass, the 6dFGS results are extremely precise compared to 
the lensing results. 
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Chapter 5 

Summary and Conclusion 



In this thesis we analyse data from the 6-degree Field Galaxy Survey. The three 
main results of this thesis are: 

1. We used the baryon acoustic oscillation signal in 6dFGS at the effective redshift 
of z e g = 0.106 to constrain the Hubble constant H with a small dependence on 
other cosmological parameters. We found H = 67.0±3.2kms _1 Mpc" 1 , which 
represents one of the best constraints on this parameter ever obtained, with 
very different systematic uncertainties compared to more standard methods 
like the distance ladder technique. 

2. The imprint of redshift-space distortions on the two-dimensional correlation 
function measures the amplitude of the gravitational force acting between 
galaxies and therefore allows us to test General Relativity. We combine our 
measurement of the parameter combination /cr 8 = 0.423 ±0.055 with WMAP- 
7 to constrain the gravitational growth index 7 = 0.547 ± 0.088. Our result is 
in very good agreement with the prediction of General Relativity (7 ~ 0.55). 

3. We use the clustering properties of 6dFGS galaxies as a function of stellar 
mass to explore how they populate dark matter haloes and hence place new 
constraints on semi-analytic models. 

In chapter [2] we use the baryon acoustic oscillation signal to make an absolute 
distance measurement to the effective redshift of 6dFGS at z e s = 0.106. We measure 
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a distance of Dv{z e s) = 456 ± 27Mpc and a distance ratio of r s (zd)/ Dv(z e s) = 
0.336 ± 0.015, where r s (zd) is the sound horizon at the drag epoch z^. An absolute 
distance measurement at such a low redshift constraints mainly the Hubble constant 
Hq, almost independently of other cosmological parameters. Our measurement of 
the Hubble constant is Hq = 67.0±3.2kms _1 Mpc -1 (4.8% precision), which is 1.7a 
smaller than the latest measurement reported using the distance ladder technique 
(H = 73.8 ± 2.4km/s/Mpc; Riess et al. 2011). This tension has been confirmed 



by several other publications (e.g. Mehta et al. 2012). While this is an interesting 
potential indication of systematic errors in either technique, we have to keep in mind 
that a 1.7a deviation corresponds to a probability of ~ 10% and hence is statistically 
not significant. Our result represents a competitive and independent alternative to 
Cepheids and low-z supernovae in measuring the Hubble constant, which depends 
on very different systematic errors and hence allows a very interesting test of the 
current standard cosmological model. 

Instead of r s (z d )/ D v (z e s) we can constrain the parameter A(z eS ) = 100D v (z e s)\/Ti 
which can be used to directly constrain the matter density fl m . We found fl m = 
0. 287+0. 039(l+w) + 0.039f2fc±0. 027 where the weak dependence on the dark energy 
equation of state parameter w and the curvature parameter Q k occurs because our 
measurement is not exactly at redshift zero. 

The 6dFGS BAO detection occupies a redshift very different to other galaxy 
surveys which report such a detection, and hence represents a new data point on 
the BAO Hubble diagram. While at low redshift we are not directly sensitive to 
dark energy, 6dFGS can help to break degeneracies arising when fitting the CMB 
data alone which then leads to improved constraints on dark energy. We show 



that combining the SDSS DR7 BAO measurement from Percival et al. 2010 with 



6dFGS and WMAP7 Komatsu et al. 2011 results in a constraint on the dark 
energy equation of state parameter of w = —0.97 ± 0.13, which represents a 24% 
improvement compared to the measurement without 6dFGS. 

The 6dFGS BAO result has been combined with other datasets in several subse- 



Blake et al. 


2011c 


Mehta et al. 


2012 


Anderson et al. 


2012 



Sanchez et al. 2012). 
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In chapter [3] we use the redshift-space distortions signal in the 2D correlation 
function of 6dFGS, £(r p ,7r), to measure the parameter combination f(z e s)as(z e s) = 
0.423 ± 0.055, where / ~ is the growth rate of cosmic structure and a 8 is 

the r.m.s. deviation of matter fluctuations in 8h~ l Mpc spheres. The parameter 
7 incorporates the gravitational force and is predicted to be 7 ~ 0.55 in General 
Relativity with a ACDM background cosmology. Measuring this parameter allows 
us to test models of Gravity. Since only the growth rate / is directly related to 
7, we have to combine our measurement of fas with another dataset, to break the 



degeneracy between / and as- We combine our result with WMAP7 Komatsu et 



al. 2011 , which yields 7 = 0.547 ± 0.088, consistent with the prediction of General 



Relativity. Because of the low redshift of 6dFGS, this measurement is independent 
of the fiducial cosmological model (Alcock-Paczynski effect). We also show that our 
result is not sensitive to the model adopted for non-linear redshift-space distortions 
or the fitting range. 

Furthermore we discuss a method to directly constrain / and erg with low redshift 
data only, where we combine 6dFGS with a prior on the Hubble constant from |Riess 
With this technique we find / 



et al. 2011 



^1 r 



0.58 ± 0.11, a 8 = 0.76 ± 0.11 and 
= 0.250 ± 0.022, again consistent with the standard model of cosmology. 

The 6dFGS measurement of fas has been combined with other datasets in sev- 



eral subsequent publications (e.g. Rapetti et al. 2012 Hudson & Turnbull 2012 



Samushia et al. 2012). 



In chapter [4] we calculate the 6dFGS projected correlation function w p {r p ) for 
four volume limited sub-samples with thresholds in redshift and stellar mass in 
order to investigate how galaxies populate dark matter halos. We found that the 
halo model fits the data better than empirical power-law fits yielding smaller values 
of x 2 , especially in the case of our high stellar mass samples for which, the power-law 
fits give a reduced \ 2 ^ 1-5- 

We find that the typical halo mass Mi, where the average dark matter halo 
has its first satellite, as well as the satellite power law index a, increase with stellar 
mass. This indicates that galaxies with higher stellar mass sit in more massive dark 
matter halos, which accumulate satellites faster with growing mass compared to 
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low stellar mass galaxies. Furthermore we found a relation between Mi and the 
minimum dark matter halo mass M m i n of Mi f=s 22M m j n , in good agreement with 
studies of SDSS galaxies Zehavi et al. 2005b, 2010 . The satellite fraction declines 
with increasing stellar mass from 21% at M ste ii a r = 2.6 x lO lo /i~ 2 M to 12% at 
M ste ii ar = 5.4 x lO lo /i~ 2 M indicating that higher stellar mass galaxies are more likely 
to be central galaxies. The effective galaxy bias increases with stellar mass and can 
be approximated by a power law given by [logi O (M ste n ar //i~ 2 M )/lO] 5 ' 74 x (0.8/cr 8 ). 



We compare these results with two mock catalogues, based on the Croton et 



al. 2006 and Bower et al. 2006 semi-analytic models. We derive these mock 



catalogues from the Millennium Simulation including the 6dFGS selection effects. 



The Croton et al. 2006 model under-predicts the satellite fraction, while the Bower 



et al. 2006 model is in better agreement with our observations. Since the satellite 



fraction and the effective dark matter halo mass, M e g- are strongly linked, the Bower 



et al. 2006 model predictions of M e g are again in better agreement with our ob- 



servations. 6dFGS allows a powerful test of semi-analytic models, because of (1) 
the robust stellar mass estimates and (2) the focus on "red satellites" , which semi- 
analytic models struggled to reproduce in the past. Our results can be used as a 
new constraint on semi-analytic models in the future. 



Finally we compare our results with Mandelbaum et al. 2006 who studied halo 
occupation as a function of stellar mass using galaxy-galaxy weak lensing. We find 
overall good agreement, which represents a valuable crosscheck of these two differ- 
ent techniques for studying matter clustering. Although our analysis depends on 
slightly more assumptions and covers a smaller range in stellar mass, the 6dFGS 
results are extremely precise compared with the lensing results. 



The BAO analysis in chapter [2] as well as the RSD analysis in chapter [3] include 
predictions for the two future galaxy surveys TAIPANQand WALLAB^Q Our anal- 
ysis shows that these two very different surveys will be able to improve upon the 
6dFGS result discussed in this thesis. The TAIPAN survey has an effective volume 
three times as big as 6dFGS and will be able to constrain the Hubble Constant to 

1 TAIPAN: Transforming Astronomical Imaging surveys through Polychromatic Analysis of Neb- 
ulae 

2 http://www. atnf.csiro.au/research/ WALLABY 
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3% precision using the BAO technique. The WALLABY survey, with its very small 
galaxy bias, is much better set-up to measure the growth of structure, where it will 
be able to constrain fa 8 to about 4 - 10% precision, depending on the modelling of 
non-linear structure formation. 

In summary, in this thesis we have presented tests of the current standard 
model of cosmology using baryon acoustic oscillations as a geometrical probe of the 
expansion of the Universe as well as redshift-space distortions to test the growth 
of structure. Since geometric probes are able to measure the dark energy equation 
of state parameter w, they can test different dark energy models, which predict 
deviations from a cosmological constant at w — — 1. The measurement of the growth 
of structure, however, is able to test the laws of Gravity. A weaker gravitational force 
on cosmic scales could be an alternative explanation for the accelerated expansion 
measured by geometrical probes. 

All our measurements are in agreement with the current standard model of 
cosmology, which treats dark energy as a cosmological constant. We discussed the 
theoretical difficulties of this model in our introduction, but unfortunately our data 
is not able to solve these problems. To find a physical explanation for dark energy 
and its behaviour as a cosmological constant remains one of the biggest challenges in 
cosmology and will call for better observational data and more innovative theoretical 
models. 



5.0.1 Future galaxy surveys 

The tremendous success of galaxy redshift surveys, in complementing the CMB to 
constrain cosmological models, sparked a large interest in the field. We already 
discussed extensively the two future galaxy surveys WALLABY and TAIPAN. To 
conclude this thesis we will give a short overview of some other projects and discuss 
their scientific goals. 

Future galaxy surveys require dedicated larger telescopes, like the ongoing 
Baryon Oscillation Spectroscopic Survey (BOSS) |Schlegel, White fc Eisenstein 



2009 1, which will run until 2014, or the Dark Energy Survey (DES) jAbbott et al. 
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2005 which commenced data taking in 2012. The main motivation for these sur- 



veys is to study the physics behind the acceleration of the expansion of the Universe. 
However, the same data will allow us to learn more about many open questions in 
astrophysics and cosmology like galaxy formation, inflation and neutrino properties. 

The Dark Energy Survey is a photometric galaxy survey using the 4m Victor 
M. Blanco Telescope located at Cerro Tololo Inter-American Observatory (CTIO) 
in Chile. The main innovation of that project consists in the development of a new 
camera called DECam. The survey will image 5000 deg 2 of the southern sky and 
will take five years to complete. 



Looking even further into the future the proposed BigBOSS experiment Schlegel 



et al. , 2011 will study Baryon Acoustic Oscillations in the redshift range 0.2 < z < 



3.5 using a dedicated 4m telescope. A new 5000-fibre spectrograph covering a 3- 
degree diameter field will detect the redshifts of about 15.3 million Emission Line 
Galaxies (ELGs), 3.4 million Luminous Red Galaxies (LRGs), and 630 000 Quasars 
(QSOs). 



The ESA space mission Euclid Laureijs, 2009| is a 1.2m telescope together with 
a CCD and a spectrograph and will map about 2/3 of the sky (15 000 deg 2 ). The 
mission is scheduled to launch in 2020, measuring the redshifts of billions of galaxies 
through the detection of the Ha emission line. 

On a similar time scale there is the Large Synoptic Survey Telescope (LSST), 
which is a planned wide-field "survey" telescope, capable of covering 20 000 deg 2 
of the southern sky using an 8.4m ground-based telescope. Each patch of sky will 
be visited about 1000 times in ten years detecting roughly 250 000 supernovae per 
year. First light is scheduled for 2020 and full operations for a ten-year survey will 
commence in January 2022. 

Finally, the Square Kilometre Array (SKA) will use thousands of radio tele- 
scopes covering a collecting area of about one square kilometre to detect the signal 



of hydrogen in the Universe, up to the epoch of re-ionisation Blake et al. , 2004 



Radio galaxy surveys have the advantage that it is relatively simple to build inter- 
ferometers, which means that the physical limitations of building one big structure 
can be bypassed, by building many small telescopes and linking them together. 
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Galaxy redshift surveys are amongst the most successful projects in cosmol- 
ogy/astronomy. The Sloan Digital Sky Survey (SDSS) lead to currently ~ 5000 
refereed papers and ~ 200 000 citations (Michael Strauss, private communication, 
see also Madrid & Macchetto 2009 Savaglio & Grothkopf 2012). For comparison, 
WMAP lead to ~ 1700 refereed papers and 100 000 citations. Galaxy redshift sur- 
veys are fantastically equipped to address the big open questions in astronomy and 
cosmology. With the future survey projects outlined above the coming years will be 
very exciting for cosmology and who knows what we will find. 
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A. GENERATING LOG-NORMAL MOCK CATALOGUES 



A Generating log-normal mock catalogues 

Here we explain in detail the different steps used to derive a log-normal mock cat- 
alog 116, cLS £L useful guide for researchers in the field. We start with an input power 



spectrum, (which is determined as explained in Section 2.3.3) in units of h Mpc . 
We set up a 3D grid with the dimensions L x X Ly X L z = 1000 x 1000 x 1000/i- 1 Mpc 
with 200 3 sub-cells. We then distribute the quantity P(k)/V over this grid, where V 
is the volume of the grid and k = y / fc|~+~fc|~+~fcf with k x = n x 2ir/L x and n x being 
an integer value specifying the x coordinates of the grid cells. 

Performing a complex-to-real Fourier transform (FT) of this grid will produce a 
3D correlation function. Since the power spectrum has the property P(—k) = P(k)* 
the result will be real. 

The next step is to replace the correlation function £ (r) at each point in the 3D 
grid by ln[l + £(r)], where In is the natural logarithm. This step prepares the input 
model for the inverse step, which we later use to produce the log-normal density 
field. 

Using a real-to-complex FT we can revert to fc-space where we now have a 
modified power spectrum, P\ n (k). At this point we divide by the number of sub- 
cells N c . The precise normalisation depends on the definition of the discrete Fourier 



transform. We use the FFTW library Frigo & Johnson, 2005 , where the discrete 
FT is defined as 

N c -1 

Yi=J2 X i ex P [±2nijV=l/N c ] . (5.1) 

3=0 

The modified power spectrum P\ n (k) is not guarantied to be neither positive defined 



nor a real function, which contradicts the definition of a power spectrum. |Weinberg 
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& Cole 1992 suggested to construct a well defined power spectrum from P\ n (k) by 

P{ D (k) = max [(D, Re[P ln (£)]l . (5.2) 

— * 

We now generate a real and an imaginary Fourier amplitude 5(k) for each point on 



the grid by randomly sampling from a Gaussian distribution with r.m.s. y P{ n (k)/2. 
However, to ensure that the final over-density field is real, we have to manipulate 
the grid, so that all sub-cells follow the condition S(— k) = 5(k)*. 

Performing another FT results in an over-density field S(x) from which we cal- 
culate the variance o'q. The mean of 5(x) should be zero. The log-normal density 
field is then given by 

fi L (x) = exp [5(f) - ctq/2] , (5.3) 

which is now a quantity defined on [0, oof only, while 5(x) is defined on ] — oo, oof. 

Since we want to calculate a mock catalogue for a particular survey we have to 
incorporate the survey selection function. If W(x) is the selection function with the 
normalisation ^2W(x) = 1, we calculate the mean number of galaxies in each grid 
cell as 

n g (x) = N W(x) /i L (x), (5.4) 

where N is the total number of galaxies in our sample. The galaxy catalogue itself 
is than generated by Poisson sampling n g (x). 

The galaxy position is not defined within the sub-cell, and we place the galaxy 
in a random position within the box. This means that the correlation function 
calculated from such a distribution is smooth at scales smaller than the sub-cell. It 
is therefore important to make sure that the grid cells are smaller than the size of 
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B. COMPARISON OF LOG-NORMAL AND JACK-KNIFE ERROR 
ESTIMATES 
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Figure 5.1: Correlation function error for different values of Pq. The weighting with 
Po = 40 000/;, 3 Mpc~ 3 reduces the error at the BAO scale by almost a factor of 
four compared to the case without weighting. The red dashed line indicates the 
jack-knife error. 



the bins in the correlation function calculation. In the 6dFGS calculations presented 
in this paper the grid cells have a size of 5/i -1 Mpc, while the correlation function 
bins are 10/i -1 Mpc in size. 



B Comparison of log-normal and jack-knife error 
estimates 

We have also estimated jack-knife errors for the correlation function, by way of 
comparison. We divided the survey into 18 regions and calculated the correlation 
function by excluding one region at a time. We found that the size of the error- 
bars around the BAO peak varies by around 20% in some bins, when we increase 
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Figure 5.2: Correlation matrix of the jack-knife errors (upper left triangle) and 
log- normal errors (lower right triangle). 



the number of jack-knife regions from 18 to 32. Furthermore the covariance matrix 
derived from jack-knife resampling is very noisy and hard to invert. 



We show the jack-knife errors in Figure 5.1 The jack-knife error shows more 



noise and is larger in most bins compared to the log-normal error. The error shown 



in Figure 5^ is only the diagonal term of the covariance matrix and does not include 
any correlation between bins. 

The full error matrix is shown in Figure |5.2[ where we plot the correlation 
matrix of the jack-knife error estimate compared to the log- normal error. The jack- 
knife correlation matrix looks much more noisy and seems to have less correlation 
in neighbouring bins. 
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C. WIDE-ANGLE FORMALISM 



The number of jack-knife regions can not be chosen arbitrarily. Each jack-knife 
region must be at least as big as the maximum scale under investigation. Since 
we want to test scales up to almost 200/i _1 Mpc our jack-knife regions must be 
very large. On the other hand we need at least as many jack-knife regions as we 
have bins in our correlation function, otherwise the covariance matrix is singular. 
These requirements can contradict each other, especially if large scales are analysed. 
Furthermore the small number of jack-knife regions is the main source of noise (for 



a more detailed study of jack-knife errors see e.g. Norberg et al. 2008). 

Given these limitations in the jack-knife error approach, correlation function 
studies on large scales usually employ simulations or log-normal realisations to derive 
the covariance matrix. We decided to use the log-normal error in our analysis. We 
showed that the jack-knife errors tend to be larger than the log-normal error at 
larger scales and carry less correlation. These differences might be connected to 
the much higher noise level in the jack-knife errors, which is clearly visible in all 
our data. It could be, however, that our jack-knife regions are too small to deliver 
reliable errors on large scales. We use the minimum number of jack-knife regions 
to make the covariance matrix non-singular (the correlation function is measured 
in 18 bins). The mean distance of the jack-knife regions to each other is about 
200/i^ 1 Mpc at the mean redshift of the survey, but smaller at low redshift. 



C Wide-angle formalism 

The general redshift space correlation function (ignoring the plane parallel approx- 
imation) depends on <j), 9 and s. Here, s is the separation between the galaxy pair, 
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s [h" 1 Mpc] 



Figure 5.3: The half opening angle 9 as a function of separation s of the 6dFGS 
weighted catalogue. The plane parallel approximation assumes 9 = 0. The mean 
half opening angle at the BAO scale is < 10°. The colour bar gives the number of 
pairs in each bin. 



9 is the half opening angle, and is the angle of s to the line of sight (see Figure 1 



in 



Raccanelli et al. 2010). For the following calculations it must be considered that 



in this parametrisation, and 9 are not independent. 

The total correlation function model, including 0(9 2 ) correction terms, is then 



given by Papai & Szapudi 2008 



£(0, 9, s) = a 00 + 2a 02 cos(20) + a 22 cos(20) + b 22 sin 2 (20) 
+ - 4a 02 cos(20) - 4a 22 - 46 22 - 4a i0 cot 2 (0) 
+ 4an cot 2 (0) - 4a i2 cot 2 (0) cos(20) + 46 u 
-8&U cos 2 {<j))] 9 2 + 0{9 4 ) 



(5.5) 
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C. WIDE-ANGLE FORMALISM 



This equation reduces to the plane parallel approximation if 9 = 0. The factors a xy 
and b xy in this equation are given by 



o-oo — 



, 2P 2p< 
1 + — + — - 
3 15 

« 2 ^ 



£ + 



Cf02 — — 



'0 



21 

3/3 : 



C 2 2 (r) + %${r) 



+ 



14 



£(0 + 2^4 2 (0 



«22 - ^ 

ft 

022 — — Qn 



140 
28 v 

M - ^ 2 2 (r) + ^£ 4 2 (r) 



/3 2 



35-^ 



aio = 



4/3 2 
an = — 



2/J + 4 f ! ] ^''M-^W 
|? KM - 2&°(r)] 



(5.6) 



where /3 = Q m (z)°' 5A5 /b, with b being the linear bias. The correlation function 
moments are given by 



2tt 2 



dk k m P hn (k)ji(rk) 



(5.7) 



with ji(x) being the spherical Bessel function of order /. 

The final spherically averaged correlation function is given by 



£(s)= ^ r / \(<f>,9,s)N(cf>,9,s)d9dcf>, 
Jo Jo 



(5.8) 



202 



where the function N(<p, 8, s) is obtained from the data. N(4>, 8, s) counts the number 
of galaxy pairs at different 0, 8 and s and includes the areal weighting sin(0) which 
usually has to be included in an integral over 0. It is normalised such that 

it 

N{<f),8,s) d8d(j) = l. (5.9) 



o Jo 



If the angle 8 is of order 1 rad, higher order terms become dominant and eq. 5.5 



is 



no longer sufficient. Our weighted sample has only small values of 8, but growing 



with s (see figure 5.3). In our case the correction terms contribute only mildly at 



the BAO scale (red line in figure 5.4). However these corrections behave like a scale 



dependent bias and hence can introduce systematic errors if not modelled correctly. 



D Partial analytical solution for the correlation 
function moments integral 



The double integrals in eqs. 3.35 and 3.38 are difficult to solve numerically. Here we 
show an analytical solution for the integral over /i which allows for a faster numerical 

as 



solution of the full integral. First we re-write eq. 3.38 



00 r 1 k^dhdy (W 

J 1 (27T) 2 

x cos(krfi)Q xy (k)Vi(fi) 
00 k m dk 



Re 



Q xy {k) I dfx e lkr ^ k ^ Vt(fi) 
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D. PARTIAL ANALYTICAL SOLUTION FOR THE CORRELATION 
FUNCTION MOMENTS INTEGRAL 




s [h" 1 Mpc] 



Figure 5.4: The black line represents the plain correlation function without redshift 
space distortions (RSD), £(r), obtained by a Hankel transform of our fiducial ACDM 
power spectrum. The blue line includes the linear model for redshift space distortions 
(linear Kaiser factor) using /3 = 0.27. The red line uses the same value of (3 but 
includes all correction terms outlined in eq. 5.5 using the N(<j),6,s) distribution of 
the weighted 6dFGS sample employed in this analysis. 



where i is the complex number. More generally the integral over /i can be written 

as 







n ~i 



d(kr) j j_ x 
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(5.10) 



The integral on the right can be solved analytically. If we set a = kr and b = (ka v ) 2 
we obtain 



djji e 



Vb 



^ . a - 2ib\ 2iaj ^ ( a + 2ib 
Daw ( — - e Daw — 

2Vb J V 2Vb 



(5.11) 



Taking the n-th derivatives of the real part of the term above gives F n (fi). The 
Dawson integral Daw(x) can be calculated using the imaginary error function erfi(x): 



Daw(x) = e x / e y dy 
Jo 



(5.12) 



'erfi(a;). 



We can than construct eq. 3.38 for the different correlation function moments: 



^b=A,xyi r ) 



k m dk 



Qxy{k)F Q (jj) 



x (2?r)2 

k "' <lk Q xy {k)\[ZF 2 ^)-FM) 



00 1 

' " Q X y(k)^m^)-3F 1 ^)\ 



(5.13) 



(2tt)2 

00 k m dk 1 

— Q^fc)- [35F 4 (/i) - 30F 2 (ai) + 3F (/i)] 



and so on. 
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E. NORMALISATION CONSTANT IN Mandelbaum et al. 2006 



E Normalisation constant in Mandelbaum et al. 
2006 



We can determine the normalisation constant k of eq. 4.47 by matching to the 



satellite fraction /™ a , which in terms of the mass function dn(M)/dM is given by 



with 



which leads to 



IM^1[N C (M) + N S 



fma _ JO dM s ^ ' (Z\'\A\ 

(5.15) 



' f™dM d ^±{N c (M) + N s (M)} 
kH 

1 [m dM^f + kH 

J Af m i n dM 



3M mill /„_,_ dM 



H =^j— I dM^^M 2 + 

'min JAf min 



dn(M) %r 
dM ; M, 

3M • rfM 



" J «m)r 



dn(M) 
dM 



(5.16) 



fc = ; JMmi - " v , (5.i7) 



where k has units of l//i _1 M and is constant for a given set of /™ a and M n 
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